• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.589-600
• /
• 1996

An explicit 4th order time integration scheme for solving the convection-diffusion equation is discussed in this paper. A system of ordinary differential equations are derived first by discretizing the spatial derivatives of the relevant PDE using the finite difference method. The integration of the ODEs is then carried out using a 4th order scheme and a self-adaptive technique based on the spatial grid spacing. For a non-uniform spatial grid, different time step sizes are used for the integration of the ODEs defined at different spatial points, which improves the computational efficiency significantly. A numerical example is also discussed in the paper to demonstrate the implementation and effectiveness of the method.

• Structural Engineering and Mechanics
• /
• v.18 no.4
• /
• pp.411-428
• /
• 2004

In this study, a newly developed unconditionally stable explicit method is employed to solve momentum equations of motion in performing pseudodynamic tests. Due to the explicitness of each time step this pseudodynamic algorithm can be explicitly implemented, and thus its implementation is simple when compared to an implicit pseudodynamic algorithm. In addition, the unconditional stability might be the most promising property of this algorithm in performing pseudodynamic tests. Furthermore, it can have the improved properties if using momentum equations of motion instead of force equations of motion for the step-by-step integration. These characteristics are thoroughly verified analytically and/or numerically. In addition, actual pseudodynamic tests are performed to confirm the superiority of this pseudodynamic algorithm.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.03a
• /
• pp.188-191
• /
• 1997

In the recent sheet metal forming simulations, it increases to adopt the dynamic explicit method for an effective computation and the elastic-plastic formulation for stress recovery. It is inevitable in the dynamic explicit method that some noises occur and sometimes partly spoil results of simulations. It is severer when complicated contact conditions are included in simulations. An effective method to control these noises is introduction of damping effects. In this paper, the concept of contact damping is introduced in order to suppress noises due to complicated contact conditions. This is checked by analyzing a simple sheet metal forming process(U-bending). From the computational results, it is shown that the contact damping can effectively control the noises due to contacts and develop more reliable internal stress states.

• Computational Structural Engineering
• /
• v.8 no.3
• /
• pp.5-10
• /
• 1995

공학분야에서 흔히 접하게 되는 접촉충격 문제들의 해석을 위해 접촉충격의 지배방정식과 효율적인 해석방법인 pinball contact-impact 알고리즘에 대해 기술하였다. 예제로 보여준 충돌하는 두 구조계의 접촉충격 해석과 자동차의 crashworthiness 해석결과들은 전술한 알고리즘이 매우 효율적임을 입증하고 있다. 자동차의 crashworthiness 해석과 같이 해석대상 모델이 크고 복잡한 경우에는 짧은 경과시간 동안의 거동변화를 파악하고자 하더라도 엄청난 계산시간이 소요된다. 따라서 효율적인 crashworthiness 해석을 위해서는 전술한 접촉충격 알고리즘 외에도, 논외의 대상이라 본문에서 소개하지는 않았으나, 강력한 요소(one-point integration element), 효과적인 비선형 동적해석 알고리즘, 해의 정확도를 높일 수 있는 adaptive method, 효율적인 시간적분을 위한 mixed time integration method(subcycling)와 explicit method 등을 복합적으로 혼합하여 사용하여야 한다.

• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.

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

In this paper, a program for real time simulation of a vehicle is developed. The program uses relative coordinates and BEF(Backward Difference Formula) numerical integration method. Numerical tests showed that the proposed implicit method is more stable in carring out the numerical integration for vehicl dynamics than the explicit method. Hardware requirements for real time simulation are suggested. Algorithms of parallel processing is developed with DSP (digital signal processor).

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2000.10a
• /
• pp.136-140
• /
• 2000

It is desirable to predict springback quantitatively and accurately for the tool and process design in sheet stamping operations, however, it is blown very difficult. The result of springback analysis by the finite element method is sensitively influenced by numerical factors such as blank element size, number of integration point, punch velocity, contact algorithm etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on springback Quantitatively and to obtain the combination of numerical factors which yields the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by explicit FEM says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The 2D draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model.

• Bulletin of the Korean Chemical Society
• /
• v.27 no.12
• /
• pp.1976-1980
• /
• 2006

We have calculated the reaction probability and the reaction cross-section of the $N(^4S)+O_2(X^3\sum_{g}^{-})\;\rightarrow\;NO(X^2\Pi)+O(^3P)$ reaction by the quasiclassical trajectory method with the 6th-order explicit symplectic algorithm, based on a new ground potential energy surface. The advantage of the 6th-order explicit symplectic algorithm, conserving both the total energy and the total angular momentum of the reaction system during the numerical integration of canonical equations, has firstly analyzed in this work, which make the calculation of the reaction probability more reliable. The variation of the reaction probability with the impact parameter and the influence of the relative translational energy on the reaction cross-section of the reaction have been discussed in detail. And the fact is found by the comparison that the reaction probability and the reaction cross-section of the reaction estimated in this work are more reasonable than the theoretical ones determined by Gilibert et al.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.6
• /
• pp.116-122
• /
• 2011

A realtime simulator using an explicit integration method is introduced to improve the solving performance for the dynamic analysis of a wheeled vehicle. Because a full vehicle system has many parts, the development of a numerical technique for multiple d.o.f. and ground contacts has been required to achieve a realtime dynamics analysis. This study proposes an efficient realtime solving technique that considers the wheeled vehicle dynamics behavior with full degrees of freedom and wheel contact with soft ground such as sand or undersea ground. A combat vehicle was developed to verify this method, and its dynamics results are compared with commercial programs using implicit integration methods. The combat vehicle consists of a chassis, double wishbone type front and rear suspension, and drive train. Some cases of vehicle dynamics analysis are carried out to verify the realtime ratio.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.5
• /
• pp.1551-1562
• /
• 1991

최근까지 발표된 유한 요소법을 이용한 압연 해석 관련 주요 논문들을 정리해 보면 다음과 같다. Li와 Kobayashil는 강소성 유한 요소법(rigidplastic finite element method)을 여러가지 마찰조건에 대하여 해석하였다. 이때 압연롤은 강체 (rigid body)로 시편은 가공경화(workhardening)를 동반한 강소성체로 모델링하였다. Hwang과 Kobayashi는 강소성 유한 요소법을 이용한 평면 변형 압연에서 재료 손실을 최소화하는 예비 성형체(preform)의 설계에 대한 연구를 수행하였다. 이 경우에도 역시 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체와 완전 소성체로 모델링 되었으나, 고착(sticking) 마찰 조건에 대해서만 해석을 수행하였다. Mori와 Osak- ada 그리고 Oda는 약간 압축성이 있는 재료의 평면 변형 압연에 대하여 연구하였다. 이때 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체로 모델링 되었으며 경계 면에서는 Coulomb 마찰을 고려하였다. 이밖에도 오일러(Eulerian) 수식화를 이용한 Dawson과 Thompson, Berman의 해석 결과가 있으며, 또 폭 방향의 변형까지를 고려한 Li와 Kobayashi, Mori와 Osakada의 3차원 해석 결과가 있다.