• Title/Summary/Keyword: Runge Kutta Method

Search Result 502, Processing Time 0.027 seconds

Dynamic Snapping and Frequency Characteristics of 3-Free-Nodes Spatial Truss Under the Periodic Loads (주기 하중을 받는 3-자유절점 공간 트러스의 동적 불안정 현상과 주파수 특성)

  • Shon, Sudeok;Hwang, Kyung-Ju
    • Journal of Korean Association for Spatial Structures
    • /
    • v.20 no.4
    • /
    • pp.149-158
    • /
    • 2020
  • The governing equation for a dome-type shallow spatial truss subjected to a transverse load is expressed in the form of the Duffing equation, and it can be derived by considering geometrical non-linearity. When this model under constant load exceeds the critical level, unstable behavior is appeared. This phenomenon changes sensitively as the number of free-nodes increases or depends on the imperfection of the system. When the load is a periodic function, more complex behavior and low critical levels can be expected. Thus, the dynamic unstable behavior and the change in the critical point of the 3-free-nodes space truss system were analyzed in this work. The 4-th order Runge-Kutta method was used in the system analysis, while the change in the frequency domain was analyzed through FFT. The sinusoidal wave and the beating wave were utilized as the periodic load function. This unstable situation was observed by the case when all nodes had same load vector as well as by the case that the load vector had slight difference. The results showed the critical buckling level of the periodic load was lower than that of the constant load. The value is greatly influenced by the period of the load, while a lower critical point was observed when it was closer to the natural frequency in the case of a linear system. The beating wave, which is attributed to the interference of the two frequencies, exhibits slightly more behavior than the sinusoidal wave. And the changing of critical level could be observed even with slight changes in the load vector.

Two-dimensional curved panel vibration and flutter analysis in the frequency and time domain under thermal and in-plane load

  • Moosazadeh, Hamid;Mohammadi, Mohammad M.
    • Advances in aircraft and spacecraft science
    • /
    • v.8 no.4
    • /
    • pp.345-372
    • /
    • 2021
  • The analysis of nonlinear vibrations, buckling, post-buckling, flutter boundary determination and post-flutter behavior of a homogeneous curved plate assuming cylindrical bending is conducted in this article. Other assumptions include simply-supported boundary conditions, supersonic aerodynamic flow at the top of the plate, constant pressure conditions below the plate, non-viscous flow model (using first- and third-order piston theory), nonlinear structural model with large deformations, and application of mechanical and thermal loads on the curved plate. The analysis is performed with constant environmental indicators (flow density, heat, Reynolds number and Mach number). The material properties (i.e., coefficient of thermal expansion and modulus of elasticity) are temperature-dependent. The equations are derived using the principle of virtual displacement. Furthermore, based on the definitions of virtual work, the potential and kinetic energy of the final relations in the integral form, and the governing nonlinear differential equations are obtained after fractional integration. This problem is solved using two approaches. The frequency analysis and flutter are studied in the first approach by transferring the handle of ordinary differential equations to the state space, calculating the system Jacobin matrix and analyzing the eigenvalue to determine the instability conditions. The second approach discusses the nonlinear frequency analysis and nonlinear flutter using the semi-analytical solution of governing differential equations based on the weighted residual method. The partial differential equations are converted to ordinary differential equations, after which they are solved based on the Runge-Kutta fourth- and fifth-order methods. The comparison between the results of frequency and flutter analysis of curved plate is linearly and nonlinearly performed for the first time. The results show that the plate curvature has a profound impact on the instability boundary of the plate under supersonic aerodynamic loading. The flutter boundary decreases with growing thermal load and increases with growing curvature.

Development and Application of Two-Dimensional Numerical Tank using Desingularized Indirect Boundary Integral Equation Method (비특이화 간접경계적분방정식방법을 이용한 2차원 수치수조 개발 및 적용)

  • Oh, Seunghoon;Cho, Seok-kyu;Jung, Dongho;Sung, Hong Gun
    • Journal of Ocean Engineering and Technology
    • /
    • v.32 no.6
    • /
    • pp.447-457
    • /
    • 2018
  • In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

Implementation of Markerless Augmented Reality with Deformable Object Simulation (변형물체 시뮬레이션을 활용한 비 마커기반 증강현실 시스템 구현)

  • Sung, Nak-Jun;Choi, Yoo-Joo;Hong, Min
    • Journal of Internet Computing and Services
    • /
    • v.17 no.4
    • /
    • pp.35-42
    • /
    • 2016
  • Recently many researches have been focused on the use of the markerless augmented reality system using face, foot, and hand of user's body to alleviate many disadvantages of the marker based augmented reality system. In addition, most existing augmented reality systems have been utilized rigid objects since they just desire to insert and to basic interaction with virtual object in the augmented reality system. In this paper, unlike restricted marker based augmented reality system with rigid objects that is based in display, we designed and implemented the markerless augmented reality system using deformable objects to apply various fields for interactive situations with a user. Generally, deformable objects can be implemented with mass-spring modeling and the finite element modeling. Mass-spring model can provide a real time simulation and finite element model can achieve more accurate simulation result in physical and mathematical view. In this paper, the proposed markerless augmented reality system utilize the mass-spring model using tetraheadron structure to provide real-time simulation result. To provide plausible simulated interaction result with deformable objects, the proposed method detects and tracks users hand with Kinect SDK and calculates the external force which is applied to the object on hand based on the position change of hand. Based on these force, 4th order Runge-Kutta Integration is applied to compute the next position of the deformable object. In addition, to prevent the generation of excessive external force by hand movement that can provide the natural behavior of deformable object, we set up the threshold value and applied this value when the hand movement is over this threshold. Each experimental test has been repeated 5 times and we analyzed the experimental result based on the computational cost of simulation. We believe that the proposed markerless augmented reality system with deformable objects can overcome the weakness of traditional marker based augmented reality system with rigid object that are not suitable to apply to other various fields including healthcare and education area.

Image Warping Using Vector Field Based Deformation and Its Application to Texture Mapping (벡터장 기반 변형기술을 이용한 이미지 와핑 방법 : 텍스쳐 매핑에의 응용을 중심으로)

  • Seo, Hye-Won;Cordier, Frederic
    • Journal of KIISE:Computer Systems and Theory
    • /
    • v.36 no.5
    • /
    • pp.404-411
    • /
    • 2009
  • We introduce in this paper a new method for smooth foldover-free warping of images, based on the vector field deformation technique proposed by Von Funck et al. It allows users to specify the constraints in two different ways: positional constraints to constrain the position of a point in the image and gradient constraints to constrain the orientation and scaling of some parts of the image. From the user-specified constraints, it computes in the image domain a C1-continuous velocity vector field, along which each pixel progressively moves from its original position to the target. The target positions of the pixels are obtained by solving a set of partial derivative equations with the 4th order Runge-Kutta method. We show how our method can be useful for texture mapping with hard constraints. We start with an unconstrained planar embedding of a target mesh using a previously known method (Least Squares Conformal Map). Then, in order to obtain a texture map that satisfies the given constraints, we use the proposed warping method to align the features of the texture image with those on the unconstrained embedding. Compared to previous work, our method generates a smoother texture mapping, offers higher level of control for defining the constraints, and is simpler to implement.

Dynamic Instability and Multi-step Taylor Series Analysis for Space Truss System under Step Excitation (스텝 하중을 받는 공간 트러스 시스템의 멀티스텝 테일러 급수 해석과 동적 불안정)

  • Lee, Seung-Jae;Shon, Su-Deok
    • Journal of Korean Society of Steel Construction
    • /
    • v.24 no.3
    • /
    • pp.289-299
    • /
    • 2012
  • The goal of this paper is to apply the multi-step Taylor method to a space truss, a non-linear discrete dynamic system, and analyze the non-linear dynamic response and unstable behavior of the structures. The accurate solution based on an analytical approach is needed to deal with the inverse problem, or the dynamic instability of a space truss, because the governing equation has geometrical non-linearity. Therefore, the governing motion equations of the space truss were formulated by considering non-linearity, where an accurate analytical solution could be obtained using the Taylor method. To verify the accuracy of the applied method, an SDOF model was adopted, and the analysis using the Taylor method was compared with the result of the 4th order Runge-Kutta method. Moreover, the dynamic instability and buckling characteristics of the adopted model under step excitation was investigated. The result of the comparison between the two methods of analysis was well matched, and the investigation shows that the dynamic response and the attractors in the phase space can also delineate dynamic snapping under step excitation, and damping affects the displacement of the truss. The analysis shows that dynamic buckling occurs at approximately 77% and 83% of the static buckling in the undamped and damped systems, respectively.

Integrated Algorithm for Identification of Long Range Artillery Type and Impact Point Prediction With IMM Filter (IMM 필터를 이용한 장사정포의 탄종 분리 및 탄착점 예측 통합 알고리즘)

  • Jung, Cheol-Goo;Lee, Chang-Hun;Tahk, Min-Jea;Yoo, Dong-Gil;Sohn, Sung-Hwan
    • Journal of the Korean Society for Aeronautical & Space Sciences
    • /
    • v.50 no.8
    • /
    • pp.531-540
    • /
    • 2022
  • In this paper, we present an algorithm that identifies artillery type and rapidly predicts the impact point based on the IMM filter. The ballistic trajectory equation is used as a system model, and three models with different ballistic coefficient values are used. Acceleration was divided into three components of gravity, air resistance, and lift. And lift acceleration was added as a new state variable. The kinematic condition that the velocity vector and lift acceleration are perpendicular was used as a pseudo-measurement value. The impact point was predicted based on the state variable estimated through the IMM filter and the ballistic coefficient of the model with the highest mode probability. Instead of the commonly used Runge-Kutta numerical integration for impact point prediction, a semi-analytic method was used to predict impact point with a small amount of calculation. Finally, a state variable initialization method using the least-square method was proposed. An integrated algorithm including artillery type identification, impact point prediction and initialization was presented, and the validity of the proposed method was verified through simulation.

A Comic Facial Expression Method for Intelligent Avatar Communications in the Internet Cyberspace (인터넷 가상공간에서 지적 아바타 통신을 위한 코믹한 얼굴 표정의 생성법)

  • 이용후;김상운;청목유직
    • Journal of the Institute of Electronics Engineers of Korea CI
    • /
    • v.40 no.1
    • /
    • pp.59-73
    • /
    • 2003
  • As a means of overcoming the linguistic barrier between different languages in the Internet, a new sign-language communication system with CG animation techniques has been developed and proposed. In the system, the joint angles of the arms and the hands corresponding to the gesture as a non-verbal communication tool have been considered. The emotional expression, however, could as play also an important role in communicating each other. Especially, a comic expression is more efficient than real facial expression, and the movements of the cheeks and the jaws are more important AU's than those of the eyebrow, eye, mouth etc. Therefore, in this paper, we designed a 3D emotion editor using 2D model, and we extract AU's (called as PAU, here) which play a principal function in expressing emotions. We also proposed a method of generating the universal emotional expression with Avatar models which have different vertex structures. Here, we employed a method of dynamically adjusting the AU movements according to emotional intensities. The proposed system is implemented with Visual C++ and Open Inventor on windows platforms. Experimental results show a possibility that the system could be used as a non-verbal communication means to overcome the linguistic barrier.

Development of Runoff Hydrograph Model for the Derivation of Optimal Design Flood of Agricultural Hydraulic Structures(1) (농업수리구조물의 적정설계홍수량 유도를 위한 유출수문곡선모형의 개발(I))

  • 이순혁;박명근;맹승진
    • Magazine of the Korean Society of Agricultural Engineers
    • /
    • v.37 no.3_4
    • /
    • pp.34-47
    • /
    • 1995
  • It is experienced fact as a regular annual event that the structure to he designed on unreasonable flood for the agricultural structures including reservoirs have been brought not only loss of lives, but also enormous property damage. For the solution of this problem at issue, this study was conducted to develop an optimal runoff hydrograph model by comparison of the peak flows and time to peak between observed and simulated flows derived by linear time-invariant and linear time-variant models under the condition of having a short duration of heavy rainfall with uniform rainfall intensity at nine small watersheds which are within the range of 55.9 to 140.7 square kilometers in area in Han, Geum, Nagdong and Yeongsan Rivers. The results obtained through this study can be summarized as follows. 1. Storage constants and Gamma function arguments were calculated within the range of 1.2 to 6.42 and of 1.28 to 8.05 respectively by the moment method as the parameters for the analysis of runoff hydrograph based on linear time-invariant model. 2. Parameters for both linear time-invariant and linear time-variant models were calibrated with nine gaged watershed data, using a trial and error method. The resulting parameters including Gamma function argument, N and storage constant, K for linear time-invariant model were related statistically to watershed characteristic variables such as area, slope, length of main stream and the centroid length of the basin. 3. Average relative errors of the simulated peak discharge of calibrated runoff hydrographs by using linear time-variant and linear time-invariant models were shown to be 0.75 and 5.42 percent respectively to the peak of observed runoff hydrographs. Correlation coefficients for the statistical analysis in the same condition were shown to be 0.999 and 0.978 with a high significance respectively. Therefore, it can be concluded that the accuracy of a linear time-variant model is approaching more closely to the observed runoff hydrograph than that of a linear time-invariant model in the applied watersheds. 4. Average relative errors of the time to peak of calibrated runoff hydrographs by using linear time-variant and linear time-invariant models were shown to be 16.44 and 19.89 percent respectively to the time to peak of observed runoff hydrographs. Correlation coefficients in the same condition were also shown to be 0.999 and 0.886 with a high significance respectively. 5. It can be seen that the shape of simulated hydrograph based on a linear time- variant model is getting closer to the observed runoff hydrograph than that of a linear time-invariant model in the applied watersheds. 6. Two different models were verified with different rainfall-runoff events from data for the calibration by relative error and correlation analysis. Consequently, it can be generally concluded that verification results for the peak discharge and time to peak of simulated runoff hydrographs were in good agreement with those of calibrated runoff hydrographs.

  • PDF

Rotordynamic Analysis of a Dual-Spool Turbofan Engine with Focus on Blade Defect Events (블레이드 손상에 따른 이축식 터보팬 엔진의 동적 안정성 해석)

  • Kim, Sitae;Jung, Kihyun;Lee, Junho;Park, Kihyun;Yang, Kwangjin
    • Tribology and Lubricants
    • /
    • v.36 no.2
    • /
    • pp.105-115
    • /
    • 2020
  • This paper presents a numerical study on the rotordynamic analysis of a dual-spool turbofan engine in the context of blade defect events. The blades of an axial-type aeroengine are typically well aligned during the compressor and turbine stages. However, they are sometimes exposed to damage, partially or entirely, for several operational reasons, such as cracks due to foreign objects, burns from the combustion gas, and corrosion due to oxygen in the air. Herein, we designed a dual-spool rotor using the commercial 3D modeling software CATIA to simulate blade defects in the turbofan engine. We utilized the rotordynamic parameters to create two finite element Euler-Bernoulli beam models connected by means of an inter-rotor bearing. We then applied the unbalanced forces induced by the mass eccentricities of the blades to the following selected scenarios: 1) fully balanced, 2) crack in the low-pressure compressor (LPC) and high pressure compressor (HPC), 3) burn on the high-pressure turbine (HPT) and low pressure compressor, 4) corrosion of the LPC, and 5) corrosion of the HPC. Additionally, we obtained the transient and steady-state responses of the overall rotor nodes using the Runge-Kutta numerical integration method, and employed model reduction techniques such as component mode synthesis to enhance the computational efficiency of the process. The simulation results indicate that the high-vibration status of the rotor commences beyond 10,000 rpm, which is identified as the first critical speed of the lower speed rotor. Moreover, we monitored the unbalanced stages near the inter-rotor bearing, which prominently influences the overall rotordynamic status, and the corrosion of the HPC to prevent further instability. The high-speed range operation (>13,000 rpm) coupled with HPC/HPT blade defects possibly presents a rotor-case contact problem that can lead to catastrophic failure.