• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11b
• /
• pp.1249-1253
• /
• 2001

Numerical analysis is sometimes used to solve the problems in the engineering and natural science fields. On this reason, the faster, more practical system in computing the numerical solution is required. This paper deals with the numerical evaluation of various numerical integration methods which is frequently used in the engineering fields. This paper choices four integration methods such as Euler method, Heun's method, Runge-Kutta method and Gill's method for evaluating the each integration method. In numerical examples, the free vibration problem on an elastic foundation is chosen. As the numerical results, the natural frequencies and the running time are obtained, and these results are compared to examine the practicality of integration methods.

• The Mathematical Education
• /
• v.37 no.1
• /
• pp.101-114
• /
• 1998

In this thesis, We tried to introduce definite integration to the curriculum of high school mathematics with numerical integration, which had been introduced with quadrature method. For this purpose, We used new experimental mathematics approaches, so-called investigation and examination. In chapter II, We examined how much computers had been used in teaching mathematics. In chapter III, We presented the theoretical background of approximation integration within numerical integration. In chapter IV, We studied and compared various methods of numerical integration, and examined the relation between curvature of a curved line and numerical integration. In order to study more easily, We used some of computer programs. We hope that this thesis will be a turning point in developing new teaching methods and improving curriculum of mathematics in high school.

• Bulletin of the Korean Space Science Society
• /
• 2008.10a
• /
• pp.26.4-27
• /
• 2008

The completion ('initiation' de facto) of the KASI Orbit Propagator and Estimator (KASIOPEA) has been delayed for several reasons unfortunately. Due to the lack of working staffs and the Division priority rearrangement, the initial plan was dismantled and ignored for many years. However, fundamental researches regarding the essential parts of KASIOPEA has been done by author. The numerical integration module of the KASIOPEA is the most sensitive part in the precision of the final output in general. There is no silver bullet in the numerical integration in an orbit propagation as a non-stiff ODE case. Many numerical integration method like single-step methods, multi-step method, and extrapolation methods have been used in overly populated orbit propagator or estimator. In this study, several popular methods from single-step, multi-step, and extrapolation methods have been tested in numerical accuracy and stability.

• Structural Engineering and Mechanics
• /
• v.75 no.5
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.

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

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• 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.

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

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.