• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.743-749
• /
• 2007

Algebraic substructuring (AS) is a state-of-the-art method in eigenvalue computations, especially for large size problems, but, originally, it was designed to calculate only the smallest eigenvalues. In this paper, an updated version of AS is proposed to calculate the interior eigenvalues over a specified range by using a shift value, which is referred to as the shifted AS. Numerical experiments demonstrate that the proposed method has better efficiency to compute numerous interior eigenvalues for the finite element models of structural problems than a Lanczos-type method.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.10a
• /
• pp.240-240
• /
• 2003

V-노치 균열에서 열하중이 작용하는 경우는 비제차형 경계조건의 문제가 되고, 이 조건에 대한 방정식의 일반해를 구하기 위해서 재차형 연립방정식에 대한 일반해(Homogeneous solution)와 비제차형 연립방정식에 대한 특수해(Particular solution)의 두 가지 해를 구할 수 있다. 이들 해는 V-노치 균열에 대한 고유치가 되고 이 고유치가 중복근을 가지게 되는 경우에는 로그항(1n[r])이 나타나게 되고 이 항에 의해서 응력을 무한대로 발산시키므로 이를 대수응력특이성이라 한다. 열하중이 작용할 때 대수응력특이성을 나타내는 로그항의 계수가 영(0)이 되어 대수응력특이성이 사라지게 되므로 V-노치 선단에서의 응력특이성은 고유치와 그에 대한 고유벡터에 의해 결정된다. 본 논문에서는 비정상상태 열하중이 가해지는 등방성 이종재료 내의 V-노치 균열문제에서 패기 각도와 이종재료의 기계적 성질에 의해 결정되는 응력특이성지수를 구하고 이에 대한 응력강도계수를 유한요소해석 프로그램인 ANSYS와 상반일 경로 적분법(RWCIM)을 이용하여 구하였다.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.6
• /
• pp.607-613
• /
• 2009

A new formulation of NDIF method to the algebraic eigenvalue problem is introduced to efficiently extract natural frequencies of arbitrarily shaped plates with the simply supported boundary condition. NDIF method, which was developed by the authors for the free vibration analysis of arbitrarily shaped membranes and plates, has the feature that it yields highly accurate natural frequencies compared with other analytical methods or numerical methods(FEM and BEM). However, NDIF method has the weak point that it needs the inefficient procedure of searching natural frequencies by plotting the values of the determinant of a system matrix in the frequency range of interest. A new formulation of NDIF method developed in the paper doesn't require the above inefficient procedure and natural frequencies can be efficiently obtained by solving the typical algebraic eigenvalue problem. Finally, the validity of the proposed method is shown in several case studies, which indicate that natural frequencies by the proposed method are very accurate compared to other exact, analytical, or numerical methods.

• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.4
• /
• pp.419-429
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies of thick, hyperboloidal shells of revolution. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components u/sub r/, u/sub θ/, u/sub z/ in the radial, circumferential, and axial directions, respectively, we taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the r and z directions. Potential(strain) and kinetic energies of the hyperboloidal shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Convergence to four digit exactitude is demonstrated for the first five frequencies of the hyperboloidal shells of revolution. Numerical results are tabulated for eighteen configurations of completely free hyperboloidal shells of revolution having two different shell thickness ratios, three variant axis ratios, and three types of shell height ratios. Poisson's ratio (ν) is fixed at 0.3. Comparisons we made among the frequencies for these hyperboloidal shells and ones which ate cylindrical or nearly cylindrical( small meridional curvature. ) The method is applicable to thin hyperboloidal shells, as well as thick and very thick ones.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.2
• /
• pp.119-129
• /
• 2004

A three dimensional (3D) method of analysis is presented for determining the free vibration frequencies and mode shapes of thick, circular and annular plates with nonlinear thickness variation along the radial direction. Unlike conventional plate theories, which are mathematically two dimensional (2D), the present method is based upon the 3D dynamic equations of elasticity. Displacement components u/sub s/, u/sub z/, and u/sub θ/ in the radial, thickness, and circumferential directions, respectively, are taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the s and z directions. Potential (strain) and kinetic energies of the plates are formulated, and the Ritz method is used to solve the eigenvalue problem thus yielding upper bound values of the frequencies by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Convergence to four digit exactitude is demonstrated for the first five frequencies of the plates. Numerical results we presented for completely free, annular and circular plates with uniform linear, and quadratic variations in thickness. Comparisons are also made between results obtained from the present 3D and previously published thin plate (2D) data.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.197-206
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies and mode shapes of solid and hollow hemispherical shells of revolution of arbitrary wall thickness having arbitrary constraints on their boundaries. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components μ/sub Φ/, μ/sub z/, and μ/sub θ/ in the meridional, normal, and circumferential directions, respectively, are taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the Φ and z directions. Potential (strain) and kinetic energies of the hemispherical shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies obtained by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Novel numerical results are presented for solid and hollow hemispheres with linear thickness variation. The effect on frequencies of a small axial conical hole is also discussed. Comparisons are made for the frequencies of completely free, thick hemispherical shells with uniform thickness from the present 3-D Ritz solutions and other 3-D finite element ones.