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

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

진원도 오차를 고려한 스핀들 시스템의 동적 특성 해석 (Dynamic Analysis of a Tilted HDD Spindle System due to Roundness)

  • 곽규열;장건희
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2007년도 추계학술대회논문집
    • /
    • pp.840-846
    • /
    • 2007
  • This paper investigates the dynamic behavior of a HDD spindle system due to the imperfect roundness of a rotating shaft. The shaft of a spindle motor rotates with eccentricity by the unbalanced mass of the rotating part. The eccentricity generates the run-out of a spindle motor which results in the eccentric motion of a rotating part. Roundness of a shaft affects this motion which limits the memory capacity of a HDD. This research proposes a modified Reynolds equation for the coupled journal and thrust FDBs to include the variable film thickness due to the roundness. Finite element method is used to solve the Reynolds equation for the pressure distribution. Reaction forces and friction torque are obtained by integrating the pressure and shear stress, respectively. The dynamic behavior is determined by solving the equations of a motion of a HDD spindle system in six degrees of freedom with the Runge-Kutta method to characterize the motion of a rotating part. This research shows that the roundness of a rotating shaft causes the excitation frequency with integer multiple of a rotating frequency.

  • PDF

수치 파동 수조를 이용한 비선형파의 파형변화와 속도분포 해석 (Spatial Modulation of Nonlinear Waves and Their Kinematics using a Numerical Wave Tank)

  • 구원철;최가람
    • 한국해양공학회지
    • /
    • 제23권6호
    • /
    • pp.12-16
    • /
    • 2009
  • In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

가공 오차를 고려한 스핀들 시스템의 동적 특성 해석 (Dynamic Analysis of a Tilted HDD spindle system due to Manufacturing Tolerance)

  • 곽규열;김학운;장건희
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2007년도 춘계학술대회논문집
    • /
    • pp.852-858
    • /
    • 2007
  • This paper investigates the dynamic characteristics of a tilted HDD spindle system with fluid dynamic bearings (FDBs). Tilting motion of a HDD spindle system may be caused by improper manufacturing tolerance, such as imperfect cylindricity between shaft and sleeve of FDBs, imperfect perpendicularity between shaft and thrust as well as the gyroscopic moment of the unbalanced mass of the rotating part. Tilting motion may result in the instability of the HDD spindle system and it may increase the disk run-out to limit memory capacity. This research proposes a modified Reynolds equation for the coupled journal and thrust FDBs to include the variable film thickness due to the cylindricity and the perpendicularity. Finite element method is used to solve the Reynolds equation for the pressure distribution. Reaction forces and friction torque are obtained by integrating the pressure and shear stress, respectively. The dynamic behavior is determined by solving the equations of a motion of a HDD spindle system in six degrees of freedom with the Runge-Kutta method to study whirling and tilting motions. This research shows that the cylindricity and the perpendicularity increase the tilting angle and whirl radius of the rotor.

  • PDF

유압펌프용 실린더 블록의 윤활 및 동특성 해석 (Analysis of Lubrication and Dynamic Characteristics of a Cylinder Block for Hydraulic Pump)

  • 안성용;임윤철;홍예선
    • 한국윤활학회:학술대회논문집
    • /
    • 한국윤활학회 2004년도 학술대회지
    • /
    • pp.99-107
    • /
    • 2004
  • Lubrication characteristics between a cylinder block and a valve plate for high speed bent-axis type hydraulic pump play an important role in volumetric efficiency and durability of pump. In this paper, a finite element method is presented for the computation of the pressure distribution between a cylinder block and a valve plate for high speed bent-axis type hydraulic pump. Also, a Runge-Kutta method is applied to simulate the cylinder block dynamics of three-degrees of freedom motion. From the results of computation, we can draw two major conclusions. One is related to the fluid film characteristics between a cylinder block and a valve plate and the other is related to the average leakage that is determined by the pressure gradient and the clearance near the discharge port. The numerical results of cylinder block dynamics were compared with the experimental results using eddy-current type gap sensors those are fixed at a pump housing.

  • PDF

Large deflections of spatial variable-arc-length elastica under terminal forces

  • Phungpaingam, Boonchai;Athisakul, Chainarong;Chucheepsakul, Somchai
    • Structural Engineering and Mechanics
    • /
    • 제32권4호
    • /
    • pp.501-516
    • /
    • 2009
  • This paper aims to study the large deflections of variable-arc-length elastica subjected to the terminal forces (e.g., axial force and torque). Based on Kirchhoff's rod theory and with help of Euler parameters, the set of nonlinear governing differential equations which free from the effect of singularity are established together with boundary conditions. The system of nonlinear differential equations is solved by using the shooting method with high accuracy integrator, seventh-eighth order Runge-Kutta with adaptive step-size scheme. The error norm of end conditions is minimized within the prescribed tolerance ($10^{-5}$). The behavior of VAL elastica is studied by two processes. One is obtained by applying slackening first. After that keeping the slackening as a constant and then the twist angle is varied in subsequent order. The other process is performed by reversing the sequence of loading in the first process. The results are interpreted by observing the load-deflection diagram and the stability properties are predicted via fold rule. From the results, there are many interesting aspects such as snap-through phenomenon, secondary bifurcation point, loop formation, equilibrium configurations and effect of variable-arc-length to behavior of elastica.

Nonlinear vibration analysis of an electrostatically excited micro cantilever beam coated by viscoelastic layer with the aim of finding the modified configuration

  • Poloei, E.;Zamanian, M.;Hosseini, S.A.A.
    • Structural Engineering and Mechanics
    • /
    • 제61권2호
    • /
    • pp.193-207
    • /
    • 2017
  • In this study, the vibration of an electrostatically actuated micro cantilever beam is analyzed in which a viscoelastic layer covers a portion of the micro beam length. This proposed model is considered as the main element of mass and pollutant micro sensors. The nonlinear motion equation is extracted by means of Hamilton principle, considering nonlinear shortening effect for Euler-Bernoulli beam. The non-linear effects of electrostatic excitation, geometry and inertia have been taken into account. The viscoelastic model is assumed as Kelvin-Voigt model. The motion equation is discretized by Galerkin approach. The linear free vibration mode shapes of non-uniform micro beam i.e. the linear mode shape of the system by considering the geometric and inertia effects of viscoelastic layer, have been employed as comparison function in the process of the motion equation discretization. The discretized equation of motion is solved by the use of multiple scale method of perturbation theory and the results are compared with the results of numerical Runge-Kutta approach. The frequency response variations for different lengths and thicknesses of the viscoelastic layer have been founded. The results indicate that if a constant volume of viscoelastic layer is to be deposited on the micro beam for mass or gas sensor applications, then a modified configuration may be found by using the analysis of this paper.

Analytical, numerical and experimental investigation of low velocity impact response of laminated composite sandwich plates using extended high order sandwich panel theory

  • Salami, Sattar Jedari;Dariushi, Soheil
    • Structural Engineering and Mechanics
    • /
    • 제68권3호
    • /
    • pp.325-334
    • /
    • 2018
  • The Nonlinear dynamic response of a sandwich plate subjected to the low velocity impact is theoretically and experimentally investigated. The Hertz law between the impactor and the plate is taken into account. Using the Extended High Order Sandwich Panel Theory (EHSAPT) and the Ritz energy method, the governing equations are derived. The skins follow the Third order shear deformation theory (TSDT) that has hitherto not reported in conventional EHSAPT. Besides, the three dimensional elasticity is used for the core. The nonlinear Von Karman relations for strains of skins and the core are adopted. Time domain solution of such equations is extracted by means of the well-known fourth-order Runge-Kutta method. The effects of core-to-skin thickness ratio, initial velocity of the impactor, the impactor mass and position of the impactor are studied in detail. It is found that these parameters play significant role in the impact force and dynamic response of the sandwich plate. Finally, some low velocity impact tests have been carried out by Drop Hammer Testing Machine. The results are compared with experimental data acquired by impact testing on sandwich plates as well as the results of finite element simulation.

최적회된 고차-고해상도 집적 유한 차분법을 이용한 초음속 제트 스크리치 톤 수치 해석 (Numerical Analysis of Supersonic Axisymmetric Screech Tone Noise Using Optimized High-Order, High-Resolution Compact Scheme)

  • 이인철;이덕주
    • The Journal of the Acoustical Society of Korea
    • /
    • 제25권1E호
    • /
    • pp.32-35
    • /
    • 2006
  • The screech tone of underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. Fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. Adaptive nonlinear artificial dissipation model and generalized characteristic boundary condition are also used. The screech tone, generated by a closed loop between instability waves and quasi-periodic shock cells at the near field, is reasonably analyzed with present numerical methods for the underexpanded jet having Mach number 1.13. First of all, the centerline mean pressure distribution is calculated and compared with experimental and other numerical results. The instantaneous density contour plot shows Mach waves due to mixing layer convecting supersonically, which propagate downstream. The pressure signal and its Fourier transform at upstream and downstream shows the directivity pattern of screech tone very clearly. Most of all, we can simulate the axisymmetric mode change of screech tone very precisely with present method. It can be concluded that the basic phenomenon of screech tone including the frequency can be calculated by using high-order and high-resolution schemes without any specific numerical modeling for screech tone feedback loop.

Bragg 반사에 의한 비선형파의 공간적 파형변조해석 (Spatial Modulation of Nonlinear Waves due to Bragg Reflection)

  • 최가람;구원철
    • 한국해양공학회지
    • /
    • 제24권6호
    • /
    • pp.34-40
    • /
    • 2010
  • Bragg reflection of nonlinear waves is simulated by a 2D fully nonlinear numerical wave tank (NWT). The developed NWT was based on the Boundary Element Method (BEM) with potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme with Runge-Kutta 4th-order time integration. A spatial variation of wave elevations and their Fourier amplitudes of each component are compared to investigate the effect of sea bottom ripples and their relative heights. The incident waves over an undulated sea bottom are partially reflected and changed to partial standing waves due to Bragg reflection. The present results are verified with linear calculations and experimental data. It is found that the 1st-order wave component is mainly affected by Bragg reflection and its spatial modulation is significant in front of the bottom ripples.

Geometrically nonlinear analysis of sandwich beams under low velocity impact: analytical and experimental investigation

  • Salami, Sattar Jedari;Dariushi, Soheil
    • Steel and Composite Structures
    • /
    • 제27권3호
    • /
    • pp.273-283
    • /
    • 2018
  • Nonlinear low velocity impact response of sandwich beam with laminated composite face sheets and soft core is studied based on Extended High Order Sandwich Panel Theory (EHSAPT). The face sheets follow the Third order shear deformation beam theory (TSDT) that has hitherto not reported in conventional EHSAPT. Besides, the two dimensional elasticity is used for the core. The nonlinear Von Karman type relations for strains of face sheets and the core are adopted. Contact force between the impactor and the beam is obtained using the modified Hertz law. The field equations are derived via the Ritz based applied to the total energy of the system. The solution is obtained in the time domain by implementing the well-known Runge-Kutta method. The effects of boundary conditions, core-to-face sheet thickness ratio, initial velocity of the impactor, the impactor mass and position of the impactor are studied in detail. It is found that each of these parameters have significant effect on the impact characteristics which should be considered. Finally, some low velocity impact tests have been carried out by Drop Hammer Testing Machine. The contact force histories predicted by EHSAPT are in good agreement with that obtained by experimental results.