• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.17 no.3 s.120
• /
• pp.257-263
• /
• 2007

Detection of structural damage is an inverse problem in structural engineering. There are three main questions in the damage detection: existence, location and extent of the damage. In concept, the natural frequency and mode shapes of any structure must satisfy an eigenvalue problem. But, if a potential damage exists in a structure, an error resulting from the substitution of the refined analytical finite element model and measured modal data into the structural eigenvalue equation will occur, which is called the residual modal forces, and can be used as an indicator of potential damage in a structure. In this study, a useful damage detection method is proposed and compared with other two methods. Two degree-of-freedom system and Cantilever beam are used to demonstrate the approach. And the results of three introduced method are compared.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.205-216
• /
• 1998

A solution method is presented to solve the eigenvalue problem arising in the dynamic analysis of nonclassicary damped structural systems with multiple eigenvalues. The proposed method is obtained by applying the modified Newton-Raphson technique and the orthonormal condition of the eigenvectors to the linear eigenproblem through matrix augmentation of the quadratic eigenvalue problem. In the iteration methods such as the inverse iteration method and the subspace iteration method, singularity may be occurred during the factorizing process when the shift value is close to an eigenvalue of the system. However, even though the shift value is an eigenvalue of the system, the proposed method provides nonsingularity, and that is analytically proved. Since the modified Newton-Raphson technique is adopted to the proposed method, initial values are need. Because the Lanczos method effectively produces better initial values than other methods, the results of the Lanczos method are taken as the initial values of the proposed method. Two numerical examples are presented to demonstrate the effectiveness of the proposed method and the results are compared with those of the well-known subspace iteration method and the Lanczos method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.1
• /
• pp.13-20
• /
• 2009

The IPSAP which is a finite element analysis program has been developed for high parallel performance computing. This program consists of various analysis modules - stress, vibration and thermal analysis module, etc. The M orthogonal block Lanczos algorithm with shiftinvert transformation is used for solving eigenvalue problems in the vibration module. And the multifrontal algorithm which is one of the most efficient direct linear equation solvers is applied to factorization and triangular system solving phases in this block Lanczos iteration routine. In this study, the performance enhancement procedures of the IPSAP are composed of the following stages: 1) communication volume minimization of the factorization phase by modifying parallel matrix subroutines. 2) idling time minimization in triangular system solving phase by partial inverse of the frontal matrix and the LCM (least common multiple) concept.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.03a
• /
• pp.481-490
• /
• 2003

This paper presents an efficient eigensolution method for non-proportionally damped systems. The proposed method is obtained by applying the accelerated Newton-Raphson technique and the orthonormal condition of the eigenvectors to the linearized form of the quadratic eigenproblem. A step length and a selective scheme are introduced to increase the convergence of the solution. The step length can be evaluated by minimizing the norm of the residual vector using the least square method. While the singularity may occur during factorizing process in other iteration methods such as the inverse iteration method and the subspace iteration method if the shift value is close to an exact eigenvalue, the proposed method guarantees the nonsingularity by introducing the orthonormal condition of the eigenvectors, which can be proved analytically. A numerical example is presented to demonstrate the effectiveness of the proposed method.

• Journal of Korean Association for Spatial Structures
• /
• v.10 no.2
• /
• pp.87-94
• /
• 2010

Form-Finding of tensegrity structures by eigenvalue problem is presented, In ardor to maintain the structures stable, "Form-Finding" should be performed. The types of analytical methods are known to solve this phenomenon: One is to use force density method, and the other is to apply so called, generalized inverse method. In this paper, new form finding methods are presented to obtain the self-equilibrium stress of the tensegrity structures. This method is based on the equilibrium equation of the all of the joint and the governing equation is formulated as eigonvalue problem. In order to verify this approach, numerical example(tensegrity structures) are compared with others calculated by previous methods. The solution by present method is shown identical results. Furthermore, the developed process to find the results is more efficient than previous approaches.

• Journal of applied mathematics & informatics
• /
• v.11 no.1_2
• /
• pp.305-316
• /
• 2003

Recently, an iterative algorithm for finding the interior eigenvalues of a definite matrix by CG-type method has been proposed. This method compares to the inverse power method. The given matrices A, and B are assumed to be large and sparse, and SPD( Symmetric Positive Definite) The CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for smallest eigenvalue. Also, it is very amenable to parallel computations, like the CG method for the linear systems. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. But for parallel computations we need to find an efficient parallel preconditioner. Our candidates we ILU(0) in the wave-front order, ILU(0) in the multi-coloring order, Point-SSOR(Symmetric Successive Overrelaxation), and Multi-Color Block SSOR preconditioner. Wavefront order is a simple way to increase parallelism in the natural order, and Multi-coloring realizes a parallelism of order(N), where N is the order of the matrix. Another choice is the Multi-Color Block SSOR(Symmetric Successive OverRelaxation) preconditioning. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test problem was drawn from the discretizations of partial differential equations by finite difference methods. The results show that for small number of processors Multi-Color ILU(0) has the best performance, while for large number of processors Multi-Color Block SSOR performs the best.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.371-383
• /
• 1999

An efficient and numerically stable eigensolution method for structures with multiple natural frequencies is presented. The proposed method is developed by improving the well-known subspace iteration method with shift. A major difficulty of the subspace iteration method with shift is that because of singularity problem, a shift close to an eigenvalue can not be used, resulting in slower convergence. In this paper, the above singularity problem has been solved by introducing side conditions without sacrifice of convergence. The proposed method is always nonsingular even if a shift is on a distinct eigenvalue or multiple ones. This is one of the significant characteristics of the proposed method. The nonsingularity is proved analytically. The convergence of the proposed method is at least equal to that of the subspace iteration method with shift, and the operation counts of above two methods are almost the same when a large number of eigenpairs are required. To show the effectiveness of the proposed method, two numerical examples are considered.

• Nuclear Engineering and Technology
• /
• v.3 no.3
• /
• pp.121-128
• /
• 1971

A unified modal-nodal expansion of tile angular distribution of neutron flux in one spatial dimension is considered, following the proposal of Harms. Several standard nodal and/or modal methods of analysis are shown to be specializations of this technique. The modal-nodal moment from of the mono-energetic transport equation with isotropic sources and scattering is derived and the infinite medium eigenvalue problem solved. The technique is shown to yield results which approximate the exact value of the inverse diffusion length in non-multiplying media more accurately than standard methods of equal or somewhat greater computational complexity.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.3
• /
• pp.218-224
• /
• 2011

A modal method is presented for the investigation of the effects of measurement errors in damage detection for dynamic structural systems. The structural modifications to the baseline system result in the response changes of the perturbed structure, which are measured to determine a unique system in the inverse problem of damage detection. If the numerical modal data are exact, mathematical programming techniques can be applied to obtain the accurate structural changes. In practice, however, the associated measurement errors are unavoidable, to some extent, and cause significant deviations from the correct perturbed system because of the intrinsic instability of eigenvalue problem. Hence, a self-equilibrating inverse system is allowed to drift in the close neighborhood of the measured data. A numerical example shows that iterative procedures can be used to search for the damaged structural elements. A small set of selected degrees of freedom is employed for practical applicability and computational efficiency.

• Journal of Korean Society of Steel Construction
• /
• v.18 no.5
• /
• pp.575-584
• /
• 2006

In this paper, a shear-flexible finite element model is developed for the buckling analysis of axially loaded, thin-walled composite I-beams. Based on an orthogonal Cartesian coordinate system, the displacement fields are defined using the first-order shear-deformable beam theory. The derived element takes into account flexural shear deformation and torsional warping deformation. Three different types of beam elements, namely, the two-noded, three-noded, and four-noded beam elements, were developed to solve the governing equations. An inverse iteration with shift eigenvalue solution was used to solve the resulting linearized buckling problem. A parametric study was conducted to show the importance of shear flexibility and fiber orientation on the buckling behavior of thin-walled composite beams. A good agreement was obtained among the proposed shear-flexible model, other results available in literature, and the finite element solution.