• 대한기계학회논문집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.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1995년도 추계학술대회 논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제21권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.

• 한국항공우주학회지
• /
• 제34권8호
• /
• pp.33-40
• /
• 2006

비선형 외연 유한요소법에서 유한요소 병렬 처리 방안을 기술하고 코드에 구현하였다. 성능테스트 장비로 자체 구축한 520 개의 CPU를 갖는 리눅스 클러스터 슈퍼컴퓨터를 사용하였다. 대규모 모델 테스트 결과 256 개의 CPU 까지도 거의 이상적인 속도 증가를 보였다. 유한요소 계산시간 대비 통신시간 계산이 전체 성능에 미치는 영향도 검토하였다. 사용 프로세서가 증가할수록 상용코드의 병렬 성능 대비 더 좋은 성능을 보이는 것으로 나타났다.

• Structural Engineering and Mechanics
• /
• 제62권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.

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

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

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 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.

• Composites Research
• /
• 제18권5호
• /
• pp.15-20
• /
• 2005

본 논문에서는 복합재 파손에 의한 음향방출해석을 3차원 유한요소법과 외연시간적분법을 이용하여 구현하였다. 음원모델은 등가체적력 모델을 사용하였다. 계산기법의 타당성을 검증하기 위해 단일 섬유가 내재된 등방성 평판에서 섬유파손 시 발생하는 탄성파에 의한 동적변위를 시험과 비교하였다. 적층 복합재의 경우, 섬유와 기지를 각기 모델링한 방법과 균질화한 모델을 비교하여 차이점을 비교하였다. 음향방출에서 발생하는 고주파 성분을 검출하기 위해 계산시간 스텝이 매우 작아야 하며, 매우 많은 자유도의 모델이 동반되어야 한다. 이러한 대규모 문제를 효과적으로 해결하기 위해 병렬 계산 기법을 도입하였다.