• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.453-458
• /
• 2004

This study proposes a new two-level condensation scheme for the construction of a reduced system. In the first step, the candidate area is selected for the construction of the reduced system by energy estimation in element-level. In the second step, primary degrees of freedom are selected by sequential elimination from the candidate degrees of freedom linked to the selected elements. Numerical examples demonstrate that the proposed method saves the computational cost effectively and provides a reduced system which predicts the eigenvalues accurately. Moreover, the well-constructed reduced system can present the reliable behavior of the structure under arbitrary dynamic loads comparing to that of global system. Time integration in a reduced system can save the computing time remarkably. Through a few numerical examples, the efficiency and reliability of the proposed scheme are verified.

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

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

• Earthquakes and Structures
• /
• v.11 no.2
• /
• pp.297-313
• /
• 2016

A virtual parameter is introduced into the formulation of the previously published integration method to improve its stability properties. It seems that the numerical properties of this integration method are almost unaffected by this parameter except for the stability property. As a result, it can have second order accuracy, explicit formulation and controllable numerical dissipation in addition to the enhanced stability property. In fact, it can have unconditional stability for the system with the instantaneous degree of nonlinearity less than or equal to the specified value of the virtual parameter for the modes of interest for each time step.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.589-600
• /
• 1996

An explicit 4th order time integration scheme for solving the convection-diffusion equation is discussed in this paper. A system of ordinary differential equations are derived first by discretizing the spatial derivatives of the relevant PDE using the finite difference method. The integration of the ODEs is then carried out using a 4th order scheme and a self-adaptive technique based on the spatial grid spacing. For a non-uniform spatial grid, different time step sizes are used for the integration of the ODEs defined at different spatial points, which improves the computational efficiency significantly. A numerical example is also discussed in the paper to demonstrate the implementation and effectiveness of the method.

• Journal of Mechanical Science and Technology
• /
• v.20 no.3
• /
• pp.376-381
• /
• 2006

In this paper, the quaternion, the dual Euler, and the direction cosine methods are numerically compared using a non-aerodynamic 6 degree-of-freedom rigid model at all-attitude angles of an aircraft. The dual Euler method turns out to be superior to the others in the applications because it shows better numerical accuracy, stability, and robustness in integration step sizes. The dual Euler method is affordably less efficient than the quaternion method in terms of computational cost. Numerical accuracy and stability, which allow larger integration step sizes, are more critical in modern real-time applications than computational efficiency because of today's increased computational power. If the quaternion method is required because of constraints in computation time, then a suppression mechanism should be provided for algebraic constraint errors which will eventually add computational burden.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.19-26
• /
• 2004

This study constructs the reduced system by two-level condensation scheme. This scheme consists of two steps. First step selects the candidate area for the primary degrees of freedom by energy estimation in element level. In the second step, the primary degrees of freedom are selected by the sequential elimination scheme. The efficiency and reliability of this scheme is shown through the prediction of eigenvalues of a few numerical examples. Time integration in the reduced system can save the computing time effectively. The well-constructed reduced system can present the accurate behavior of the structure under arbitrary dynamic loads so much as the global system. Through the numerical example, the efficiency and reliability of the proposed scheme will be demonstrated.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

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

• Structural Engineering and Mechanics
• /
• v.41 no.5
• /
• pp.675-689
• /
• 2012

Most connections of steel structures exhibit flexible behaviour under cyclic loading. The flexible connections can be assumed as nonlinear rotational springs attached to the ends of each beam. The nonlinear behaviour of the connections can be considered by suitable moment-rotation relationship. Time-history analysis by direct integration method can be used as a powerful technique to determine the nonlinear dynamic response of the structure. In conventional numerical integration, the response is evaluated for a series of short time increments. The limitations on the size of time intervals can be removed by using Chen and Robinson improved time history analysis method, in which the integrated displacements are used as the new variables in integrated equation of motion. The proposed method permits longer time intervals and reduces the computational works. In this paper the nonlinearity behaviour of the structure is summarized on the connections, and the step by step improved time-history analysis is used to calculate the dynamic response of the structure. Several numerical calculations which indicate the applicability and advantages of the proposed methodology are presented. These calculations illustrate the importance of the effect of the nonlinear behaviour of the flexible connections in the calculation of the dynamic response of steel frames.