• Title/Summary/Keyword: Runge-kutta method

Search Result 506, Processing Time 0.029 seconds

Dynamic Behavior of Rotating Shaft System Corresponding to Operating Modes (운전모드에 따른 회전축계의 동적거동)

  • Kim, Sang-Hwan
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.20 no.9
    • /
    • pp.2744-2751
    • /
    • 1996
  • In case of limited power supply, a rotating shaft system may not reach its operating speed that is greater than its critical speed, but the speed oscillates with small ampllitude near critical speed. As a result, it is considered that the operating mode plays an important role in the smooth start of machines. In order to investigate the dynamic behaviors of the rotating shaft system at the beginning stage, one has derived the equations of motion whose degrees of freedom is three, two translations and one rotation. The simultaneous differential equations are numerically solved by using runge-Kutta method, and thus the small time step length could be required corresponding to the stability of solution. Three types of operating modes dependent upon the driving torque rate have been numerically investigated according to the maximum displacement of shaft center. The first type of relation is linear, the second type is composed of two linear curves recommended by machine manufacturer, and the last one is the proposed torque curve reflecting the frequency response curve of one degree of freedom system. For the second type of modes, it is found that the optimal range of intermediate speed to the critical speed lies between 0.8 and 0.9. In addition to that, the maximum displacement can be reduced more if the third type of mode is utilized.

Effect of Geometric Parameters on the Performance of an Automotive Scroll Compressor Using R-134a (R-134a를 사용한 자동차용 스크롤 압축기의 스크롤 형상변화가 성능에 미치는 영향)

  • Lee, Geonho;Kim, Haksoo;Cho, Keumnam;Yoo, Jungyul
    • Transactions of the Korean Society of Mechanical Engineers B
    • /
    • v.22 no.11
    • /
    • pp.1635-1646
    • /
    • 1998
  • The present study investigated the effects of geometric parameters considered on the performance of an automotive scroll compressor by assuming ideal, semi-real and real gases for R-134a. The geometric parameters were center thickness of scroll, height of scroll and the size of discharge port. Fourth-order Runge-Kutta method was applied to solve the thermodynamic equations, leaking rate equation and the equation of motion of discharge valve for ideal, semi-real, and real gases. The volumetric and adiabatic efficiencies for semi-real and real gases differed little, but those for ideal gas differed by 18% and 25% compared with those for real gas at 2,000rpm. The volumetric efficiency changed little as the design angle of scroll (${\gamma}$) changed, but the adiabatic efficiency at ${\gamma}$ of $34^{\circ}$ was higher by 2.4% than that of $147^{\circ}$ for 2,000rpm. The volumetric and adiabatic efficiencies at scroll height of 29.8mm were higher by 1.7% and 2.8% than those of 65.8mm. The volumetric efficiency changed little as the size of discharge port changed, but the adiabatic efficiency increased a little as the size of discharge port decreased.

Non-linearity Characteristics Consideration of the Disc Spring Under Shock Using Numerical Analysis (수치해석을 통한 충격 시 접시 스프링의 비선형 특성 고찰)

  • Bang, Seung-Woo;Lee, Hae-Jin;Sim, Hyun-Jin;Park, Sang-Gil;Lee, Jung-Youn;Oh, Jae-Eung
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.18 no.2
    • /
    • pp.153-159
    • /
    • 2008
  • General coil spring has linearity. However, disc spring has non-linearity so that using this non-linearity disc spring can be designed to do shock-absorbing in cases we need because shock response also has non-linearity. By changing the shape and stacking number, it is satisfactory with response of displacement, velocity and acceleration. Conventionally, disc spring was used to control the vibration against huge load and limited space. However, it is limitedly used because of difficulty of the designing guidance. Therefore, disc spring is needed to study further in order to apply it widely. Response of disc spring is compared to response of coil spring by changing $h_o/t$ ratio with computer simulation and the usage of disc spring is increased through analysis of effect of design factors. The purpose of this paper is that the shock response of disc spring is calculated through numerical simulation and effect of $h_o/t$ and stiffness is analyzed to broad usage so that design factor of disc spring is presented.

Comparison of Acceleration of Vibration Isolator and Accelerometer Attached Vibration Isolator Using Numerical Analysis (수치 해석을 이용하여 제진대와 제진대에 부착된 가속도계의 가속도 비교)

  • Shin, Dong Ho;Lee, Jung Woo;Oh, Jae-Eung;Lee, Jung Youn
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.23 no.1
    • /
    • pp.17-24
    • /
    • 2013
  • The process of producing high precision and light weight product is always exposed to impact load or shock. Because of this, isolator device is required. To measure the response of the isolator, accelerometer is practically used. However, the measured response of the accelerometer is different to the response of the isolator. To predict the response of the accelerometer and the isolator, 2-DOF damped system with an input shock is modeled using numerical analysis. 1-DOF damped system with a base excitation is also used to predict the response of the isolator. The mass ratio, damping ratio, and natural frequency ratio are then varied. The predicted responses from the two modeling approaches are compared and large errors are found.

Comparison of Three Optimization Methods Using Korean Population Data

  • Oh, Deok-Kyo
    • Korean System Dynamics Review
    • /
    • v.13 no.2
    • /
    • pp.47-71
    • /
    • 2012
  • The purpose of this research is the examination of validity of data as well as simulation model, i.e. to simulate the real data in the SD model with the least error using the adjustments for the faithful reflection of real data to the simulation. In general, SD programs (e.g. VENSIM) utilize the Euler or Runge-Kutta method as an algorithm. It is possible to reflect the trend of real data via these two estimation methods however can cause the validity problem in case of the simulation requiring the accuracy as they have endogenous errors. In this article, the future population estimated by the Korea National Statistical Office (KNSO) to 2050 is simulated by the aging chain model, dividing the population into three cohorts, 0-14, 15-64, 65 and over cohorts by age and offering the adjustments to them. Adjustments are calculated by optimization with three different methods, optimization in EXCEL, manual optimization with iterative calculation, and optimization in VENSIM DSS, the results are compared, and at last the optimal adjustment set with the least error are found among them. The simulation results with the pre-determined optimal adjustment set are validated by methods proposed by Barlas (1996) and other alternative methods. It is concluded that the result of simulation model in this research has no significant difference from the real data and reflects the real trend faithfully.

  • PDF

Elasto-Plastic Finite Element Analysis in Consideration of Phase Transformations (상변태를 고려한 탄소성 유한요소 해석)

  • Lee, M.G.;Kim, S.J.;Jeong, W.C.
    • Proceedings of the Korean Society for Technology of Plasticity Conference
    • /
    • 2009.05a
    • /
    • pp.334-336
    • /
    • 2009
  • An elastic-plasticity model during the austenitic decomposition was derived and implemented to incorporate the two important deformation behaviors observed during the phase transformations: the volumetric strain and transformation induced plasticity due to the temperature change and phase transformation. To obtain transformed phase volume fractions during cooling, the fourth order Runge-Kutta method was used to solve the Kirkaldy's phase kinetics model which is function of temperature, austenitic grain size and chemical composition. The volumetric strain was calculated by considering the densities of constituent phases, while the transformation induced plasticity was based on the micro-plasticity due to the volume mismatch between soft austenitic phase and other harder phases. The constitutive equations were implemented into the implicit finite element software and a simple boundary value problem was chosen as a model problem to validate the effect of transformation plasticity on the deformation behavior of steel under cooling from high temperature. It was preliminary concluded that the transformation plasticity plays a critical role in relaxing the developed stress during forming and thus reducing the magnitude of springback.

  • PDF

Three-Dimensional Trajectory of a Fluid Particle in Air with Wind Effects and Air Resistance (공기 저항과 바람의 영향을 고려한 대기에서의 유체입자의 3차원 궤적)

  • 이동렬
    • Journal of Advanced Marine Engineering and Technology
    • /
    • v.25 no.4
    • /
    • pp.797-808
    • /
    • 2001
  • Three-dimensional trajectory of fluid particle is simulated by a particle motion, which is able to examine the influences of changes in the several parameters. To calculate the trajectory of a particle, the Runge-Kutta method was utilized. The use of a projectile of particles for the trajectory of liquid jet has been shown to be useful to estimate the influence of different operating parameters such as best particle diameter, density of liquid body, initial take-off velocity, wind velocity, cross wind velocity, take-off angle, and base angle for a released flow from the nozzle. The results give the trajectories of various types of particle of body and at different elevations, base angles, wind velocities and densities of liquid body. The trajectories in a vacuum show that air resistances decreases both the distance and the maximum height of a projectile, and also explain that the termination time is also reduced in air. In addition, the maximum distance in the x direction was obtained with take-off angles from 30 degrees to 45 degrees in still air and the projectile of particles was highly effected by wind and cross wind. Clearly, a particle has to be so positioned as to take the optimum possible advantage of the wind if the maximum distances is requested. The wind astern increased the maximum distances of x direction compared with the wind ahead. Finally, it is possible to optimize the design of pump by using these results.

  • PDF

A Simulation Study of Impedance Plethysmography for Diagnosing Deep Vein Thrombosis (Deep Vein Thrombosis 진단을 위한 Impedance Plethysmography의 시뮬레이션 연구)

  • Lee, Jeon;Lee, Kyoung-Joung
    • The Transactions of the Korean Institute of Electrical Engineers D
    • /
    • v.50 no.10
    • /
    • pp.494-501
    • /
    • 2001
  • In this study, the effects of vascular parameter changes and electrodes on VOP measurement based on IPG were simulated mathematically. For the evaluation of the effects of hemodynamic changes on VOP, a mathematical model, which consists of cardiovascular system model and venous occlusion model, was developed and the model solution representing the blood flow and pressure in measuring point was found by 2nd order Runge-Kutta method. And, with sensitivity coefficients obtained from finite element solution of electric field in measuring point, the effects of electrode system on measurement were evaluated. As increasing the resistance, the venous capacitance was not changed but the venous outflows were decreased and the decreased compliance reduced the venous capacitance. And, for several configurations of round electrodes and band electrodes, the sensitivity coefficients were computed using the electric field distribution along deep vein. In conclusion, the proposed mathematical cardiovascular model could be applied to the simulation study on the effects of hemodynamic parameters on DVT diagnosis with IPG. And, also the sensitivity coefficients could provide effective electrode configuration for exact measurement of VOP.

  • PDF

System Representation for the Control System of the Follow-up Mechanism on the Marine Gyro Compass

  • Sang-Jib Lee
    • Journal of the Korean Institute of Navigation
    • /
    • v.4 no.1
    • /
    • pp.31-50
    • /
    • 1980
  • It does not seem necessarily practicable to keep the system always in optimal condition, athough the control system of the follow-up mechanism on the most marine gyro compasses is to be adjusted by the operator through the gain adjustment. Sometimes a sustained oscillation or an incorrect gyro reading occurs to the system. For such a system any systematical research or theoretical basis of the guide for the optimal gain adjustment has not been reported yet. As a basic investigation of the theoretical system analysis to solve the problems concerned, the author attempts in this paper to express the system in a mathematical model deduced from the results of the theoretical approach and the experimental observation of each element contained in the follow-up mechanism of Hokshin D-1 gyro compass, and to constitute an over-all closed loop transfer function. This funciton being reverted to a fourth orderlinear differential equation, the first order simultaneous differential equations are obtained by means of the state-variables. The latter equations are solved by the Runge-Kutta method with digital computer. By comparing the characteristic of the simulated over-all output with that of the experimental result, it is shown that both outputs are nearly consistent with each other. It is also expected that the system representation proposed by this paper is valid and will be a prospective means in a further study on the design and optimal adjustment of the system.

  • PDF

Optimization of the Earthing Resistance and Research on the Electrical Characteristics of New Catalyst for the Quality Earthing (접지품질 개선을 위한 접지저항 최적화 및 접지충진제의 전기적 특성 연구)

  • Park, C.G.;Ahn, S.J.;Woo, J.W.;Ahn, S.J.;Yu, Y.J.;Ahn, S.J.
    • Journal of the Korea Academia-Industrial cooperation Society
    • /
    • v.7 no.6
    • /
    • pp.1078-1085
    • /
    • 2006
  • In earthing technologies, the intrusion of the surge wave can be protected when the earthing resistance is as low as several ohms. However, the cost fer realizing such a low resistance is quite high. Therefore it is important to determine appropriate target value of the earthing resistance to install cost-effective grounding system. In this work, we have proposed an optimization method of the earthing resistance according to the various environmental parameters by using a numerical simulator We have also developed new catalyst composed of water-soluble alkaline elements and activated carbon and investigated its electrical characteristics for enhancement of the earthing qualities.

  • PDF