• Structural Engineering and Mechanics
• /
• v.37 no.5
• /
• pp.561-573
• /
• 2011

For efficient calculation of natural modes of structures, a numerical scheme which accelerates convergence of the subspace iteration method by employing accelerated starting Lanczos vectors was proposed in 2005. This paper is an extension of the study. The previous study simply showed feasibility of the proposed method by analyzing structures with smaller degrees of freedom. While, the present study verifies efficiency of the proposed method more rigorously by comparing closeness of conventional and accelerated starting vectors to genuine eigenvectors. This study also analyzes an example structure with larger degrees of freedom and more complex constraints in order to investigate applicability of the proposed method.

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

It is shown that a new Krylov subspace method for solving symmetric indefinite systems of linear equations can be obtained. We call the method as the projection method in this paper. The residual vector of the projection method is maintained at each iteration, which may be useful in some applications.

• Interaction and multiscale mechanics
• /
• v.1 no.2
• /
• pp.211-230
• /
• 2008

Owing to the growing size of the eigenvalue problem and the growing number of eigenvalues desired, solution methods of iterative nature are becoming more popular than ever, which however suffer from low efficiency and lack of proper convergence criteria. In this paper, three efficient iterative eigenvalue algorithms are considered, i.e., subspace iteration method, iterative Ritz vector method and iterative Lanczos method based on the cell sparse fast solver and loop-unrolling. They are examined under the mode error criterion, i.e., the ratio of the out-of-balance nodal forces and the maximum elastic nodal point forces. Averagely speaking, the iterative Ritz vector method is the most efficient one among the three. Based on the mode error convergence criteria, the eigenvalue solvers are shown to be more stable than those based on eigenvalues only. Compared with ANSYS's subspace iteration and block Lanczos approaches, the subspace iteration presented here appears to be more efficient, while the Lanczos approach has roughly equal efficiency. The methods proposed are robust and efficient. Large size tests show that the improvement in terms of CPU time and storage is tremendous. Also reported is an aggressive shifting technique for the subspace iteration method, based on the mode error convergence criteria. A backward technique is introduced when the shift is not located in the right region. The efficiency of such a technique was demonstrated in the numerical tests.

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

The Helmholtz equation is very important in physics and engineering. However, solution of the Helmholtz equation is in general known as a very difficult phenomenon. For if the ${\omega}$ is negative, the FDM discretized linear system becomes indefinite, whose solution by iterative method requires a very clever preconditioner. In this paper we assume that ${\omega}$ is nonnegative, and determine the optimal ${\rho}$ parameter for the three dimensional ADI iteration for the Helmholtz equation. The ADI(Alternating Direction Implicit) method is also getting new attentions due to the fact that it is very suitable to the vector/parallel computers, for example, as a preconditioner to the Krylov subspace methods. However, classical ADI was developed for two dimensions, and for three dimensions it is known that its convergence behaviour is quite different from that in two dimensions. So far, in three dimensions the so-called Douglas-Rachford form of ADI was developed. It is known to converge for a relatively wide range of ${\rho}$ values but its convergence is very slow. In this paper we determine the necessary conditions of the ${\rho}$ parameter for the convergence and optimal ${\rho}$ for the three dimensional ADI iteration of the Peaceman-Rachford form for the real Helmholtz equation with nonnegative ${\omega}$. Also, we conducted some experiments which is in close agreement with our theory. This straightforward extension of Peaceman-rachford ADI into three dimensions will be useful as an iterative solver itself or as a preconditioner to the the Krylov subspace methods, such as CG(Conjugate Gradient) method or GMRES(m).

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.677-680
• /
• 2004

This study investigates the buckling evaluation of connecting rods used in the diesel engine through finite element analysis. The Rankine formula, which is modified from classical Euler‘s formula, has been widely accepted in automotive industry to evaluate the buckling of connecting rods. Apparently, this formula is most suitable for the straight and idealized rod shape, and over-simplifies the geometric complexity associated with connecting rods. The subspace iteration method in FEA is used to predict the critical buckling stress of a connecting rod with certain slenderness ratio. To create models with various slenderness ratios for shank portion in the rod, the automatic meshing preprocessor was implemented. Results from FEA were verified by the experiments, in which the embedded strain gages measured for the connecting rod running at 4000rpm. The result indicates that the buckling prediction curve through FEA and experiment is effectively different from the curve of classical Rankine formula.

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

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

• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.5
• /
• pp.210-218
• /
• 2003

MDO (multidisciplinary design optimization) technology has been proposed and applied to solve large and complex optimization problems where multiple disciplinaries are involved. In this research. an MDO problem is defined for automobile design which has crashworthiness analyses. Crash model which are consisted of airbag, belt integrated seat (BIS), energy absorbing steering system .and safety belt is selected as a practical example for MDO application to vehicle system. Through disciplinary analysis, vehicle system is decomposed into structure subspace and occupant subspace, and coupling variables are identified. Before subspace optimization, values of coupling variables at given design point must be determined with system analysis. The system analysis in MDO is very important in that the coupling between disciplines can be temporary disconnected through the system analysis. As a result of system analysis, subspace optimizations are independently conducted. However, in vehicle crash, system analysis methods such as Newton method and fixed-point iteration can not be applied to one. Therefore, new system analysis algorithm is developed to apply to crashworthiness. It is conducted for system analysis to determine values of coupling variables. MDO algorithm which is applied to vehicle crash is MDOIS (Multidisciplinary Design Optimization Based on Independent Subspaces). Then, structure and occupant subspaces are independently optimized by using MDOIS.

• Nuclear Engineering and Technology
• /
• v.54 no.11
• /
• pp.4280-4295
• /
• 2022

In this paper, a new multi-group neutron-gamma transport calculation code system STRAUM-MATXST for complicated geometrical problems is introduced and its development status including numerical tests is presented. In this code system, the MATXST (MATXS-based Cross Section Processor for S_N Transport) code generates multi-group neutron and gamma cross sections by processing MATXS format libraries generated using NJOY and the STRAUM (S_N Transport for Radiation Analysis with Unstructured Meshes) code performs multi-group neutron-gamma coupled transport calculation using tetrahedral meshes. In particular, this work presents the recent implementation and its test results of the Krylov subspace methods (i.e., Bi-CGSTAB and GMRES(m)) with preconditioners using DSA (Diffusion Synthetic Acceleration) and TSA (Transport Synthetic Acceleration). In addition, the Krylov subspace methods for accelerating the energy-group coupling iteration through thermal up-scatterings are implemented with new multi-group block DSA and TSA preconditioners in STRAUM.

• ETRI Journal
• /
• v.46 no.3
• /
• pp.421-431
• /
• 2024

In this study, a direct position tracking method for non-circular (NC) signals using distributed passive arrays is proposed. First, we calculate the initial positions of sources using a direct position determination (DPD) approach; next, we transform the tracking into a compensation problem. The offsets of the adjacent time positions are calculated using a first-order Taylor expansion. The fusion calculation of the noise subspace is performed according to the NC characteristics. Because the proposed method uses the signal information from the previous iteration, it can realize automatic data associations. Compared with traditional DPD and two-step localization methods, our novel process has lower computational complexity and provides higher accuracy. Moreover, its performance is better than that of the traditional tracking methods. Numerous simulation results support the superiority of our proposed method.