• Earthquakes and Structures
• /
• 제13권3호
• /
• pp.279-288
• /
• 2017

Having the ability to keep on yielding stable solutions in problems involving high potential of instability, composite time integration methods have become very popular among scientists. These methods try to split a time step into multiple sub-steps so that each sub-step can be solved using different time integration methods with different behaviors. This paper proposes a new composite time integration in which a time step is divided into two sub-steps; the first sub-step is solved using the well-known Newmark method and the second sub-step is solved using Simpson's Rule of integration. An unconditional stability region is determined for the constant parameters to be chosen from. Also accuracy analysis is perform on the proposed method and proved that minor period elongation as well as a reasonable amount of numerical dissipation is produced in the responses obtained by the proposed method. Finally, in order to provide a practical assessment of the method, several benchmark problems are solved using the proposed method.

• 대기
• /
• 제33권4호
• /
• pp.331-341
• /
• 2023

A numerical forecasting models usually predict future states by performing time integration considering fixed static time-steps. A time-step that is too long can cause model instability and failure of forecast simulation, and a time-step that is too short can cause unnecessary time integration calculations. Thus, in numerical models, the time-step size can be determined by the CFL (Courant-Friedrichs-Lewy)-condition, and this condition acts as a necessary condition for finding a numerical solution. A static time-step is defined as using the same fixed time-step for time integration. On the other hand, applying a different time-step for each integration while guaranteeing the stability of the solution in time advancement is called an adaptive time-step. The adaptive time-step algorithm is a method of presenting the maximum usable time-step suitable for each integration based on the CFL-condition for the adaptive time-step. In this paper, the adaptive time-step algorithm is applied for the Korean Integrated Model (KIM) to determine suitable parameters used for the adaptive time-step algorithm through the monthly verifications of 10-day simulations (during January and July 2017) at about 12 km resolution. By comparing the numerical results obtained by applying the 25 second static time-step to KIM in Supercomputer 5 (Nurion), it shows similar results in terms of forecast quality, presents the maximum available time-step for each integration, and improves the calculation efficiency by reducing the number of total time integrations by 19%.

• 대한기계학회논문집A
• /
• 제21권10호
• /
• pp.1712-1717
• /
• 1997

Stability and accuracy for the trapezoidal rule of the Newmark time integration method are analyzed when variable time step sizes are adopted. A new analytic approach to stability and accuracy analysis is also proposed for time integration methods with variable time step sizes. The trapezoidal rule with variable time step sizes has the "actual" unconditional stability which is the same as that of the method with constant time step sizes. However, the method with variable time step sizes is first-order accurate while the method with constant time step sizes is second-order accurate. accurate.

• Earthquakes and Structures
• /
• 제8권6호
• /
• pp.1325-1348
• /
• 2015

Real-time hybrid testing (RTHT) involves virtual splitting of the structure into two parts: physical substructure that contains the key region of interest which is tested in a laboratory and numerical substructure that contains the remaining part of the structure in the form of a numerical model. This paper numerically assesses four step-by-step integration methods (Central difference method (CDM), Operator splitting method (OSM), Rosenbrock based method (RBM) and CR-integration method (CR)) which are widely used in RTHT. The methods have been assessed in terms of stability and accuracy for various realistic damping ratios of the physical substructure. The stability is assessed in terms of the spectral radii of the amplification matrix while the accuracy in terms of numerical damping and period distortion. In order to evaluate the performance of the methods, five carefully chosen examples have been studied - undamped SDOF, damped SDOF, instantaneous softening, instantaneous hardening and hysteretic system. The performance of the methods is measured in terms of a non-dimensional error index for displacement and velocity. Based on the error indices, it is observed that OSM and RBM are robust and performs fairly well in all the cases. CDM performed well for undamped SDOF system. CR method can be used for the system showing softening behaviour. The error indices indicate that accuracy of OSM is more than other method in case of hysteretic system. The accuracy of the results obtained through time integration methods for different damping ratios of the physical substructure is addressed in the present study. In the presence of a number of integration methods, it is preferable to have criteria for the selection of the time integration scheme. As such criteria are not available presently, this paper attempts to fill this gap by numerically assessing the four commonly used step-by-step methods.

• 대한기계학회논문집A
• /
• 제24권11호
• /
• pp.2804-2811
• /
• 2000

A Nordsick form opf the multirate integration scheme has been proposed for flexible multibody dynamic systems. It is assumed that vibrational modal coordinates in the equations of motion are treated as fast variables, whereas the relative joint coordinates are treated as slow variables. In the multirate integration, the fast variables are integrated with small step-size, and the slow variables are integrated with larger step-size. The proposed multirate integration method is based on the Adams-Bashforth-Moulton predictor-corrector method and implemented in the Nordsieck vector form. The Nordsieck form of multrate integration method provides effective step-size control and at the same time, inherits the efficiency from the Adams integration method. Simulations of a flexible gun and turret system of the military tank have been carried out to show the effectiveness and efficiency of the proposed method.

• Smart Structures and Systems
• /
• 제14권6호
• /
• pp.1269-1289
• /
• 2014

With the sub-stepping technique, the numerical analysis in real-time dynamic hybrid testing is split into the response analysis and signal generation tasks. Two target computers that operate in real-time may be assigned to implement these two tasks, respectively, for fully extending the simulation scale of the numerical substructure. In this case, the integration time-step of solving the dynamic response of the numerical substructure can be dozens of times bigger than the sampling time-step of the controller. The time delay between the real and desired feedback forces becomes more striking, which challenges the well-developed delay compensation methods in real-time dynamic hybrid testing. This paper focuses on displacement prediction and force correction for delay compensation in the real-time dynamic hybrid testing with a large integration time-step. A new displacement prediction scheme is proposed based on recently-developed explicit integration algorithms and compared with several commonly-used prediction procedures. The evaluation of its prediction accuracy is carried out theoretically, numerically and experimentally. Results indicate that the accuracy and effectiveness of the proposed prediction method are of significance.

• Journal of Positioning, Navigation, and Timing
• /
• 제8권4호
• /
• pp.201-208
• /
• 2019

Numerical integration is necessary for satellite orbit determination and its prediction. The numerical integration algorithm can be divided into single-step and multi-step method. There are lots of single-step and multi-step methods. However, the Runge-Kutta method in single-step and the Adams method in multi-step are generally used in global navigation satellite system (GNSS) satellite orbit. In this study, 4th and 8th order Runge-Kutta methods and various order of Adams-Bashforth-Moulton methods were used for GLObal NAvigation Satellite System (GLONASS) orbit integration using its broadcast ephemeris and these methods were compared with international GNSS service (IGS) final products for 7days. As a result, the RMSE of Runge-Kutta methods were 3.13m and 4th and 8th order Runge-Kutta results were very close and also 3rd to 9th order Adams-Bashforth-Moulton results. About result of computation time, this study showed that 4th order Runge-Kutta was the fastest. However, in case of 8th order Runge-Kutta, it was faster than 14th order Adams-Bashforth-Moulton but slower than 13th order Adams-Bashforth-Moulton in this study.

• 한국자동차공학회논문집
• /
• 제4권4호
• /
• pp.190-200
• /
• 1996

Step-by-step time integration methods are widely used for solving structural dynamics problem. One difficult yet critical choice an analyst must make is to decide an appropriate time step size. The choice of time step size has a significant effect on solution accuracy and computational expense. The objective of this research is to derive error estimate for newly developed time integration method and develop automatic time step size control algorithm for structural dynamics. A formula for computing error tolerance is derived based on desired period resolution. An automatic time step size control strategy is proposed based on a normalized local error estimate for the generalized-α method. Numerical examples demonstrate the developed strategy satisfies general design criteria for time step size control algorithm for dynamic problem.

• 한국측량학회지
• /
• 제28권6호
• /
• pp.573-578
• /
• 2010

위성항법시스템에서 정확한 위성궤도결정 기술은 측위 정확도 향상의 필수적인 조건이다. 이 연구에서는 GLONASS의 방송궤도력과 4차 Runge-Kutta 수치적분법을 이용하여 위성좌표를 결정하였으며, 적분간격과 적분시간에 따른 위성궤도의 정확도를 비교하였다. 적분간격에 따른 위성궤도 정확도분석결과, 적분간격이 l초일 때와 300초일 때의 3차원 RMS 오자의 차이가 3cm에 불과한 반면 처리시간은 100배 이상 향상되었다. 적분시간에 따른 위성좌표의 3차원 RMS 오차는 적분시간이 30분, 150분, 300분일 때 각각 8.3m, 187.3m, 661.5m로 나타났으며, 이를 통해 적분시간을 짧게 할수록 정확도가 향상되는 것을 확인하였다. 따라서 이 연구에서는 GLONASS 측위를 위한 위성좌표 결정의 정확도 향상을 위해 적분시간을 최소화할 수 있는 Forward와 Backward 적분을 적용하는 방안을 제안하였으며, 이와 같은 방법을 사용할 경우 5m이하의 위성좌표 산출 정확도를 확보할 수 있다.

• Structural Engineering and Mechanics
• /
• 제13권5호
• /
• pp.569-589
• /
• 2002

In performing the dynamic analysis, the step size used in a step-by-step integration method might be much smaller than that required by the accuracy consideration in order to capture the rapid chances of dynamic loading or to eliminate the linearization errors. It was first found by Chen and Robinson that these difficulties might be overcome by integrating the equations of motion with respect to time once. A further study of this technique is conducted herein. This include the theoretical evaluation and comparison of the capability to capture the rapid changes of dynamic loading if using the constant average acceleration method and its integral form and the exploration of the superiority of the time integration to reduce the linearization error. In addition, its advantage in the solution of the impact problems or the wave propagation problems is also numerically demonstrated. It seems that this time integration technique can be applicable to all the currently available direct integration methods.