• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.37-49
• /
• 1999

In this paper we develop a general Homotopy method called the Group Homotopy method to solve the symmetric eigenproblem. The Group Homotopy method overcomes notable drawbacks of the existing Homotopy method, namely, (i) the possibility of breakdown or having a slow rate of convergence in the presence of clustering of the eigenvalues and (ii) the absence of any definite criterion to choose a step size that guarantees the convergence of the method. On the other hand, We also have a good approximations of the largest eigenvalue of a Symmetric matrix from Lanczos algorithm. We apply it for the largest eigenproblem of a very large symmetric matrix with a good initial points.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.283-290
• /
• 1998

A solution method is presented to solve the eigenproblem arising in tile dynamic analysis of non-proportional damping systems with symmetric matrices. The method is based on tile use of Lanczos method to generate a Krylov subspace of trial vectors, witch is then used to reduce a large eigenvalue problem to a much smaller one. The method retains the η order quadratic eigenproblem, without the need to the method of matrix augmentation traditionally used to cast the problem as a linear eigenproblem of order 2n. In the process, the method preserves tile sparseness and symmetry of the system matrices and does not invoke complex arithmetics, therefore, making it very economical for use in solving large problems. Numerical results are presented to demonstrate the efficiency and accuracy of the method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.7 no.1
• /
• pp.25-31
• /
• 2003

The concepts of stability of algorithms for solving the symmetric and generalized symmetric-definite eigenproblems are discussed. An algorithm for solving the symmetric eigenproblem $Ax={\lambda}x$ is stable if the computed solution z is the exact solution of some slightly perturbed system $(A+E)z={\lambda}z$. We use both normwise approach and componentwise way of measuring the size of the perturbations in data. If E preserves symmetry we say that an algorithm is strongly stable (in a normwise or componentwise sense, respectively). The relations between the stability and strong stability are investigated for some classes of matrices.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.2
• /
• pp.99-110
• /
• 2000

This paper is devoted to the parallelism of a numerical matrix eigenvalue problem. The eigenproblem arises in a variety of applications, including engineering, statistics, and economics. Especially we try to approach the industrial techniques from mathematical modeling. This paper has developed a parallel algorithm to find all eigenvalues. It is contributed to solve a specific practical problem, a vibration problem in the industry. Also we compare the runtime between the serial algorithm and the parallel algorithm for the given problems.

• Computational Structural Engineering : An International Journal
• /
• v.1 no.1
• /
• pp.31-38
• /
• 2001

In the case of the non-proportionally damped system such as the soil-structure interaction system, the structural control system and composite structures, the eigenproblem with the damping matrix should be necessarily performed to obtain the exact dynamic response. However, most of the eigenvalue analysis methods such as the subspace iteration method and the Lanczos method may miss some eigenvalues in the required ones. Therefore, the eigenvalue analysis method must include a technique to check the missed eigenvalues to become the practical tools. In the case of the undamped or proportionally damped system the missed eigenvalues can easily be checked by using the well-known Sturm sequence property, while in the case of the non-proportionally damped system a checking technique has not been developed yet. In this paper, a technique of checking the missed eigenvalues for the eigenproblem with the damping matrix is proposed by applying the argument principle. To verify the effectiveness of the proposed method, two numerical examples are considered.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.11 s.242
• /
• pp.1527-1533
• /
• 2005

Three methods for design sensitivity analysis such as finite difference method(FDM), direct differentiation method(DDM) and adjoint variable method(AVM) are well known. FDM and DDM for design sensitivity analysis cost too much when the number of design variables is too large. An AVM is required to compute adjoint variables from the simultaneous linear system equation, the so-called adjoint equation. Because the adjoint equation is independent of the number of design variables, an AVM is efficient for when number of design variables is too large. In this study, AVM has been extended to the eigenproblem of damped systems whose eigenvlaues and eigenvectors are complex numbers. Moreover, this method is implemented into a commercial finite element analysis program by means of the semi-analytical method to show applicability of the developed method into practical structural problems. The proposed_method is compared with FDM and verified its accuracy for analytical and practical cases.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.64-73
• /
• 2000

We present finite element equations in the Laplace-domain for linear viscoelastic and viscoelstically damped structures governed by a constitutive equation involving factional order derivative opeartors. These equations yield a nonstandard eigenproblem consisted of frequency dependent stiffness matrix. To solve this nonstandard eigenproblem we suggest an eigenvalue iteration procedure in the Laplace-domain. Improved Zenor and GHM material function type constitutive equations in the Laplace-domain are also available for this procedure. From above equations, complex eigenvalues and complex eigenvectors are obtained. Using obtained eigenvalues and eigenvectors, time domain analysis is performed by means of mode superposition. Finally, finite element solutions of viscoelastic and viscoeleastically damped sandwich beam are presented as an example.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.467-474
• /
• 2006

We present and study the stability and convergence of a deflation-preconditioned conjugate gradient(PCG) scheme for the interior generalized eigenvalue problem $Ax = {\lambda}Bx$, where A and B are large sparse symmetric positive definite matrices. Numerical experiments are also presented to support our theoretical results.

• Journal of applied mathematics & informatics
• /
• v.9 no.1
• /
• pp.331-339
• /
• 2002

A preconditioned conjugate gradient(PCG) scheme with the aid of deflation for computing a few of the smallest eigenvalues arid their corresponding eigenvectors of the large generalized eigenproblems is considered. Topically there are two types of deflation techniques, the deflation with partial shifts and an arthogonal deflation. The efficient way of determining partial shifts is suggested and the deflation-PCG schemes with various partial shifts are investigated. Comparisons of theme schemes are made with orthogonal deflation-PCG, and their asymptotic behaviors with restart operation are also discussed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.133-140
• /
• 2002

This paper investigates the applicability of the modified Lanczos method using the power technique, which was developed in the field of quantum physics, to the eigenproblem in the field of engineering mechanics by introducing matrix-powered Lanczos recursion and numerically evaluating the suitable power value. The matrix-powered Lanczos method has better convergence and less operation count than the conventional Lanczos method. By analyzing four numerical examples, the effectiveness of the matrix-powered Lanczos method is verified and the appropriate matrix power is also recommended.