• 제목/요약/키워드: Runge-kutta method

Search Result 506, Processing Time 0.029 seconds

Active control of a nonlinear and hysteretic building structure with time delay

  • Liu, Kun;Chen, Long-Xiang;Cai, Guo-Ping
    • Structural Engineering and Mechanics
    • /
    • v.40 no.3
    • /
    • pp.431-451
    • /
    • 2011
  • Time delay inevitably exists in active control systems, and it may cause the degradation of control efficiency or instability of the systems. So time delay needs to be compensated in control design in order to eliminate its negative effect on control efficiency. Today time delay in linear systems has been more studied and some treating methods had been worked out. However, there are few treating methods for time delay in nonlinear systems. In this paper, an active controller for a nonlinear and hysteretic building structure with time delay is studied. The nonlinear and hysteretic behavior of the system is illustrated by the Bouc-Wen model. By specific transformation and augmentation of state parameters, the motion equation of the system with explicit time delay is transformed into the standard state space representation without any explicit time delay. Then the fourth-order Runge-Kutta method and instantaneous optimal control method are applied to the controller design with time delay. Finally, numerical simulations and comparisons of an eight-story building using the proposed time-delay controller are carried out. Simulation results indicate that the control performance will deteriorate if time delay is not taken into account in the control design. The simulations also prove the proposed time delay controller in this paper can not only effectively compensate time delay to get better control effectiveness, but also work well with both small and large time delay problems.

Control of Plume Interference Effects on a Missile Body Using a Porous Extension (다공확장벽을 이용한 미사일 동체에 대한 플룸간섭 현상의 제어)

  • Young-Ki Lee;Heuy-Dong Kim
    • Journal of the Korean Society of Propulsion Engineers
    • /
    • v.7 no.4
    • /
    • pp.33-38
    • /
    • 2003
  • The Physics of the Plume-induced shock and separation Particularly at a high Plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with and without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG $\kappa$-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The control methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated The computational results show the main effect of the porous extension on plume-afterbody interactions is to restrain the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect rut plume interference.

A Study on Dynamic Characteristics of Hydraulic Motor Brake System with Counter Balance Valve (카운터 밸런스 밸브를 내장한 유압 모터 브레이크 시스템의 동특성)

  • Yun, So-Nam;Lee, Ill-Yeong
    • Journal of the Korean Society of Fisheries and Ocean Technology
    • /
    • v.29 no.3
    • /
    • pp.214-219
    • /
    • 1993
  • Counter balance valve is used as one part of hydraulic motor brake system. The function of this valve is to protect over-run or free falling of inertia load. But occasionally the brake system with counter balance valve makes some undesirable problems such as pressure surges or vibrations. These problems may hurt system safety and driver's conformability. Nevertheless, studies on dynamic characteristics of hydraulic system including counter balance valve are very rare, so further accumulation of research results are required. In this study, for the purpose of easy estimation about dynamic characteristics of hydraulic system including counter balance valve, precise formulation describing fluid dynamics and valve dynamics under various boundary conditions were made. The equations obtained in the preceding process include some parameters that must be got experimentally. Flow coefficients of valve and choke are the most significant ones among the parameters. So these parameters are obtained experimentally in this study, and experimental equations obtained from the experimental data were used for numerical calculation. The equations were analysed by numerical integration using Runge-Kutta method, because the equations contain various nonlinear terms. From the numerical analysis, it was verified that the dynamic response of counter balance valve and pressure variation at each elements can be estimated very easily. So the analysing method developed in this study enabled very easy estimating the relation between the performances of counter balance valve and various physical parameters related to the valve. Conclusively, it is said that the results obtained in this study can be used very usefully to develop a new type counter balance valve or to apply the valve to actual hydraulic system for various industrial equipments.

  • PDF

Development of 2D Tight-fitting Collar Pattern from 3D Scan Data of Various Types of Men's Dressform (남성 체형별 인대의 3차원 형상 데이터와 칼라 패턴 개발)

  • Jeong Yeon-Hee;Kim So-Young;Hong Kyung-Hi
    • Journal of the Korean Society of Clothing and Textiles
    • /
    • v.30 no.5 s.153
    • /
    • pp.722-732
    • /
    • 2006
  • The pattern making of the tight-fitting collars which often used in diving suits, dance wear, or cycle wear has not been fully established. To develop tight-fitting collar pattern directly from 3D images from the representative somatotypes, dressforms developed by Jaeun Jung were used. The 3D scan data of the four male dressforms were obtained using Exyma-1200. Triangle Simplification and the Runge-Kutta method were applied to reduce the 3D scan data points and to make the segmented triangular patches in a plane from 3D data. As results, apparent differences between the tight-fitting collar patterns obtained from the 3D scan data and the ordinary 2D collar patterns were found around the center back line. The curvatures of the center back line were higher in all types of the tight-fitting collar than in the ordinary collar pattern. Relative differences in the shape of collar lines among four representative Korean men were reported. To fit the curved shape of the back neckline, 1.8 cm should be reduced from the upper neckline in average. We suggested the direct pattern making method for the 2D tight-fitting collar patterns considering the 3D shape of various types of men's dressform.

Three-dimensional dynamics of vortex-induced vibration of a pipe with internal flow in the subcritical and supercritical regimes

  • Duan, Jinlong;Chen, Ke;You, Yunxiang;Wang, Renfeng;Li, Jinlong
    • International Journal of Naval Architecture and Ocean Engineering
    • /
    • v.10 no.6
    • /
    • pp.692-710
    • /
    • 2018
  • The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.

Modeling and Application of Chlorine Bulk Decay in Drinking Water Distribution System (배급수계통에서 잔류염소 감소 특성 및 적용연구)

  • Ahn, Jae-Chan;Park, Chang-Min;Koo, Ja-Yong
    • Journal of Korean Society of Water and Wastewater
    • /
    • v.19 no.4
    • /
    • pp.487-496
    • /
    • 2005
  • Chlorine bulk decay tests were carried out by bottle test under controlled conditions in a laboratory. Experiments were performed at different temperatures: $5^{\circ}C$, $15^{\circ}C$, $25^{\circ}C$, and the water temperatures when samples were taken from the effluent just before entering to its distribution system. 38 bulk tests were performed for water of Al (water treatment plant), 4 bulk tests for A2 (large service reservoir), and A3(pumping station). Residual chlorine concentrations in the amber bottles were measured over time till about 100 hours and bulk decay coefficients were evaluated by assuming first-order, parallel first-order, second-order. and $n^{th}-order$ reaction. The $n^{th}-order$ coefficients were obtained using Fourth-order Runge-Kutta Method. A good-fit by the average coefficient of determination ($R^2$) was first-order ($R^2=0.90$) < parallel first-order ($R^2{_{fast}}=0.92$, $R^2{_{slow}}=0.95$) < second-order ($R^2=0.95$) < $n^{th}-order$ ($R^2=0.99$). But if fast reaction of parallel first-order bulk decay were applied to the effluent of large service reservoir with ca. 20 hours of travel time and slow reaction in the water distribution system following the first 20 hours, parallel first-order bulk decay would be best and easy for application of water quality modeling technique.

Simulation Study on the Effects of Heating Rate and Particle Size Distribution for the Formation of the Agglomerate During CWM Combustion (CWM 연료의 연소시 입자 가열속도와 입자 크기가 CWM 응집물 형성에 미치는 영향에 관한 전산모사 연구)

  • Kim, Soo-Ho;Kim, Young-Hwan;Hwang, Kap-Sung;Hong, Song-Sun
    • Applied Chemistry for Engineering
    • /
    • v.9 no.3
    • /
    • pp.364-371
    • /
    • 1998
  • A theoretical model of particle agglomeration was developed to investigate the adhesive force between contiguous coal particles in CWM agglomerate. While heating bituminous coal to about $400^{\circ}C$ or above, the transient occurrence of plastic behavior of coal particles can be observed. The adhesive force in the process of agglomeration of coal particles was found to be proportional to the duration of plasticity of the particles. In the research, how the heating rate and the particle size distribution of CWM fuel influenced the formation of the agglomerate in CWM fuel at the heat-up stage was investigated by the model of particle agglomeration. Simulation program used to this experiment was RKG method and was programmed by Fortran. It was represented that by the model of particle agglomeration, the adhesive force in the process of the particle agglomeration in CWM fuel was inversely proportional to the heating rate but proportional to particle size.

  • PDF

Basic Study on the Heat Transfer During Rapid Freezing of Syobean Seed by Liquid Nitrogen

  • Kawano, Toshio;Nakano, Kohei;Murata, Satoshi
    • Proceedings of the Korean Society for Agricultural Machinery Conference
    • /
    • 1993.10a
    • /
    • pp.443-452
    • /
    • 1993
  • Direct freezing tests of soybean seed by liquid nitrogen were carried out at various moisture contents and the following important conclusions were drawn from the results of temperature measurements of soybean seed and photographs of bubbles generated on its surface : 1) Assuming that the temperature gradient in a soybean seed is negligible because of its small seed size and the freezing ratio is followed the Heiss's formula, and a differential equation based on the heat energy balance was introduced . The equation was easily solved by the Runge-Kutta-Gill method and the predicted values of the temperature were in good agreement with the observed data. 2) The photographs of bubble generation during freezing showed the boiling mode was nucleate, and then the most suitable formula on the nucleate boiling heat transfer was introduced from many formulate proposed up to now by fitting the calculated values based on the formula to the observed data. The formula used for the predict on of the seed temperature was as follows: $\frac{{\partial}T_s}{\partial\theta}\;=\;-\frac{{\alpha}(T_s\;-\;T_L)^{3.3}}{W(C_s\;-\;\frac{{\delta}m(CT_s\;+\;{\sigma})}{T_s^2})}$ where C = difference of the specific heat between pure ice and water m=moisture content of soybean seed $T_s$ = seed temperature $T_L$ = Temperature of liquid nitrogen W = mass of soybean seed $\alpha$ = proportional constant $\delta$ = constant depends on variety or the type of seed $\theta$ = time $\sigma$ = latent heat of melting of pure ice This study will give important information in the hydro-freezing technique by liquid nitrogen, available as a new technique of processing agricultural products in the near future.

  • PDF

2D Pattern Development of Body Surface from 3D Human Scan Data Using Standing and Cycling Postures (3D 스캔을 이용한 사이클 동작 전후 체표 변화 고찰 및 2D 전개 패턴의 비교)

  • Jeong, Yeonhee;Lee, Yejin
    • Korean Journal of Human Ecology
    • /
    • v.21 no.5
    • /
    • pp.975-988
    • /
    • 2012
  • Although the pattern development for tight-fitting clothing has been carried out using 3D data on humans, the pattern development using 3D scan data obtained for various postures still remains an interesting subject. In this study, we have developed the 2D pattern using the 3D human body reflecting standing and cycling postures. The 3D scan data of a subject was obtained using Cyberware. 2C-AN program(Triangle simplification and the Runge-Kutta method) was used in the system to reduce the 3D scan data points and to make segmented triangular patches in a plane from 3D data. As results, surface distance and area of each body part of standing and cycling postures were also provided for the future application of the functional clothing construction. The area of center piece on the front (c.front) decreased by $106.45cm^2$(-13.08%) and that of lateral piece(s.back) on the back increased by $144.96cm^2$(18.69%) in the patterns of cycling posture. The girth of neck and waist for the cycling posture increased by 0.88cm (3.92%) and 1.56cm(4.40%) respectively, and the that of thigh decreased by 1.01cm(-2.24%). The differences between the area in the 2D pattern obtained from the 3D scan data and that in the 3D scan surface data for standing and cycling postures were very small($-10.34cm^2$(-0.32%) and $-44.33cm^2$(-1.32%)).

A Study on the Influence of Nonlinearity Coefficients in Air-Bearing Spindle Parametric Vibration

  • Chernopyatov, Y.A.;Lee, C.M.;Chung, W.J.;Dolotov, K.S.
    • International Journal of Precision Engineering and Manufacturing
    • /
    • v.6 no.1
    • /
    • pp.51-58
    • /
    • 2005
  • The development of the high-efficiency machine-tools equipment and new cutting tool materials with high hardness, heat- and wear-resistance has opened the way to application of high-speed cutting process. The basic argument of using of high-speed cutting processes is the reduction of time and the respective increase of machining productivity. In this sense, the spindle units may be regarded as one of the most important units, directly affecting many parameters of high-speed machining efficiency. One of the possible types of spindle units for high-speed cutting is the air-bearing type. In this paper, we propose the mathematical model of the dynamic behavior of the air-bearing spindle. To provide the high-level of speed capacity and spindle rotation accuracy we need the adequate model of "spindle-bearings" system. This model should consider characteristics of the interactions between system components and environment. To find the working characteristics of spindle unit we should derive the equations of spindle axis movement under the affecting factors, and solve these equations together with equations which describe the behavior of lubricant layer in bearing (bearing stiffness equations). In this paper, the three influence coefficients are introduced, which describe the center of spindle mass displacement, angle of shaft rotation around the axes under the unit force application and that under the unit torque application. These coefficients are operated in the system of differential equations, which describes the spindle axis spatial movement. This system is solved by Runge-Kutta method. Obtained trajectories and amplitude-frequency characteristics were then compared to experimental ones. The analysis shows good agreement between theoretical and experimental results, which confirms that the proposed model of air-bearing spindle is correctis correct