• Journal of KSNVE
• /
• v.9 no.3
• /
• pp.525-534
• /
• 1999

The present paper proposes a generalized modal analysis procedure for non-uniform, distributed-parameter rotor-bearing systems. An exact element matrix is derived for a Timoshenko shaft model which contains rotary inertia, shear deformation, gyroscopic effect and internal damping. Complex coordinates system is adopted for the convenience in formulation. A generalized orthogonality condition is provided to make the modal decomposition possible. The generalized modal analysis by using a modal decomposition delivers exact and closed form solutions both for frequency and time responses. Two numerical examples are presented for illustrating the proposed method. The numerical study proves that the proposed method is very efficient and useful for the analysis of distributed-parameter rotor-bearing systems.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.13 no.10
• /
• pp.813-820
• /
• 2003

Conventional modal analysis techniques are known to be inappropriate for asymmetrical rotor systems, when the equations of motion are written in the stationary coordinates, due to the presence of time varying parameters. This paper presents a generalized modal analysis method for asymmetrical rotor systems in the stationary coordinates, employing the modulated coordinates and the lambda matrix formulation. A numerical example with a flexible asymmetric rotor model is provided to demonstrate the effectiveness of the proposed modal analysis method. As an application of the proposed method, modal analysis is also performed with an open cracked rotor system.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2003.05a
• /
• pp.526-531
• /
• 2003

Conventional modal analysis techniques are known to be inappropriate for asymmetrical rotor systems. when the equations of motion are written in the stationary coordinates, due to the presence of time varying parameters. This paper presents a generalized modal analysis method for asymmetrical rotor systems in the stationary coordinates, employing the modulated coordinates and the lambda matrix formulation. A numerical example with a flexible asymmetric rotor model is provided to demonstrate the effectiveness of the proposed modal analysis method. As an application of the proposed method, modal analysis is also performed with an open cracked rotor system.

• Journal of KSNVE
• /
• v.9 no.1
• /
• pp.103-112
• /
• 1999

A new modeling and analysis technique for one-dimensional distributed parameter systems is presented. First. discretized equations of motion in Laplace domain are derived by applying discretization methods for partial differential equations of a one-dimensional structure with respect to spatial coordinate. Secondly. the z and inverse z transformations are applied to the discretized equations of motion for obtaining a dynamic matrix for a uniform element. Four different discretization methods are tested with an example. Finally, taking infinite on the number of step for a uniform element leads to an exact dynamic matrix for the uniform element. A generalized modal analysis procedure for eigenvalue analysis and modal expansion is also presented. The resulting element dynamic matrix is tested with a numerical example. Another application example is provided to demonstrate the applicability of the proposed method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.195-201
• /
• 1996

This paper is mainly forcused on the application of modal analysis In analyze the geometrically non-linear buckling behaviors of large span spatial structures, and the evaluation of each eigen mode affected post-buckling behaviors and buckling loads. Modal analysis is applied . to derivation of the system matrices transforming actual displacement space into generalized coordinates space represented by coefficients multiplied in the linear combination of eigen modes which are independent and orthogonal each other. By using modal analysis method, it will be expected to save the calculating time by computer extremely. For example, we can obtain the satisfactorily good results by using about 7% of total eigen modes only in case of single layer latticed dome. And we can decrease the possibility of divergence on the bifurcation point in the calculation of post-buckling path. Arc-length method and Newton-Raphson iteration method are used to calculate the nonlinear equilibrium path.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.271-276
• /
• 1995

This paper presents a systematic method for the dynamic analysis of flexible mechanical systems containing closed kinematic loops. Kinematics between pairs of contiguous flexible bodies is described with the joint coordinates and the deformation modal coordinates. The cut-joint constraint equations associated with the closed kinematic loops are derived, simply using the geometric conditions. The equations of motions are initially written in terms of the joint and modal coordinates using the velocity transformation technique. Lagrange multipliers associated with the cut-joint constraints for closed-loop systems are then eliminated systematically using the generalized coordinate partitioning method, resulting to a minimal set of equations of motion.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.6-15
• /
• 2008

This paper deals with the aeroelastic instability of vibrating multiple blade rows under aerodynamic coupling with each other. A model composed of three blade rows, e.g., rotor-stator-rotor, where blades of the two rotor cascades are simultaneously vibrating, is considered. The displacement of a blade vibrating under aerodynamic force is expanded in a modal series with the natural mode shape functions, and the modal amplitudes are treated as the generalized coordinates. The generalized mass matrix and the generalized stiffness matrix are formulated on the basis of the finite element concept. The generalized aerodynamic force on a vibrating blade consists of the component induced by the motion of the blade itself and those induced not only by vibrations of other blades of the same cascade but also vibrations of blades in another cascade. To evaluate the aerodynamic forces, the unsteady lifting surface theory for the model of three blade rows is applied. The so-called k method is applied to determine the critical flutter conditions. A numerical study has been conducted. The flutter boundaries are compared with those for a single blade row. It is shown that the effect of the aerodynamic blade row coupling substantially modifies the critical flutter conditions.

• Structural Engineering and Mechanics
• /
• v.73 no.3
• /
• pp.271-286
• /
• 2020

Passive control technologies are commonly used in several areas to suppress structural vibrations by the addition of supplementary damping, and some modal damping may be heavy beyond critical damping even for regular structures with energy dissipation devices. The design of passive control structures is typically based on (complex) mode superposition approaches. However, the conventional mode superposition approach is predominantly applied to cases of under-critical damping. Moreover, when any modal damping ratio is equal or close to 1.0, the system becomes defective, i.e., a complete set of eigenvectors cannot be obtained such that some well-known algorithms for the quadratic eigenvalue problem are invalid. In this paper, a generalized complex mode superposition method that is suitable for under-critical, critical and over-critical damping is proposed and expressed in a unified form for structural displacement, velocity and acceleration responses. In the new method, the conventional algorithm for the eigenvalue problem is still valid, even though the system becomes defective due to critical modal damping. Based on the modal truncation error analysis, modal corrected methods for displacement and acceleration responses are developed to approximately consider the contribution of the truncated higher modes. Finally, the implementation of the proposed methods is presented through two numerical examples, and the effectiveness is investigated. The results also show that over-critically damped modes have a significant impact on structural responses. This study is a development of the original complex mode superposition method and can be applied well to dynamic analyses of non-classically damped systems.

• Current Optics and Photonics
• /
• v.5 no.5
• /
• pp.491-499
• /
• 2021

A scalar Fourier modal method for the numerical analysis of the scalar wave equation in inhomogeneous space with an arbitrary permittivity profile, is proposed as a novel theoretical embodiment of Fourier optics. The modeling of devices and systems using conventional Fourier optics is based on the thin-element approximation, but this approach becomes less accurate with high numerical aperture or thick optical elements. The proposed scalar Fourier modal method describes the wave optical characteristics of optical structures in terms of the generalized transmittance function, which can readily overcome a current limitation of Fourier optics.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.789-817
• /
• 2004

This paper deals with the modal analysis of rotational shell structures by means of the numerical solution technique known as the Generalized Differential Quadrature (G. D. Q.) method. The treatment is conducted within the Reissner first order shear deformation theory (F. S. D. T.) for linearly elastic isotropic shells. Starting from a non-linear formulation, the compatibility equations via Principle of Virtual Works are obtained, for the general shell structure, given the internal equilibrium equations in terms of stress resultants and couples. These equations are subsequently linearized and specialized for the rotational geometry, expanding all problem variables in a partial Fourier series, with respect to the longitudinal coordinate. The procedure leads to the fundamental system of dynamic equilibrium equations in terms of the reference surface kinematic harmonic components. Finally, a one-dimensional problem, by means of a set of five ordinary differential equations, in which the only spatial coordinate appearing is the one along meridians, is obtained. This can be conveniently solved using an appropriate G. D. Q. method in meridional direction, yielding accurate results with an extremely low computational cost and not using the so-called "delta-point" technique.