• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1998년도 가을 학술발표회 논문집
• /
• pp.431-438
• /
• 1998

An efficient solution method is presented to solve the eigenvalue problem arising in tile dynamic analysis of non-proportionally damped structural systems with multiple or close eigenvalues. The proposed method is obtained by applying the modified Newton-Raphson technique and the orthonormal condition of the eigenvectors to the quadratic eigenvalue problem. Even if the shift value is an eigenvalue of the system, the proposed method guarantees nonsingularity, which is analytically proved. The initial values of the proposed method can be taken as the intermediate results of iteration methods or results of approximate methods. Two numerical examples are also 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.

• Coupled systems mechanics
• /
• 제1권2호
• /
• pp.183-204
• /
• 2012

A method for the analysis of wave scattering in 3-D elastic full space is developed by means of the coupled boundary-volume integral equation, which takes into account the effects of both the boundary of inclusions and the uctuation of the wave field. The wavenumber domain formulation is used to construct the Krylov subspace by means of FFT. In order to achieve the wavenumber domain formulation, the boundary-volume integral equation is transformed into the volume integral equation. The formulation is also focused on this transform and its numerical implementation. Several numerical results clarify the accuracy and effectiveness of the present method for scattering analysis.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.111-116
• /
• 2007

This paper investigates and compares the natural modes and static reponses of moduled and one-bodied floating structures. Equations for calculating natural modes and static responses are formulated by finite element method and the natural modes are solved by subspace iteration method. A floating parking place whose length is 120 m and width 60 m is considered as an example structure.

• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1114-1124
• /
• 2017

Acceleration/preconditioning strategies available in the SCEPTRE radiation transport code are described. A flexible transport synthetic acceleration (TSA) algorithm that uses a low-order discrete-ordinates ($S_N$) or spherical-harmonics ($P_N$) solve to accelerate convergence of a high-order $S_N$ source-iteration (SI) solve is described. Convergence of the low-order solves can be further accelerated by applying off-the-shelf incomplete-factorization or algebraic-multigrid methods. Also available is an algorithm that uses a generalized minimum residual (GMRES) iterative method rather than SI for convergence, using a parallel sweep-based solver to build up a Krylov subspace. TSA has been applied as a preconditioner to accelerate the convergence of the GMRES iterations. The methods are applied to several problems involving electron transport and problems with artificial cross sections with large scattering ratios. These methods were compared and evaluated by considering material discontinuities and scattering anisotropy. Observed accelerations obtained are highly problem dependent, but speedup factors around 10 have been observed in typical applications.

• 전자공학회논문지
• /
• 제51권5호
• /
• pp.177-187
• /
• 2014

본 논문에서는 고속 부공간 추적 기법인 FAPI (Fast Approsimated Power Iteration)에 GVFF RLS (Gradient-based Variable Forgetting Factor Recursive Least Square Error)를 적용한 GVFF FAPI 를 제안한다. 기존의 FAPI는 신호의 공분산 행렬을 추정하기 위해 고정 망각 인자를 사용하기에, 부공간이 지속적으로 변하는 비정재 환경에 적용하기 여려운 단점이 있다. 이러한 문제점을 해결하기 위해, GVFF FAPI는 개선된 MSE (Mean Square Error)의 분석으로부터 유도된 MSE의 기울기 기반의 시변 망각 인자를 사용한다. 또한 GVFF RLS의 망각 인자 업데이트 식을 개선하여 부공간이 지속적으로 변하는 비정재 환경에서 부공간 에러를 줄인다. 개선된 망각 인자 업데이트 식은 MSE의 기울기가 양수이면 망각 인자를 빠르게 감소하게 하고 MSE의 기울기가 음수이면 망각 인자를 천천히 증가시킨다. 모의실험을 통해서 도래각이 지속적으로 변하는 환경에서 GVFF FAPI 알고리즘이 기존의 FAPI 알고리즘보다 작은 부공간 에러를 가지는 것을 보이고, 추적된 부공간을 도래각 추정기법에 적용하였을 때 추적된 도래각의 RMSE (Root Mean Square Error)가 더 작은 것을 확인한다.

• International Journal of Aeronautical and Space Sciences
• /
• 제17권2호
• /
• pp.184-194
• /
• 2016

In existing identification methods for on-orbit spacecraft, such as eigensystem realization algorithm (ERA) and subspace method identification (SMI), singular value decomposition (SVD) is used frequently to estimate the modal parameters. However, these identification methods are often used to process the linear time-invariant system, and there is a lower computation efficiency using the SVD when the system order of spacecraft is high. In this study, to improve the computational efficiency in identifying time-varying modal parameters of large spacecraft, a faster recursive algorithm called fast approximated power iteration (FAPI) is employed. This approach avoids the SVD and can be provided as an alternative spacecraft identification method, and the latest modal parameters obtained can be applied for updating the controller parameters timely (e.g. the self-adaptive control problem). In numerical simulations, two large flexible spacecraft models, the Engineering Test Satellite-VIII (ETS-VIII) and Soil Moisture Active/Passive (SMAP) satellite, are established. The identification results show that this recursive algorithm can obtain the time-varying modal parameters, and the computation time is reduced significantly.

• 전산구조공학
• /
• 제4권1호
• /
• pp.73-82
• /
• 1991

판의 탄소성 좌굴문제는 판 구조의 해석과 설계시의 중요성으로 인하여 상당한 관심이 모아져 온 분야이다. 본 연구에서는 유한요소법에 의한 효율적인 탄소성 좌굴해석 프로그램을 개발하였다. 탄소성 강성행렬을 구성하기 위한 소성이론으로는 실험결과와 잘 일치하는 Stowell의 변형이론을 사용하였으며, 좌굴하중을 해석하기 위해서는 고유치해석에 의한 반복기법을 사용하였다. 고유치해석에서는 불필요한 고유치의 계산을 피할 수 있는 subspace반복기법을 사용하였다. 해석결과를 Stowell이 제시한 이론해와 Pride에 의한 실험결과와 비교하여 프로그램의 타당성을 보이고, 이를 이용하여 단순지지, 또는 고정된 경계조건에 대하여 일축 또는 이축응력이 작용되는 여러 경우에 대하여 좌굴하중을 구하였다. 또한, 탄소성 좌굴에 미치는 형상비의 영향을 검토하였다.

• Structural Engineering and Mechanics
• /
• 제65권1호
• /
• pp.69-79
• /
• 2018

This paper is devoted to two new techniques for free vibration analysis of concrete arch dam-reservoir systems. The proposed schemes are quadratic ideal-coupled eigen-problems, which can solve the originally non-symmetric eigen-problem of the system. To find the natural frequencies and mode shapes, a new special-purpose eigen-value solution routine is developed. Moreover, the accuracy of the proposed approach is thoroughly assessed, and it is confirmed that the new scheme is very accurate under all practical conditions. It is also concluded that both decoupled and ideal-coupled strategy proposed in the previous works can be considered as special cases of the current more general procedure.

• Structural Engineering and Mechanics
• /
• 제14권3호
• /
• pp.245-262
• /
• 2002

The free vibration analysis of stiffened laminated composite plates has been performed using the layered (zigzag) finite element method based on the first order shear deformation theory. The layers of the laminated plate is modeled using nine-node isoparametric degenerated flat shell element. The stiffeners are modeled as three-node isoparametric beam elements based on Timoshenko beam theory. Bilinear in-plane displacement constraints are used to maintain the inter-layer continuity. A special lumping technique is used in deriving the lumped mass matrices. The natural frequencies are extracted using the subspace iteration method. Numerical results are presented for unstiffened laminated plates, stiffened isotropic plates, stiffened symmetric angle-ply laminates, stiffened skew-symmetric angle-ply laminates and stiffened skew-symmetric cross-ply laminates. The effects of fiber orientations (ply angles), number of layers, stiffener depths and degrees of orthotropy are examined.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제13권4호
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.