• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.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.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.10a
• /
• pp.223-226
• /
• 1995

A three-dimensional elastic-plastic finite element analysis using the explicit time integration method has been performed for the characterization of theimpact forming machines. The block upsetting using a forging hammer has been analyzed. The effects of machine type, work capacity of equipment and the mass ratio in an anvil-type hammer have been studied through the analysis.

• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.10-18
• /
• 2016

Numerical boundary conditions are as important as the governing equations when analyzing the fluid flows numerically. An explicit boundary condition method updates the solutions at the boundaries with extrapolation from the interior of the computational domain, while the implicit boundary condition method in conjunction with an implicit time integration method solves the solutions of the entire computational domain including the boundaries simultaneously. The implicit boundary condition method, therefore, is more robust than the explicit boundary condition method. In this paper, steady compressible 2-Dimensional Navier-Stokes solver is developed. We present the implicit boundary condition method coupled with LU-SGS(Lower Upper Symmetric Gauss Seidel) method. Also, the explicit boundary condition method is implemented for comparison. The preconditioning Navier-Stokes equations are solved on unstructured meshes. The numerical computations for a number of flows show that the implicit boundary condition method can give accurate solutions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.8
• /
• pp.33-40
• /
• 2006

A parallel computing strategy for finite element(FE) processing is described and implemented in nonlinear explicit FE code and its parallel performances are evaluated. A self-made linux-cluster supercomputer with 520 CPUs is used as a bench mark test bed. It is observed that speed-up is increased almost idealy even up to 256 CPUs for a large scale model. A communication over head and its effect on the parallel performance is also examined. Parallel performance is compare with the commercial code and developed code shows superior performance as the number of CPUs used are increased.

• Structural Engineering and Mechanics
• /
• v.62 no.2
• /
• pp.151-162
• /
• 2017

Two $Chang-{\alpha}$ dissipative family methods and two $KR-{\alpha}$ family methods were developed for time integration recently. Although the four family methods are in the category of the dissipative structure-dependent integration methods, their performances may be drastically different due to the detrimental property of weak instability or overshoot for the two $KR-{\alpha}$ family methods. This weak instability or overshoot will result in an adverse overshooting behavior or even numerical instability. In general, the four family methods can possess very similar numerical properties, such as unconditional stability, second-order accuracy, explicit formulation and controllable numerical damping. However, the two $KR-{\alpha}$ family methods are found to possess a weak instability property or overshoot in the high frequency responses to any nonzero initial conditions and thus this property will hinder them from practical applications. Whereas, the two $Chang-{\alpha}$ dissipative family methods have no such an adverse property. As a result, the performances of the two $Chang-{\alpha}$ dissipative family methods are much better than for the two $KR-{\alpha}$ family methods. Analytical assessments of all the four family methods are conducted in this work and numerical examples are used to confirm the analytical predictions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2245-2252
• /
• 2002

Design sensitivity analysis scheme is proposed in an elasto -plastic finite element method with explicit time integration using a direct differentiation method. The direct differentiation is concerned with large deformation, the elasto-plastic constitutive relation, shell elements with reduced integration and the contact scheme. The design sensitivities with respect to the process parameter are calculated with the direct analytical differentiation of the governing equation. The sensitivity results obtained from the present theory are compared with that obtained by the finite difference method in a class of sheet metal forming problems such as hemi-spherical stretching and cylindrical cup deep-drawing. The result shows good agreement with the finite difference method and demonstrates that the preposed sensitivity calculation scheme is a pplicable in the complicated sheet metal forming analysis and design.

• Kyungpook Mathematical Journal
• /
• v.56 no.2
• /
• pp.311-327
• /
• 2016

In this paper, we propose an error embedded Runge-Kutta method to improve the traditional embedded Runge-Kutta method. The proposed scheme can be applied into most explicit embedded Runge-Kutta methods. At each integration step, the proposed method is comprised of two equations for the solution and the error, respectively. These solution and error are obtained by solving an initial value problem whose solution has the information of the error at each integration step. The constructed algorithm controls both the error and the time step size simultaneously and possesses a good performance in the computational cost compared to the original method. For the assessment of the effectiveness, the van der Pol equation and another one having a difficulty for the global error control are numerically solved. Finally, a two-body Kepler problem is also used to assess the efficiency of the proposed algorithm.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• 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.

• Composites Research
• /
• v.18 no.5
• /
• pp.15-20
• /
• 2005

In this paper, damage induced acoustic emission in the composite plate in numerically simulated by using the three dimensional finite element method and explicit time integration. Acoustic source is modeled by equivalent volume source. To verify the proposed method, dynamic displacements due to the elastic wave are compared with the experiment when the fiber is broken in the single fiber embedded isotropic plate. For the laminated composite plates, the results are compared between homogenized model and DNS approach which models fibers and matrix separately. To capture high frequencies in the elastic wave, small time step size and a large number of meshes are required. The parallel computing technology is introduced to solve a large scale problem efficiently.