• Journal of KSNVE
• /
• v.7 no.5
• /
• pp.773-780
• /
• 1997

Using Generalized Inverse Method presented by Udwadia and Kalaba in 1992, we can obtain equations to exactly describe the motion of constrained systems. When the differential equations are numerically integrated by any numerical integration scheme, the numerical results are generally found to veer away from satisfying constraint equations. Thus, this paper deals with the numerical integration of the differential equations describing constrained systems. Based on Baumgarte method, we propose numerical methods for reducing the errors in the satisfaction of the constraints.

• Journal of Positioning, Navigation, and Timing
• /
• v.8 no.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.

• Geosystem Engineering
• /
• v.21 no.6
• /
• pp.309-317
• /
• 2018

We employ a three-dimensional indirect boundary element method (BEM) to simulate temperature change around an underground liquefied natural gas storage cavern. The indirect BEM (IBEM) uses fictitious heat source strength on boundary elements as basic variables which are solved from equations of boundary conditions and then used to compute the temperature change at other points in the considered problem domain. The IBEM requires evaluation of singular integration for temperature change due to heat conduction from a constant heat source on a planar (triangular) region. The singularity can be eliminated by a semi-analytical integration scheme. However, it is found that the semi-analytical integration scheme yields sharp temperature gradient for points close to vertices of triangle. This affects the accuracy of heat flux, if they are evaluated by finite difference method at these points. This difficulty can be overcome by a combination of using a direct numerical integration for these points and the semi-analytical scheme for other points distance away from the vertices. The IBEM and the hybrid integration scheme have been verified with an analytic solution and then used to the application of the underground storage.

• Atmosphere
• /
• v.33 no.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%.

• Smart Structures and Systems
• /
• v.14 no.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.

• 한국기계연구소 소보
• /
• s.18
• /
• pp.55-66
• /
• 1988

Many direct integration methods are used for numerical analyses of dynamic motion. In these methods, the governing equations of a dynamic system are integrated successively using a step-by-step numerical integration procedure. Time derivatives in the equations are generally approximated using difference formulas involving one or more increments of the time. Time increment has closely relationship with the accuracy of the motion analysis. In this paper, a 4th order Houbolt direct integration method is derived. For a spring-mass system, the motion of the system are analyzed from the 3rd order Houbolt and the 4th order Houbolt approaches respectively. Finally the paper proposes the optimal time-increment based on the accuracy of numerical analyses.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1111-1117
• /
• 2008

During the dynamic analysis of a system, the Coulomb friction law is emploved to calculate the friction force. Since the static friction coefficient is only employed during the zero relative velocity, it is impractical to employ the coefficient during the dynamic analysis. To calculate the static friction force, therefore, some friction models have been developed. In this study, the integration stability and the accuracy of the models are investigated with some numerical examples. The effect of time step size during the numerical integration is also investigated. The numerical study shows that the friction model employed for most commercial codes is not as good as the one proposed in this study.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1719-1730
• /
• 1993

This paper presents the theoretical natural frequencies of out-of-plane deformable ring based on the variables such as out-of-plane deflection, torsional rotation and shear rotation. Based on the same variables, a finite element eigen analysis is carried out by using the $C^0$-continuous, isoparametric element which has three nodes per element and three degrees-of-freedom at each node. Numerical experiments are peformed to find the integration scheme which produces accurate natural frequencies, natural modes and correct rigid body motion. The uniformly reduced integration and the selective reduced integration give more accurate numerical frequencies than the uniformly full integration, but the uniformly reduced integration produces incorrect rigid body motion while selective reduced integration does correct one. Therefore, the ring element based on the three variables which employes selective reduced integration is recommended to avoid spurious modes, to alleviate the error due to shear locking and to produce correct rigid body motion, simultaneously.

• Earthquakes and Structures
• /
• v.13 no.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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.10
• /
• pp.1605-1612
• /
• 2001

Least-squares meshfree method is presented. Conventional meshfree methods based on the Galerkin formulation suffer from inaccurate numerical integration. Least-squares formulation exhibits rather different integration-related characteristics. It is demonstrated through numerical examples that least-squares formulation is much more robust to integration errors than the Galerkin's. Therefore efficient meshfree methods can be devised by combining very simple integration algorithms and least-squares formulation.