• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.389-394
• /
• 2006

Fully flexible cell preserves Hamiltonian in structure, so the symplectic time integrator is applied to the equations of motion. Primarily, generalized leapfrog time integration (GLF) is applicable, but the equations of motion by GLF have some of implicit formulas. The implicit formulas give rise to a complicate calculation for coding and need an iteration process. In this paper, the time integration formulas are obtained for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term, so the simple and completely explicit recursion formula was obtained. The explicit formulas are easy to implementation for coding and may be reduced the integration time because they are not need iteration process. We are going to compare the resulting splitting time integration with the implicit generalized leapfrog time integration.

• 대한기계학회논문집A
• /
• 제31권8호
• /
• pp.826-832
• /
• 2007

Fully flexible cell preserves Hamiltonian in structure so that the symplectic time integrator is applicable to the equations of motion. In the direct formulation of fully flexible cell N-Sigma-T ensemble, a generalized leapfrog time integration (GLF) is applicable for fully flexible cell simulation, but the equations of motion by GLF has structure of implicit algorithm. In this paper, the time integration formula is derived for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term. Thus the simple and completely explicit recursion formula was obtained. We compare the performance and the result of present splitting time integration with those of the implicit generalized leapfrog time integration.

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

• 대한기계학회논문집A
• /
• 제27권11호
• /
• pp.1907-1916
• /
• 2003

During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

• Structural Engineering and Mechanics
• /
• 제28권5호
• /
• pp.549-570
• /
• 2008

Numerical integration is an efficient approach for nonlinear dynamic analysis. In this paper, general category of the implicit integration errors will be discussed. In order to decrease the errors, Dynamic Relaxation method with modified time step (MFT) will be used. This procedure leads to an alternative algorithm which is very general and can be utilized with any implicit integration scheme. For numerical verification of the proposed technique, some single and multi degrees of freedom nonlinear dynamic systems will be analyzed. Moreover, results are compared with both exact and other available solutions. Suitable accuracy, high efficiency, simplicity, vector operations and automatic procedures are the main merits of the new algorithm in solving nonlinear dynamic problems.

• 한국지진공학회논문집
• /
• 제10권5호
• /
• pp.99-106
• /
• 2006

이 논문에서 소개되는 실시간 하이브리드 실험 시스템은 유사동적실험법을 근거로 하고 있으며, 실제 실험과수학적 모델을 이용한 수치해석을 결합한 실험법이다. 이 시스템은 종래의 유사동적해석법에 비해 지진하중을 받는 구조물의 실시간 반응에 근접하는 현저히 높은 하중재하 속도를 고려할 수 있도록 설계되었다. 또한 다자유도 구조물에 대해 안정적인 해석환경을 제공하기 위해 이 시스템은 implicit 시간적분법을 이용하여 수치해석을 수행한다. 이 논문은 연구를 통해 개발된 시스템의 전반적인 개요와 구성요소 그리고 이 시스템에서 사용하는 수치 해석법의 성능을 평가하기 위해 수행된 수치해석을 소개한다. 연구 결과 개발된 시스템에 적용된 수치해석법은 성능이 매우 우수하다는 것이 증명되었다.

• Structural Engineering and Mechanics
• /
• 제48권5호
• /
• pp.711-736
• /
• 2013

In this paper new implicit time integration called N-IHOA is presented for dynamic analysis of high damping systems. Here, current displacement and velocity are assumed to be functions of the velocities and accelerations of several previous time steps, respectively. This definition causes that only one set of weighted factors is calculated from the Taylor series expansion which leads to a simple approach and reduce the computational efforts. Moreover a comprehensive study on stability of the proposed method i.e., N-IHOA compared with IHOA integration which is performed based on amplification matrices proves the ability of the N-IHOA in high damping vibrations such as control systems. Also, wide range of numerical examples which contains single/multi degrees of freedom, damped/un-damped, free/forced vibrations from finite element/finite difference demonstrate that the accuracy and efficiency of the proposed time integration is more than the common approaches such as the IHOA, the Wilson-${\theta}$ and the Newmark-${\beta}$.

• Wind and Structures
• /
• 제34권1호
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회A
• /
• pp.512-517
• /
• 2007

In the previous development of the recursive thermostat chained fully flexible cell molecular dynamics simulation, implicit time integration method such as generalized leapfrog integration is used. The implicit algorithm is very much complicated and not easy to show time reversibility because it is solved by the nonlinear iterative procedure. Thus we develop simple, explicit symplectic time integration formula for the recursive thermostat chained fully flexible unit cell simulation. Uniaxial tension test is performed to verify the present explicit algorithm. We check that the present simulation satisfies the ergodic hypothesis for various values of fictitious mass and coefficient of multiple thermostat system. The proposed method should be helpful to predict mechanical and thermal behavior of nano-scale structure.

• 한국항공우주학회지
• /
• 제31권7호
• /
• pp.1-8
• /
• 2003

본 논문은 극초음속 유동과 공력가열 해석에서 나타나는 공기 반응의 수치해석적 특징을 다루고 공기반응을 효율적으로 해석하는 부분 내재적 적분법을 공기반응에 적용하였다. 안정적 계산을 위해 화학반응 자코비안이 필수적임을 밝혔으며 자코비안의 양의 실수 특성치로 인한 수치기법의 경직성은 일반적인 연소반응에 비하여 미약하였다. 공기반응에서 부분 내재적 적분법은 화학종 순서의 종속성이 없었으며 완전 내재적 적분법과 동일한 수렴율과 계산 결과를 보였다. 극초음속 유동해석에 부분 내재적 적분법을 적용하면 전체 연산 시간이 감소되었다.