• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.103-112
• /
• 2000

Numerical simulations of unsteady opposed-flow flames are performed using an adaptive time integration method designed for differential-algebraic systems. The compressibility effect is considered in deriving the system of equations, such that the numerical difficulties associated with a high-index system are alleviated. The numerical method is implemented for systems with detailed chemical mechanisms and transport properties by utilizing the Chemkin software. Two test simulations are performeds hydrogen/air diffusion flames with an oscillatory strain rate and transient ignition of methane against heated air. Both results show that the rapid transient behavior is successfully captured by the numerical method.

• Journal of the Korean Society of Industry Convergence
• /
• v.12 no.3
• /
• pp.149-155
• /
• 2009

A new finite element equation is derived by applying quadratic and cubic time integration scheme to the variational formulation in time-integral for the analysis of the transient elastodynamic problems to increase the numerical accuracy and stability. Emphasis is focused on methodology for cubic time integration scheme procedure which are never presented before. In this semidiscrete approximations of the field variables, the time axis is divided equally and quadratic and cubic time variation is assumed in those intervals, and space is approximated by the usual finite element discretization technique. It is found that unconditionally stable numerical results are obtained in case of the cubic time variation. Some numerical examples are given to show the versatility of the presented formulation.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.29 no.10
• /
• pp.89-96
• /
• 2015

Explicit numerical integration methods for power system transient stability simulation require very small time steps to avoid numerical instability. The EXST1 exciter model is a primary source of fast dynamics in power system transients. In case of the EXST1, the required small integration time step for entire system simulation increases the computational demands in terms of running time and storage. This paper presents a practical exciter model reduction approach which allows the increase of the required step size and thus the method can decrease the computational demands. The fast dynamics in the original EXST1 are eliminated in the reduced exciter model. The use of a larger time step improves the computational efficiency. This paper describes the way to eliminate the fast dynamics from the original exciter model based on linear system theory. In order to validate the performance of the proposed method, case studies with the GSO-37 bus system are provided. Comparisons between the original and reduced models are made in simulation accuracy and critical clearing time.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.9
• /
• pp.674-684
• /
• 2002

In order to decrease remarkably the computation time and storage used in the direct integration method without the loss of accuracy, authors suggest a new transient analysis algorithm. This algorithm is derived from the combination of three techniques, that is, the transfer technique of the transfer stiffness coefficient method, the modeling technique of the finite element method, and the numerical integration technique of the Newmark method. In this paper, the transient analysis algorithm of a frame structure is formulated by the proposed method. The accuracy and computation efficiency of the proposed method are demonstrated through the comparing with the computation results by the direct integration method for three computation models under various excitations.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1993.04a
• /
• pp.52-56
• /
• 1993

An efficient method is presented for determining transient responses of locally nonlinear structures using substructure eigenproperties and Lagrange multiplier technique. The method is based upon the mode synthesis formulation procedure, but does not construct the equations of motion of the combined whole structure compared with the conventional methods. For modal bases of each linear substructure, either fixed or free interface modes can be employed. The transient analysis is based upon the recurrence discrete-time state equations and offers the simplicity of the Euler integration method without requiring small time increment and iterative solution procedure. Numerical examples reveal that the method is very accurate and efficient in calculating transient responses compared with the direct numerical integration method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2457-2466
• /
• 1993

A simple method is presented for determining transient responses of locally nonlinear structures using substructure eigenproperties and Lagrange multiplier technique. Although the method is based upon the mode synthesis formulation procedure, the equations of the combined whole structure are not constructed compared with the conventional methods. Lagrange multi-pliers are used to enforce the conditions of geometric compatibility between the substructure interfaces and they are treated as external forces on each substructure itself. Substructure eigenvalue problem is defined with the substructure interface free of fixed. The transient analysis is based upon the recurrence discrete-time state equations and offers the simplicity of the Euler integration method without requiring small time increment and iterative solution procedure. Numerical examples reveal that the method is very accurated and efficient in calculating transient responses compared with the direct numerical integration method.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.357-367
• /
• 2020

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.11a
• /
• pp.242-247
• /
• 2004

In this study was proposed a transient response analysis technique of a rotor system, applying the generalized FE modeling method of a rotor-bearing system considering a base-transferred shock force and together the state-space Newmark method of direct time integration scheme based on the average velocity concept. Experiments were performed to a test rig of a mock-up rotor-bearing system with series of half-sine shock waves imposed by an electromagnetic shaker, and quantitative error analyses between analytical and experimental results were carried out. The transient reponses of the rotor were sensitive to duration times and shape-qualities of the shock waves, and overally the analytical results agreed quite well with the experimental ones. Particularly, in cases that the frequencies, $1/(2{\times}duration\;time)$, of the shock waves were close to the critical speed of the rotor-bearing system, resonances occurred and the transient responses of the rotor were amplified.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.407-412
• /
• 1998

In this work, a new time integration method is proposed using the generalized derivative concept to simulate the dynamic phenomena having sudden constraint occurring in dynamic contact/impact problems. By the adoption of the generalized derivative concept and jump assumption, discontinuity can be incorporated in time integration and as a result, the algorithm does not need any other special consideration of jumps in dynamic field variables due to sudden constraint like dynamic contact-release conditions. To observe the characteristics of the proposed time integration method, the stability and convergence analyses are carried out. In numerical tests, several dynamic contact/impact problems are analyzed by straightforward application of the proposed time integration method with the exterior penalty method.

• Structural Engineering and Mechanics
• /
• v.38 no.4
• /
• pp.429-441
• /
• 2011

In large-scale problem, a huge size of computational resources is needed for a reliable solution which represents the detailed description of dynamic behavior. Recently, eigenvalue reduction schemes have been considered as important technique to resolve computational resource problems. In addition, the efforts to advance an efficiency of reduction scheme leads to the development of the multi-level system condensation (MLSC) which is initially based on the two-level condensation scheme (TLCS). This scheme was proposed for approximating the lower eigenmodes which represent the global behavior of the structures through the element-level energy estimation. The MLSC combines the multi-level sub-structuring scheme with the previous TLCS for enhancement of efficiency which is related to computer memory and computing time. The present study focuses on the implementation of the MLSC on the direct time response analysis and the frequency response analysis of structural dynamic problems. For the transient time response analysis, the MLSC is combined with the Newmark's time integration scheme. Numerical examples demonstrate the efficiency of the proposed method.