• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.245-255
• /
• 1994

A two-dimensional Lagrangian finite-difference computer program is developed for the wave propagation analysis of impact phenomena. The numerical scheme is the standard method originally proposed by Von Neuman and Richtmyer, using artificial viscosity to smooth shock fronts. The material model used in the study is the standard hydrodynamic-elastic-plastic relations with Von-Mises yield criterion. A test configuration consisted of a target and a projectile were calculated to understand the response of a colliding event. However, the computer code is in plane strain, the calculations were intended for generating the qualitative features of the model behaviors. Nevertheless, the computational results were consistent with the experimental observations and provided a rational basis to interpret the modes of failures.

• Geomechanics and Engineering
• /
• v.21 no.5
• /
• pp.489-501
• /
• 2020

This paper conducts a complicated coupled effective stress analysis of X-section-in-place concrete (XCC) pile installation and consolidation processes using the dual-stage Eulerian-Lagrangian (DSEL) technique incorporating the modified Cam-clay model. The numerical model is verified by centrifuge data and field test results. The main objective of this study is to investigate the shape effect of XCC pile cross-section on radial total stress, excess pore pressure and time-dependent strength. The discrepancies of the penetration mechanism and set-up effects on pile shaft resistance between the XCC pile and circular pile are discussed. Particular attention is placed on the time-dependent strength around the XCC pile shaft. The results show that soil strength improved more significantly close to the flat side compared with the concave side. Additionally, the computed ultimate shaft resistance of XCC pile incorporating set-up effects is 1.45 times that of the circular pile. The present findings are likely helpful in facilitating the incorporation of set-up effects into XCC pile design practices.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.1-10
• /
• 1991

Two-dimensional explicit finite element code was developed. The transient dynamics code can analyse large deformations of non-linear materials subjected to extremely high strain rates. The Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion, thus the stiffness matrix is not introduced. Cylinder upsetting and ring compression problems are simulated to check the effects of friction and inertia force. It is shown that (1) calculated results agree very well with experimental results, (2) constant shear friction method overestimates the decrease of inner ring radius and then underestimates after on in comparison with the Coulomb friction method, and (3) the effect of the increase in initial strain rate is similar to the effect of higher frictional coefficient.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2113-2122
• /
• 1994

The 3-D FEM analysis for simulating the stamping operation of planar anisotropic sheet metals with arbitrarily-shaped tools is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The consistent full set of governing relations, comprising equilibrium equation and mesh-normal geometric constraints, is appropriately linearized. The linear triangular elements are used for depicting the formed sheet, based on membrane approximation. Barlat's non-quadratic anisotropic yield criterion(strain-rate potential) is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and non-quadratic function parameter. The planar anisotropic finite element formulation is tested with the numerical simulations of the stamping of an automotive hood inner panel and the drawing of a hemispherical punch. The in-plane anisotropic effects on the formability of both mild steel and aluminum alloy sheet metals are examined.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• Journal of Korean Society for Atmospheric Environment
• /
• v.23 no.4
• /
• pp.478-486
• /
• 2007

Numerical analysis had been performed to understand flow characteristics of the flue gas and slurry in the absorber of a flue gas desulphurization (FGD) system using computational fluid dynamics (CFD) technique. Two-fluid(Euler-Lagrangian) model had been employed to simulate physical phenomenon, which slurry particles injected through slurry spray nozzles fall down and bump into the flue gas inflowing through inlet duct. It was not necessary to adopt pre-defined pressure drop inside the absorber because interaction between flue gas and slurry particles was considered. Hundreds of slurry spray nozzles were considered with the spray velocity at the nozzles, swirl velocity and spreading angle. The results note that the flow disturbance of flue gas is found at the bottom of the absorber, and the current rising with high speed stream is observed in the opposite region of the inflow duct. The high speed stream is reduced as the flue gas goes up, because the high speed stream of flue gas dumps falling slurry particles due to momentum exchange between flue gas and slurry particles. In spite of some disproportion in slurry distribution inside the absorber, escape of slurry particles from the absorber facility is not observed. The pressure drop inside the absorber is mainly occurred at the bottom section.

• The Korean Journal of Rheology
• /
• v.8 no.3_4
• /
• pp.215-225
• /
• 1996

평판형의 열가소성 수지를 유리전이온도 이상으로 가열한 다음 압력을 가함으로써 원하는 형상의 제품을 성형하는 열성형공정은 대상 수지가 큰 변형을 일으킬 뿐만 아니라 비선형적 거동을 보이게 된다. 따라서 수지의 변형거동 예측과 최적성형조건의 설정에 많은 어려움과 시행착오를 거치게 되는 바, 열성형 공정의 최적화를 위한 연구의 일환으로 원형 평판위 수지를 대상으로 수지의 부풀림 거동과 이에 따른 두께 분포를 예측할수 있는유한요 소법의 수치모사 알고리듬을 개발하고자 하였다. Piola-Kirchhoff 응력 텐서와 Green 변형 텐서 및 lagrangian 변형 텐서를 사용하여 평판상의 응력-변형에 대한 비선형의 에너지 수 지식을 수립하고 Newton-Raphson 반복수렴법을 이용하여 근사적으로 해석하였으며 수지의 유변학적 구성방정식으로는 neo-Hookean 모델, Mooney-Rivlin 모델 및 Ogden 모델등의 초탄성 모델을 사용하여 그결과를 비교하였다. 수치모사에는 두께가 매우 얇기 때문에 두께 방향의 응력변화를 무시할수 있는 membrane 가정을 도입한 2차원적 해석과 두께 방향의 응력 변화를 고려하는 3차원적 해석을 모두 수행하고 그 차이를 비교하였으며 3차원적 해석 의 경우에는 penalty법을 이용하여 비 압축성을 만족하였다. 일차적으로 내압을 받는 두꺼 운 원통계에 대한 수치모사 해석을 수행하고 완전해와 비교함으로써 개발된 수치모사 알고 리듬의 열성형 공정에의 적용 타당성을 검증하였으며 이를 이용하여 원형 평판의 자유부풀 림거동을 예측한 결과 Treloar 등의 실험결과와 잘 부합함을 확인하였다. 또 간단한 형상의 금형이 있는 경우와 반지름 방향으로의 온도변화에 따른 수지의 변형거동을 해석함으로써 실제 열성형 공정에서 요구되고 있는 성형품의 두께 분포를 균일하게 하기 위한 방안을 제 시하였다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.5
• /
• pp.425-434
• /
• 2007

An integrated adsorption/catalytic process has been considered to treat dioxin and $NO_x$ simultaneously. The process consists of a cyclone and a reactor with nine bag filters. In this study, numerical analysis has been performed to understand flow characteristics with inflow-duct types in the reactor. To consider flue gas and activated carbon particles simultaneously, Euler-Lagrangian model was employed. Fundamental flow patterns of flue gas and activated carbon particles, pressure distribution and distribution of activated carbon have been obtained from the numerical analysis. Also trace length and residence time of flue gas, residence time of activated carbon particles have been calculated directly. Flow patterns of flue gas and activated carbon particles in the reactor were very complicated and they moved along very various paths. Therefore, their residence time in the reactor was also various. The flow characteristics in the reactor were strongly influenced by inflow-duct types. The results obtained would be effectively used to estimate the removal efficiency in the reactor once the residence time is combined with the reaction equation.

• Computers and Concrete
• /
• v.25 no.3
• /
• pp.255-266
• /
• 2020

This study presents comparation of fixed and viscos boundary condition effects on three-dimensional earthquake response and performance of a RCC dam considering linear and non-linear response. For this purpose, Cine RCC dam constructed in Aydın, Turkey, is selected in applications. The Drucker-Prager material model is considered for concrete and foundation rock in the nonlinear time-history analyses. Besides, hydrodynamic effect was considered in linear and non-linear dynamic analyses for both conditions. The hydrodynamic pressure of the reservoir water is modeled with the fluid finite elements based on the Lagrangian approach. The contact-target element pairs were used to model the dam-foundation-reservoir interaction system. The interface between dam and foundation is modeled with welded contact for both fixed and viscos boundary conditions. The displacements and principle stress components obtained from the linear and non-linear analyses are compared each other for empty and full reservoir cases. Seismic performance analyses considering demand-capacity ratio criteria were also performed for each case. According to numerical analyses, the total displacements and besides seismic performance of the dam increase by the effect of the viscous boundary conditions. Besides, hydrodynamic pressure obviously decreases the performance of the dam.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.