• Journal of the Korean Society for Railway
• /
• v.15 no.1
• /
• pp.1-8
• /
• 2012

Wheelset dynamic analysis is a key element to determine the degree of accuracy of railway vehicle dynamics. In this study, a three-dimensional wheelset dynamic analysis is presented in such a way that the precise wheel-rail contact analysis in three-dimension is implemented into the dynamic equations of a wheelset. A numerical procedure that can be used for the analysis of a wheelset dynamics when the wheel-rail two point contact occurs in a cornering maneuver is developed. Numerical solutions of the constraint equations and the dynamics equations of a wheelset are achieved by using Runge-Kutta method. The proposed wheelset dynamic analysis is validated by comparison against results obtained from VI-RAIL analysis.

• Journal of Astronomy and Space Sciences
• /
• v.17 no.1
• /
• pp.107-116
• /
• 2000

This paper presents the preliminary survey and simulation results of the prediction of Leonid stream's orbital motion. Based on the model survey on eject velocity and perturbation of meteoroid particles, a simulation program was developed and applied to orbital motion of Leonid stream. The Jones ejection distribution model was used to describe the particle's eject velocity and the orbital dynamic model includes perturbations of major planet's gravity. DE405 ephemeris file generated by Solar System Dynamics Group at Jet Propulsion Laboratory in NASA was used for the planet's ephemeris calculations. Solar radiation pressure were also considered in the simulation and 8th order Runge-Kutta algorithm was used a numerical integration method.

• Journal of Astronomy and Space Sciences
• /
• v.26 no.2
• /
• pp.199-210
• /
• 2009

Even though there are several Global Navigation Satellite Systems under development, only GPS and GLONASS are currently available for satellite positioning. In this study, GLONASS orbits were predicted from broadcast ephemeris using the 4th-order Runge-Kutta numerical integration. For accuracy validation, predicted orbits were compared with precise ephemeris. The RMS(Root Mean Square) and maximum 3-D errors were 14.3 km and 17.4 km for one-day predictions. In case of 7-day predictions, the RMS and maximum 3-D errors were 15.7 and 40.1 km, respectively. Also, the GLONASS satellite visibility predictions were compared with real observations, and they agree perfectly except for several epochs when the satellite signal was blocked by nearby buildings.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.12
• /
• pp.5689-5695
• /
• 2012

In this study, vibration properties of seat and human body are analyzed through test and numerical analysis methods by taking into account the viscoelastic characteristics of polyurethane foam as seat material which is applied for vehicle. These viscoelastic characteristics which show nonlinear and quasi-static behavior are obtained by compression test. In addition, the viscous elastic property of polyurethane foam is modelled mathematically by using convolution integral and nonlinear stiffness model. In order to analyze the performance on ride comfort of seat, vertical vibration model is established by using dynamic model of seat and vertical vibration model of human body at ISO5982, and so the related motion equations are derived. A numerical analysis simulation is applied by using the nonlinear motion equation with Runge-Kutta integral method. The dynamic responses of seat and human body on the input of vibration acceleration measured at the floor of the railway vehicle are examined. The variation of the index value at ride comfort on seat design parameters is analyzed and the methodology on seat design is suggested.

• Journal of the Institute of Electronics Engineers of Korea SC
• /
• v.46 no.6
• /
• pp.73-79
• /
• 2009

For the shake of control for movement object, control theory like neural network, nonlinear model predictive control(NMPC) is realized on digital high speed computer. Predictor of flight control system(FCS) based nonlinear model predictive control has to be satisfied with response for hard real-time to perform applications on each module in the FCS. Simultaneously, It gives a serious consideration accuracy to give full play to FCS's performance. Error of mathematical aspect affects realization of whole algorithm. But factors of bring mathematical error is not considered to calculate final accuracy on parameter of predictor. In this paper, Predictor was made using load control model on the digital computer for design FCS at hard real-time and is shown response time on realization algorithm. And is shown realization algorithm of high effective predictor over the accuracy. The predictor was realized on the load control model using Euler method, Heun method, Runge-Kutta and Taylor method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.608-614
• /
• 1997

2차원 압축성 오일러 방정식을 지배 방정식으로 하여 음향학적 문제를 직접 수치 계산하였다. Roe의 개략적인 리만해(Approximate Riemann Solver)에 의거하여 대류항을 구성하였고 높은 정확도를 유지하기 위해 Vanleer의 MUSCL 방법을 사용하여 공간을 이산화하였다. 시간 적분으로는 2차의 Runge-Kutta 방법을 사용하였다. Resonance tube 내의 음향학적 문제 해석을 위해 최근까지 제시된 음향 경계 조건들을 비교하였고 각 경계 조건들이 갖는 물리적 의미를 고찰하였다. 이들을 통하여 압력파의 물리적인 행태를 모사하였고 실험과 비교하였다. 계산결과들을 통하여 사용된 방법은 전형적인 압축성 유동 해석 문제 뿐만 아니라 선형, 비선형의 음향학적 문제들을 적절한 경계 조건과 결합하여 해석할 수 있음을 알 수 있었고 resonance tube와 같은 공학적으로 실제적인 문제들에 대한 응용이 기대된다.

• Journal of Astronomy and Space Sciences
• /
• v.2 no.1
• /
• pp.19-33
• /
• 1985

The singularities of differential Newtonian equation of motion in three body problem cause the loss of accuracy and the considerable increase of the computer time. These singularities could be eliminated during the process of regularization to transform the independent variables and the coordinate of Newtonian equations of motion. In this study, we calculated the positions and velocities of three body along the time scale to find out the unique solution of regularized Newtonian equations of motion with the $5^{th}$ order Runge-Kutta method by assuming the suitable initial velocities and positions. As the results of these calculations it is shown that the tripe stellar system eventually distintegrated, two of them formed a binary, and the last one escaped from this system with a hyperbolic orbit. This may suggest one possible explanation for the binary formation.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.38-52
• /
• 1991

In this paper, a semi-Lagrangian method is used to solve the nonlinear hydrodynamics of a three-dimensional body beneath the free surface in the time domain. The boundary value problem is solved by using the boundary integral method. The geometries of the body and the free surface are represented by the curved panels. The surfaces are discretized into the small surface elements using a bi-cubic B-spline algorithm. The boundary values of $\phi$ and $\frac{\partial{\phi}}{\partial{n}}$ are assumed to be bilinear on the subdivided surface. The singular part proportional to $\frac{1}{R}$ are subtracted off and are integrated analytically in the calculation of the induced potential by singularities. The far field flow away from the body is represented by a dipole at the origin of the coordinate system. The Runge-Kutta 4-th order algorithm is employed in the time stepping scheme. The three-dimensional form of the integral equation and the boundary conditions for the time derivative of the potential Is derived. By using these formulas, the free surface shape and the equations of motion are calculated simultaneously. The free surface shape and fille forces acting on a body oscillating sinusoidally with large amplitude are calculated and compared with published results. Nonlinear effects on a body near the free surface are investigated.

• Computational Structural Engineering
• /
• v.7 no.3
• /
• pp.95-103
• /
• 1994

The safety related components in the nuclear power plant should be designed to withstand the seismic load. Among these components the integrity of reactor internals under earthquake load is important in stand points of safety and economics, because these are classified to Seismic Class I components. So far the modelling methods of reactor internals have been investigated by many authors. In this paper, the dynamic behaviour of reactor internals of Yong Gwang 1&2 nuclear power plants under SSE(Safe Shutdown Earthquake) load is analyzed by using of the simpled Global Beam Model. For this, as a first step, the characteristic analysis of reactor internal components are performed by using of the finite element code ANSYS. And the Global Beam Model for reactor internals which includes beam elements, nonlinear impact springs which have gaps in upper and lower positions, and hydrodynamical couplings which simulate the fluid-filled cylinders of reactor vessel and core barrel structures is established. And for the exciting external force the response spectrum which is applied to reactor support is converted to the time history input. With this excitation and the model the dynamic behaviour of reactor internals is obtained. As the results, the structural integrity of reactor internal components under seismic excitation is verified and the input for the detailed duel assembly series model could be obtained. And the simplicity and effectiveness of Global Beam Model and the economics of the explicit Runge-Kutta-Gills algorithm in impact problem of high frequency interface components are confirmed.