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

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.407-412
• /
• 1998

In this work, a new time integration method is proposed using the generalized derivative concept to simulate the dynamic phenomena having sudden constraint occurring in dynamic contact/impact problems. By the adoption of the generalized derivative concept and jump assumption, discontinuity can be incorporated in time integration and as a result, the algorithm does not need any other special consideration of jumps in dynamic field variables due to sudden constraint like dynamic contact-release conditions. To observe the characteristics of the proposed time integration method, the stability and convergence analyses are carried out. In numerical tests, several dynamic contact/impact problems are analyzed by straightforward application of the proposed time integration method with the exterior penalty method.

• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.145-162
• /
• 1995

A damped trapezoidal rule method for numerical time-integration is presented, and its application in analyses of dynamic response of damped structures is discussed. It is shown that the damped trapezoidal rule method has features that make it an attractive approach for applications in dynamic analyses of structures. Accuracy and stability analyses are developed for the damped single-degree-of-freedom systems. Error analyses are also performed for the Newmark beta method and compared with the damped trapezoidal rule method as a basis for discussion of the relative merits of the proposed method. The procedure is fully explicit and easy to implement. However, since the method is an explicit method, it is conditionally stable. The methodology is applied to several example problems to illustrate its strengths, limitations and inherent simplicity.

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

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.756-761
• /
• 2004

The effects of two important numerical procedures on the high precision structural analysis are investigated in this study. The two numerical procedures include continuous variable approximation and time integration. For the continuous variable approximation, polynomial mode functions generated by the Gram-Schmidt process are introduced and the numerical results obtained by employing the polynomial mode functions are compared to those obtained by classical beam mode functions. The effect of the time integration procedure on the analysis precision is also investigated. It is found that the two procedures affect the precision of structural analysis significantly.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.28 no.6
• /
• pp.573-578
• /
• 2010

Precise determination of satellite positions is necessary to improve positioning accuracy in GNSS. In this study, GLONASS orbits were predicted from broadcast ephemeris using the 4th-order Runge-Kutta numerical integration method and their accuracy dependence on the integration step and the integration time was analyzed. The 3D RMS (Root Mean Square) differences between the results from I-second integration step and 300-second integration step was about 3 cm, but the processing time was one hundred times less for the I-second integration time case. For trials of different integration times, the 3D RMS errors were 8.3 m, 187.3 m, and 661.5 m for 30-, 150-, and 300-minutes of integration time, respectively. Though this integration-time analysis, we concluded that the accuracy gets higher with a shorter integration time. Thus we suggest forward and backward integration methods to improve GLONASS positioning accuracy, and with this method we can achieve a 5-meter level of 3-D orbit accuracy.

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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.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.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.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.