• 제목/요약/키워드: Runge Kutta Method

검색결과 502건 처리시간 0.024초

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

  • 신동호;이정우;오재응;이정윤
    • 한국소음진동공학회논문집
    • /
    • 제23권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

  • 오덕교
    • 한국시스템다이내믹스연구
    • /
    • 제13권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)

  • 이명규;김성준;정우창
    • 한국소성가공학회:학술대회논문집
    • /
    • 한국소성가공학회 2009년도 춘계학술대회 논문집
    • /
    • 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

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

  • 이동렬
    • Journal of Advanced Marine Engineering and Technology
    • /
    • 제25권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

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

  • 이전;이경중
    • 대한전기학회논문지:시스템및제어부문D
    • /
    • 제50권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
    • 한국항해학회지
    • /
    • 제4권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)

  • 박철근;안성준;우제욱;안수준;유용정;안승준
    • 한국산학기술학회논문지
    • /
    • 제7권6호
    • /
    • pp.1078-1085
    • /
    • 2006
  • 시설물이나 장비로 surge wave가 유입되는 것을 방지하기 위해서는 접지시스템의 접지저항을 수 $\Omega$ 정도로 낮추어야 하지만 시공비를 낮추기 위하여 적절한 접지저항 값을 결정하는 것이 중요한 문제이다. 본 연구에서는 접지저항 simulator를 이용하여 최적화된 경제적 접지저항(임피던스) 값을 설정하였으며 수용성 알칼리족 원소와 활성탄소를 주성분으로 하는 새로운 접지충진제를 개발하여 접지품질을 개선하기 위한 전기적 특성을 평가하였다.

  • PDF

Nonlinear model to predict the torsional response of U-shaped thin-walled RC members

  • Chen, Shenggang;Ye, Yinghua;Guo, Quanquan;Cheng, Shaohong;Diao, Bo
    • Structural Engineering and Mechanics
    • /
    • 제60권6호
    • /
    • pp.1039-1061
    • /
    • 2016
  • Based on Vlasov's torsional theory of open thin-walled members and the nonlinear constitutive relations of materials, a nonlinear analysis model to predict response of open thin-walled RC members subjected to pure torsion is proposed in the current study. The variation of the circulatory torsional stiffness and warping torsional stiffness over the entire loading process and the impact of warping shear deformation on the torsion-induced rotation of the member are considered in the formulation. The torque equilibrium differential equation is then solved by Runge-Kutta method. The proposed nonlinear model is then applied to predict the behavior of five U-shaped thin-walled RC members under pure torsion. Four of them were tested in an earlier experimental study by the authors and the testing data of the fifth one were reported in an existing literature. Results show that the analytical predictions based on the proposed model agree well with the experimental data of all five specimens. This clearly shows the validity of the proposed nonlinear model analyzing behavior of U-shaped thin-walled RC members under pure torsion.

Non-stationary vibration and super-harmonic resonances of nonlinear viscoelastic nano-resonators

  • Ajri, Masoud;Rastgoo, Abbas;Fakhrabadi, Mir Masoud Seyyed
    • Structural Engineering and Mechanics
    • /
    • 제70권5호
    • /
    • pp.623-637
    • /
    • 2019
  • This paper analyzes the non-stationary vibration and super-harmonic resonances in nonlinear dynamic motion of viscoelastic nano-resonators. For this purpose, a new coupled size-dependent model is developed for a plate-shape nano-resonator made of nonlinear viscoelastic material based on modified coupled stress theory. The virtual work induced by viscous forces obtained in the framework of the Leaderman integral for the size-independent and size-dependent stress tensors. With incorporating the size-dependent potential energy, kinetic energy, and an external excitation force work based on Hamilton's principle, the viscous work equation is balanced. The resulting size-dependent viscoelastically coupled equations are solved using the expansion theory, Galerkin method and the fourth-order Runge-Kutta technique. The Hilbert-Huang transform is performed to examine the effects of the viscoelastic parameter and initial excitation values on the nanosystem free vibration. Furthermore, the secondary resonance due to the super-harmonic motions are examined in the form of frequency response, force response, Poincare map, phase portrait and fast Fourier transforms. The results show that the vibration of viscoelastic nanosystem is non-stationary at higher excitation values unlike the elastic ones. In addition, ignoring the small-size effects shifts the secondary resonance, significantly.

수중 폭발 시뮬레이션을 위한 경계 요소법 기반의 수치 해석 기법 연구 (A Study on BEM-Based Numerical Simulation Technique for Underwater Explosions)

  • 정준모;이재빈
    • 대한조선학회논문집
    • /
    • 제57권5호
    • /
    • pp.271-277
    • /
    • 2020
  • Recoverability and vulnerability of navy ships under underwater explosion are critical verification factors in the acquisition phase of navy ships. This paper aims to establish numerical analysis techniques for the underwater explosion of navy ships. Doubly Asymptotic Approach (DAA) Equation of Motion (EOM) of primary shock wave and secondary bubble pulse proposed by Geers-Hunter was introduced. Assuming a non-compressive fluid, reference solution of the DAA EOM of Geers-Hunter using Runge-Kutta method was derived for the secondary bubble pulse phase with an assumed charge conditions. Convergence analyses to determine fluid element size were performed, suggesting that the minimum fluid element size for underwater explosion analysis was 0.1 m. The spherical and cylindrical fluid domains were found to be appropriate for the underwater explosion analyses from the fluid domain shape study. Because the element size of 0.1 m was too small to be applied to the actual navy ships, a very slender beam with the square solid section was selected for the study of fluid domain existence effect. The two underwater explosion models with/without fluid domain provided very similar results in terms of the displacement and stress processes.