• 대한기계학회논문집A
• /
• 제22권1호
• /
• pp.111-120
• /
• 1998

This paper is the first one of two-parted research efforts focusing on the modeling of rubber pad forming process. The rubber pad, driven by the pressurized fluid during the forming process, pushes the sheet metal to solid tool half and forms a part to final shape. In this part of the paper, a numerical procedure for the FE analysis of the rubber pad deformation is presented. The developed three-dimensional FE model is based on the total Lagrangian description of rubber maerial characterized by nearly incompressible hyper-elastic behavior under a large deformation assumption. Validity of the model as well as effects of different algorithms corresponding to incompresibility constraints and time integration methods on numerical solution responses are also demonstrated.

• 대한의용생체공학회:의공학회지
• /
• 제5권1호
• /
• pp.9-14
• /
• 1984

The asymptotic solution using the Tailor series has been given explicit form for the solute concentration and overall solute removal in hemodiafilter using one dimensional model. The numerical solutions have been calculated within 0.001% error by the Romberg integration method. Compared with the numerical solutions, the oneterm asymptotic solutions were found to be within 3% error for the condition > 3.0 and three-terms asymtotic solutions were required for the condition >0.7 where denotes measure of convection over diffusional transport and a the ratio of blood flow rate over dialysate flow rate.

• 대한기계학회논문집A
• /
• 제27권10호
• /
• pp.1793-1799
• /
• 2003

An optimal design of a multi-layered plate structure to endure high-velocity impact has been suggested by using size optimization after numerical simulations. The NET2D, a Lagrangian explicit time-integration finite element code for analyzing high-velocity impact, was used to find the parameters for the optimization. Three different materials such as mild steel, aluminum for a multi-layered plate structure and die steel for the pellet, were assumed. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, Johnson-Cook model and Phenomenological Material Model were used as constitutive models for the simulation. It was carried out with several different gaps and thickness of layers to figure out the trend in terms of those parameters' changes under the constraint, which is against complete penetration. Also, the measuring domain has been shrunk with several elements to reduce the analyzing time. The response surface method based on the design of experiments was used as optimization algorithms. The optimized thickness of each layer in which perforation does not occur has been obtained at a constant velocity and a designated total thickness. The result is quite acceptable satisfying both the minimized deformation energy and the weight criteria. Furthermore, a conceptual idea for topology optimization was suggested for the future work.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2005년도 춘계학술발표대회 논문집
• /
• pp.172-175
• /
• 2005

The numerical simulation is performed for the acoustic emission and the wave propagation due to fiber breakage in single fiber composite plates by the finite element transient analysis. The acoustic emission and the following wave motions from a fiber breakage under a static loading is simulated to investigate the applicability of the explicit finite element method and the equivalent volume force model as a simulation tool of wave propagation and a modeling technique of an acoustic emission. For such a simple case of the damage event under static loading, various parameters affecting the wave motion are investigated for reliable simulations of the impact damage event. The high velocity and the small wave length of the acoustic emission require a refined analysis with dense distribution of the finite element and a small time step. In order to fulfill the requirement for capturing the exact wave propagation and to cover the 3-D simulation, we utilize the parallel FE transient analysis code and the parallel computing technology.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.138-144
• /
• 1998

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

• 한국재난정보학회 논문집
• /
• 제10권1호
• /
• pp.123-140
• /
• 2014

회전쉘의 정적문제 및 동적문제에 대해서 선점법에 의한 공간영역의 이산화에 대한 유효성을 비교 검토하였고, 정식화에 있어서 이산화 후의 연립방정식 및 시간에 관한 연립상미분방정식의 계수행렬에 대한 영향을 검토하였다. 수치해석법에는 동적특성의 문제를 간단하고 효과적인 방법을 위하여 선점법을 제안하였으며, 수치해석예로서 평균2승잔차법과의 비교검토로 고정도의 해석을 위한 정확성을 검증하였다.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2003년도 춘계학술대회논문집
• /
• pp.258-261
• /
• 2003

In this study, finite element analysis for multi-stage deep drawing process of rectangular configuration with extreme aspect ratio is carried out especially for the blank design. The analysis of rectangular deep drawing process with extreme aspect ratio is likewise very difficult with respect to the design process parameters including the intermediate die profile. In order to solve the difficulties, numerical approach using finite element method is performed in the present analysis and design. A series of experiments for multi-stage rectangular deep drawing process are conducted and the deformed configuration is investigated by comparing with the results of the finite element analysis. Additionally, to minimize amount of removal material after trimming process, finite element simulation is applied for the blank modification. The analysis incorporates brick elements for a rigid-plastic finite element method with an explicit time integration scheme using LS-DYNA3D.

• Interaction and multiscale mechanics
• /
• 제2권2호
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.

• International Journal of Automotive Technology
• /
• 제1권1호
• /
• pp.35-41
• /
• 2000

The crashworthiness of vehicles with finite element methods depends on the geometry modeling and the material properties. The vehicle body structures are generally composed of various members such as frames, stamped panels and deep-drawn parts from sheet metals. In order to ensure the impact characteristics of auto-body structures, the dynamic behavior of sheet metals must be examined to provide the appropriate constitutive relation. In this paper, high strain-rate tensile tests have been carried out with a tension type split Hopkinson bar apparatus specially designed for sheet metals. Experimental results from both static and dynamic tests with the tension split Hopkinson bar apparatus are interpolated to construct the Johnson-Cook and a modified Johnson-Cook equation as the constitutive relation, that should be applied to simulation of the dynamic behavior of auto-body structures. Simulation of auto-body structures has been carried out with an elasto-plastic finite element method with explicit time integration. The stress integration scheme with the plastic predictor-elastic corrector method is adopted in order to accurately keep track of the stress-strain relation for the rate-dependent model accurately. The crashworthiness of the structure with quasi-static constitutive relation is compared to the one with the rate-dependent constitutive model. Numerical simulation has been carried out for frontal frames and a hood of an automobile. Deformed shapes and the Impact energy absorption of the structure are investigated with the variation of the strain rate.

• Structural Engineering and Mechanics
• /
• 제47권5호
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.