• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.5 no.4
• /
• pp.243-248
• /
• 1993

In this research, unsteady groundwater flow in unconfined and homogeneous three layer aquifers is studied theoretically and experimentally. Numerical solutions are obtained by Runge Kutta and Runge Kutta Gill method after transforming the governing nonlinear partial differential equations to nonlinear ordinary differential equations. Experimental apparatus includes a test section filled with fine, medium and coarse sands. Experimental results are compared with the numerical solutions and both experimental and numerical results correspond well with each other. This numerical approach may be also applied to the cases which have more aquifers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.2
• /
• pp.101-107
• /
• 2012

In this study, Constraint Force Equation Methodology (CFEM) is used to develop a multi-body dynamic analysis program with constraints. Seven constraint models are implemented to analyze constraint motions of multiple bodies. The augmented equations with the constraints are solved with the 4th order Runge-Kutta method for higher degree of accuracy. The analysis code is verified by comparing the analysis results of the motion of bodies with various constraints to published results.

• Steel and Composite Structures
• /
• v.17 no.5
• /
• pp.721-734
• /
• 2014

In the current study, we consider a new class of analytical periodic solution for free nonlinear vibration of mechanical systems. Hamiltonian approach is applied to analyze nonlinear problems which occur in dynamics. The proposed method doesn't have the limitations of the classical methods and leads us to a high accurate solution by only one iteration. Two well known examples are studied to show the convenience and effectiveness of this approach. Runge-Kutta's algorithm is also applied and the results of it are compared with the Hamiltonian approach. High accuracy of the proposed approach reveals that the Hamiltonian approach can be very useful for other nonlinear practical problems in engineering.

• Applied Chemistry for Engineering
• /
• v.9 no.7
• /
• pp.1079-1084
• /
• 1998

Reaction kinetics between 4,4'-dihexyl methane diisocyanate($H_{12}MDI$) and n-hexanol in toluene with dibutyltin dilaurate(DBTDL) as catalyst was studied by experimental measurements and mathematical modeling. Experiments were carried out at various temperatures, catalyst concentrations and [NCO]/[OH] ratios, and the reaction kinetics were described by two second-order reactions, the one between NCO and OH leading to urethane and the other between urethane and NCO leading to allophanate. The rate constants were estimated by the Runge-Kutta 4th-order method. Experiments and mathematical simulations showed a good agreement for various experimental conditions. The [allophanate]/[urethane] ratios at 90% conversion of initial NCO were estimated to be over 20% for most conditions employed in the present study, indicating that allophanate formation might significantly affect the properties of urethane polymers.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.2
• /
• pp.95-101
• /
• 2017

In this study, we derived theoretical equations for the zigzag movement of a leading vehicle, which is the most frequent behavior in train accidents, by using a simplified spring-mass model for the rolling stock. In order to solve the equations of motion, we applied the Runge-Kutta method, which is the typical numerical analysis method used for differential equations. Furthermore, the lateral displacement of the wheel-set at the wheel-rail interface was estimated using kinetic energy. In order to verify the derived equations, we compared the theoretical and simulated results under various collision conditions. The maximum relative deviations of the lateral displacements were 0.8 [%] ~ 4.7 [%] in light collisions and 0.6 [%] ~ 5.1 [%] under derailment conditions. When an accident is simulated, these theoretical equations can be used to predict the overall behavior and obtain the offset of the body-to-body link as the initial perturbation.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.25 no.4
• /
• pp.364-371
• /
• 2022

Surface state change of ammunition generated during the storage period increases the surface roughness and this affects the flight of ammunition, but there are no research results quantitatively indicating this. In this study, the drag force for each Mach number of howitzer shells was calculated through CFD to which the surface data of the howitzer shell was applied, and analysis of trajectory was performed using drag force values as an input of the 4th Runge-Kutta method, and the degree of decrease of the maximum range caused by the surface roughness of the howitzer shell was estimated. As a result, it was confirmed that the maximum range of howitzer shell with high surface roughness was 1.12 % shorter than that of howitzer shell without roughness. It was confirmed that the effect of surface roughness on the trajectory is not negligible.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.5 no.2
• /
• pp.86-96
• /
• 1995

In this paper the fuzzy extension for the classical engineering mechanics problems is studied. The governing differential equation is derived for the buckling loads of the columns with uncertain mediums: the their own weight and the flexural rigidity. The columns with one typical end constraint(hinged1 clarnped/free) and the other finite rotational spring with fuzzy constant are considered in numerical examples. The vertex method is used to evaluate the fuzzy functions. The Runge-Kutta method and Determinant Search method are used to solve the differential equation and determine the buckling loads, respectively. The membership functions of the buckling load are calculated. The index of fuzziness to quantitatively describe the propagation of fuzziness is defined. According to the fuzziness of governing factors, the varlation of index of fuzziness for buckling load is investigated, and the sensitivity for the end constraints is analyzed.

• Structural Engineering and Mechanics
• /
• v.50 no.4
• /
• pp.487-500
• /
• 2014

An analytical solution for the nonlinear in-plane free oscillations of the suspended cable which contains the quadratic and cubic nonlinearities is investigated via the homotopy analysis method (HAM). Different from the existing analytical technique, the HAM is indeed independent of the small parameter assumption in the nonlinear vibration equation. The nonlinear equation is established by using the extended Hamilton's principle, which takes into account the effects of the geometric nonlinearity and quasi-static stretching. A non-zero equilibrium position term is introduced due to the quadratic nonlinearity in order to guarantee the rule of the solution expression. Therefore, the mth-order analytic solutions of the corresponding equation are explicitly obtained via the HAM. Numerical results show that the approximate solutions obtained by using the HAM are in good agreement with the numerical integrations (i.e., Runge-Kutta method). Moreover, the HAM provides a simple way to adjust and control the convergent regions of the series solutions by means of an auxiliary parameter. Finally, the effects of initial conditions on the linear and nonlinear frequency ratio are investigated.

• Tribology and Lubricants
• /
• v.33 no.6
• /
• pp.309-314
• /
• 2017

This paper presents a theoretical investigation of the dynamic characteristics of externally pressurized air pad bearings with closed loop grooves. These grooves are made on the surface of bearings to reduce the number of supply holes so that manufacturing costs can be reduced. The semi-implicit method is applied to calculate the time varying pressure profile on the air bearing surface owing to the advantages of numerical stability and fast time tracing characteristics. The static pressure of the groove bearings is much higher than that without grooves, so the groove bearings can provide high load carrying capacity. The equation of motion considering vertical motion and tilting motion are also solved using the Runge-Kutta 4th order method. By combining the semi-implicit method and the Runge-Kutta method, fast calculations of the dynamic behavior of the air bearing can be achieved. The variations of bearing reaction force, air film reaction moment, height, and tilting angle are investigated for the step force input, which is 20% higher than the bearing reaction, when the nominal clearance is 6 mm. The effect of the groove width and the groove depth are investigated by calculating the dynamic behavior. The possibility of the air hammering with the depth of the groove is found and discussed.

• Journal of the Korean Magnetics Society
• /
• v.6 no.5
• /
• pp.281-286
• /
• 1996

For the accurate analysis of motional characteristics of electrical machines, it is needed to solve the motion equations together with the electromagnetic field equations. In this paper the sequential coupling of systems, the spring mass system and the electromagnetic system, is adopted. The induced current and the magnetic fields are calculated by finite element method(FEM) with given speed. And then, with the computed elec-tromagnetic force, the mechanical equations are solved by the Runge-Kutta method. The above two processes are repeated sequentially to obtain the time domain solutions. The resultant values are applied to the energy conservation law to prove the usefulness of the proposed sequential method.