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

A simplified method for the computation of first second and higher order derivatives of eigenvalues and eigenvectors derivatives associated with repeated eigenvalues is presented. Adjacent eigenvectors and orthonormal conditions are used to compose an algebraic equation whose order is (n+m)x(n+m), where n is the number of coordinates and m is the number of multiplicity of the repeated eigenvalues. The algebraic equation developed can be used to compute derivatives of both eigenvalues and eigenvectors simultaneously. Since the coefficient matrix in the proposed algebraic equation is non-singular, symmetric and based on N-space it is numerically stable and very efficient compared to previous methods. This method can be consistently applied to structural systems with structural design parameters and mechanical systems with lumped design parameters. To verify the effectiveness of the proposed method, the finite element model of the cantilever beam is considered.

• Journal of Korean Society of Steel Construction
• /
• v.8 no.3 s.28
• /
• pp.21-34
• /
• 1996

A numerical method is presented for computation of eigenvector derivatives used an iterative procedure with guaranteed convergence. An approach for treating the singularity in calculating the eigenvector derivatives is presented, in which a shift in each eigenvalue is introduced to avoid the singularity. If the shift is selected properly, the proposed method can give very satisfactory results after only one iteration. A criterion for choosing an adequate shift, dependent on computer hardware is suggested ; it is directly dependent on the eigenvalue magnitudes and the number of bits per numeral of the computer. Another merit of this method is that eigenvector derivatives with repeated eigenvalues can be easily obtained if the new eigenvectors are calculated. These new eigenvectors lie "adjacent" to the m (number of repeated eigenvalues) distinct eigenvectors, which appear when the design parameter varies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The results are compared with those of Nelson's method which can find the exact eigenvector derivatives. For the case of repeated eigenvalues, a cantilever beam is considered. The results are compared with those of Dailey's method which also can find the exact eigenvector derivatives. The design parameter of the cantilever plate is its thickness, and that of the cantilever beam its height.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.117-124
• /
• 2001

A simplified method fur the eigenpair sensitivities of damped system with multiple eigenvalues is presented. This approach employs a reduced equation to determine the sensitivities of eigenpairs of the damped vibratory systems with multiple natural frequencies. In the proposed method, adjacent eigenvectors and orthonormal conditions are used to compute an algebraic equation whose order is (n+m)x(n+m), where n is the number of coordinates and m the number of multiplicity of multiple natural frequencies. The proposed method is an improved Lee and Jung's method which was developed previously. Two equations are used to find eigenvalue derivatives and eigenvector derivatives in Lee and Jung's method. A significant advantage of this approach over Lee and Jung's method is that one algebraic equation newly developed is enough to compute such eigenvalue derivatives and eigenvector derivatives. This method can be consistently applied to both structural systems with structural design parameters and mechanical systems with lumped design parameters. To demonstrate the theory of the proposed method and its possibilities in the case of multiple eigenvalues, the finite element model of the cantilever beam and 5-DOF mechanical system in the case of a non-proportionally damped system are considered as numerical examples. The design parameter of the cantilever beam is its height. and that of the 5-DOF mechanical system is a spring.

• Structural Engineering and Mechanics
• /
• v.4 no.3
• /
• pp.277-301
• /
• 1996

A finite element model of a beam element with flexible connections is used to investigate the effect of the randomness in the stiffness values on the modal properties of the structural system. The linear behavior of the connections is described by a set of random fixity factors. The element mass and stiffness matrices are function of these random parameters. The associated eigenvalue problem leads to eigenvalues and eigenvectors which are also random variables. A second order perturbation technique is used for the solution of this random eigenproblem. Closed form expressions for the 1st and 2nd order derivatives of the element matrices with respect to the fixity factors are presented. The mean and the variance of the eigenvalues and vibration modes are obtained in terms of these derivatives. Two numerical examples are presented and the results are validated with those obtained by a Monte-Carlo simulation. It is found that an almost linear statistical relation exists between the eigenproperties and the stiffness of the connections.

• Structural Engineering and Mechanics
• /
• v.36 no.1
• /
• pp.37-55
• /
• 2010

For large-scale structures, the calculation of the eigensolution and the eigensensitivity is usually very time-consuming. This paper develops the Kron's substructuring method to compute the first-order derivatives of the eigenvalues and eigenvectors with respect to the structural parameters. The global structure is divided into several substructures. The eigensensitivity of the substructures are calculated via the conventional manner, and then assembled into the eigensensitivity of the global structure by performing some constraints on the derivative matrices of the substructures. With the proposed substructuring method, the eigenvalue and eigenvector derivatives with respect to an elemental parameter are computed within the substructure solely which contains the element, while the derivative matrices of all other substructures with respect to the parameter are zero. Consequently this can reduce the computation cost significantly. The proposed substructuring method is applied to the GARTEUR AG-11 frame and a highway bridge, which is proved to be computationally efficient and accurate for calculation of the eigensensitivity. The influence of the master modes and the division formations are also discussed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.225-233
• /
• 1995

This paper presents an efficient numerical method for computation of eigenpair derivatives for the real symmetric eigenvalue problem with distinct and multiple eigenvalues. The method has very simple algorithm and gives an exact solution. Furthermore, it saves computer storage and CPU time. The algorithm preserves the symmetry and band of the matrices, allowing efficient computer storage and solution techniques. Thus, the algorithm of the proposed method will be inserted easily in the commercial FEM codes. Results of the proposed method for calculating the eigenpair derivatives are compared with those of Rudisill and Chu's method and Nelson's method which is efficient one in the case of distinct natural frequencies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The design parameter of the cantilever plate is its thickness. For the eigenvalue problem with multiple natural frequencies, the adjacent eigenvectors are used in the algebraic equation as side conditions, they lie adjacent to the m (multiplicity of multiple natural frequency) distinct eigenvalues, which appear when design parameter varies. As an example to demonstrate the efficiency of the proposed method in the case of multiple natural frequencies, a cantilever beam is considered. Results of the proposed method fDr calculating the eigenpair derivatives are compared with those of Bailey's method (an amendation of Ojalvo's work) which finds the exact eigenvector derivatives. The design parameter of the cantilever beam is its height. Data is persented showing the amount of CPU time used to compute the first ten eigenpair derivatives by each method. It is important to note that the numerical stability of the proposed method is proved.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.5
• /
• pp.707-718
• /
• 1997

This paper presents an efficient numerical method for the computation of eigenpair derivatives for a real symmetric eigenvalue problem with distinct and multiple eigenvalues. The method has a very simple algorithm and gives an exact solution. Furthermore, it saves computer sotrage and CPU time. The algorithm preserves not only the symmetricity but also the band width of the matrices, allowing efficient computer storage and solution techniques. Results from the proposed method for calculating the eigenpair derivatives are compared with those from Rudisill and Chu's method and Nelson's method which is known efficient one in the case of distinct natural frequencies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The design parameter of the cantilever plate is its thickness. For the eigenvalue problem with multiple natural frequencies, the adjacent eigenvectors are used in the algebraic equation as side conditions, lying adjacent to the multiplicity of multiple natural frequency distinct eigenvalues, which appear when design parameter varies. A cantilever beam is used to demonstrate the efficiency of the proposed method in the case of multiple natural frequencies. Results form the proposed method for calculating the eigenpair derivatives are compared with those from Dailey's method(an amendation of Ojalvo's work) which finds the exact eigenvector derivatives. The design parameter of the cantilever beam is its height. Data is presented showing the amount of CPU time used to compute the first ten eigenpair derivatives by each method. It is important to note that the numerical stability of the proposed method is proved.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.10a
• /
• pp.515-522
• /
• 2001

A simplified method is presented for the computation of eigenvalue and eigenvector derivatives associated with repeated eigenvalues. In the proposed method, adjacent eigenvectors and orthonormal conditions are used to compose an algebraic equation whose order is (n+m)x(n+m), where n is the number of coordinates and m is the number of multiplicity of the repeated eigenvalue. One algebraic equation developed can be computed eigenvalue and eigenvector derivatives simultaneously. Since the coefficient matrix of the proposed equation is symmetric and based on N-space, this method is very efficient compared to previous methods. Moreover the numerical stability of the method is guaranteed because the coefficient matrix of the proposed equation is non-singular, This method can be consistently applied to both structural systems with structural design parameters and mechanical systems with lumped design parameters. To verify the effectiveness of the proposed method, the finite element model of the cantilever beam and a 5-DOF mechanical system in the case of a non-proportionally damped system are considered as numerical examples. The design parameter of the cantilever beam is its width, and that of the 5-DOF mechanical system is a spring.

• Computational Structural Engineering
• /
• v.7 no.2
• /
• pp.131-138
• /
• 1994

Derivatives of eigenvalues and eigenvectors for statistical properties is presented. Dynamic response of random system including uncertainties for the design variable is calculated with the first order perturbation method to original governing equation. In optimal design methods, there is fundamental requirement for design gradients. A method for calculating the sensitivity coefficients is developed using the governing equation and first order perturbed equation.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.796-801
• /
• 2002

The goal of this research is the shape design of the valve using a computer simulation. For an analysis a basic mathematical model describing compression cycle is considered as consisting of five sets of coupled equations. These are the volume equation (kinematics), valve dynamic equation (dynamics), ideal gas equation (thermodynamics), Bernoulli equation (fluid dynamics), and dynamic equation of fluid particle based on Helmholtz equation (acoustics). Valve motion is made by the superposition of free vibration modes obtained by the finite element method. That is, the eigenvalues and eigenvectors are the sufficient modeling factors fur the valve in the simulation program. Thus, to design a shape of the valve, shape design sensitivity through chain-ruled derivatives is considered from two sensitivity coefficients, one is the design sensitivity of the capability of compressor with respect to the eigenvalues of the valve, and the other is the design sensitivity of the eigenvalue with respect to the shape change of the valve. In this research, the continuum design sensitivity analysis concepts are used for the latter.