• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.3
• /
• pp.337-344
• /
• 2000

The wavelet analysis technique is applied in the spectral domain to efficiently represent the multi-scale features of the impedance matrices. In this scheme, the 2-D quadtree decomposition (applying the wavelet transform to only the part of the matrix) method often used in image processing area is applied for a sparse moment matrix. CG(Conjugate-Gradient) method is also applied for saving memory and computation time of wavelet transformed moment matrix. Numerical examples show that for rectangular cylinder case the non-zero elements of the transformed moment matrix grows only as O($N^{1.6}$).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.4B
• /
• pp.377-384
• /
• 2009

In order to reduce computing time significantly for a large matrix in EFEM of linear waves propagation over ripple beds, each of which is approximated to a multi-step topography, a partition method is presented to calculate reflection coefficients. By use of 10 evanescent modes in the model, the most accurate numerical solutions have been obtained up to date, which show different behaviors of computed reflection coefficient in some cases against the existing results. Both computing time and memory of the present partition model for solving a large matrix are still so much demanding that it is needed to develop an efficient method.

• Communications for Statistical Applications and Methods
• /
• v.2 no.2
• /
• pp.404-413
• /
• 1995

2수준계 요인실험법에서 Kim(1992) 에 의해 균형배열을 이용하여 설계된 resolution V 포화균형부분실시법에서 추정량들의 공분산행렬을 계산하여 통계적 특성을 연구하였다. 이러한 부분실시법은 최소의 처리조합수를 가지고 주효과와 2인자 교호작용까지 분석할 수 있는 특징이 있다. 특히 본 논문에서는 인자의 수에 따라 설계가능한 8개의 부분실시법들간의 유사성과 통계적 효율성, 그리고 index number들의 변화에 따른 공분산행렬의 특성을 살펴보았다.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.129-136
• /
• 2005

The damped least-squares inversion has become a most popular method in finding the solution in geophysical problems. Generally, the least-squares inversion is to minimize the object function which consists of data misfits and model constraints. Although both the data misfit and the model constraint take an important part in the least-squares inversion, most of the studies are concentrated on what kind of model constraint is imposed and how to select an optimum regularization parameter. Despite that each datum is recommended to be weighted according to its uncertainty or error in the data acquisition, the uncertainty is usually not available. Thus, the data weighting matrix is inevitably regarded as the identity matrix in the inversion. We present a new inversion scheme, in which the data weighting matrix is automatically obtained from the analysis of the data resolution matrix and its spread function. This approach, named 'extended active constraint balancing (EACB)', assigns a great weighting on the datum having a high resolution and vice versa. We demonstrate that by applying EACB to a two-dimensional resistivity tomography problem, the EACB approach helps to enhance both the resolution and the stability of the inversion process.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1990.10a
• /
• pp.199-203
• /
• 1990

본 연구에서는 일반적인 형상을 갖는 회전축에서 축 방향 하중과 토오크를 각각 또는 동시에 고려할 수 있는 전달 행렬을 수식화하였고, 이에 의한 진 동 특성을 해석하기 위하여 강제진동 해석 및 자유진동 해석을 수행하였으 며, 해석결과에 대한 신빙성을 얻기 위하여 과도진동 해석 및 이산화 FFT를 통하여 검산기능을 추가하였다. 그리고 전달 행렬법에 새로운 고유치 수치해 법과 수치 적분법의 적용을 통하여 `수렴성과 안정성을 개선하였다.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.193-200
• /
• 1997

Even though there are many analysis methods of circular tanks on elastic foundation, the finite element method is widely used for that purpose. But the finite element method requires a number of memory spaces, computation time to solve large stiffness equations. In this study many the simplified methods(Analogy of Beam on Elastic Foundation, Foundation Stiffness Matrix, Finite Element Method and Transfer Matrix Method) are applied to analyze a circular tank on elastic foundation. By the given analysis methods, BEF analogy and foundation matrix method, the circular tank was transformed into the skeletonized frame structure. The frame structure was divided into several finite elements. The stiffness matrix of a finite element is related with the transfer matrix of the element. Thus, the transfer matrix of each finite element utilized the transfer matrix method to simplify the analysis of the tank. There were no significant difference in the results of two methods, the finite element method and the transfer matrix method. The transfer method applied to a circular tank on elastic foundation resulted in four simultaneous equations to solve completely.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.77-83
• /
• 1993

A modified Lanczos method combined with a substructure analysis technique was used for calculating natural frequencies and mode shapes of large structural systems. The method does not require generation and storage of stiffness and mass matrices of the entire structure. It only uses the stiffness and mass matrices of each substucture. No approximating assumptions are required other than the usual assumption of linear elastic system modelled by finite elements. Thus, natural frequencies and mode shapes for the finite element model employed are the same as those with or without the suhstructuring algorithm. To check the efficiency of the proposed method, first ten natural frequencies and the corresponding mode shapes of an open truss helicopter tail-boom structure are calculated by using it.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.4 no.3
• /
• pp.83-98
• /
• 2000

유연한 액체 저장탱크 내 유체의 부가질량 및 슬러싱 강성행렬을 도출하는 새로운 방법을 제시하였다. 비점성, 비압축성 이상유체를 표면 출렁임을 고려하여 경계요소법에 의하여 모델링하였다. 유체의 표면과 저장탱크 벽체의 접촉면과 같은 불연속 경계를 다루기 위해 특별한 과정을 도입하였다. 원통형 액체저장탱크의 지진응답해석에 적용하여 우수한 결과를 얻을 수 있음을 확인하였다.

• Proceedings of the KIEE Conference
• /
• 2000.07b
• /
• pp.825-827
• /
• 2000

자성체를 포함하는 자기 시스템을 해석하는데 있어 비선형과 히스테리시스(Hysteresis)는 매우 중요한 역할을 한다. 특히 재질의 히스테리시스 특성을 유한요소법(FEM)을 이용하여 계산하기 위해서 많은 방법들이 소개되었다. 단순 반복법이나 Fixed Point Technique(FPT), M-iteration 법. 뉴튼 랍슨 (Newton-Raphson) 법 등이 그 예이다. 이 방법들 중에서 뉴튼 랍슨법은 빠른 수렴 특성으로 가장 많이 사용되고 있다. 하지만 뉴튼-랍슨법을 이용하여 히스테리시스 재질을 해석할 때는 매 반복 계산 때마다 계 계수행렬(System Stiffness matrix)이 변화하기 때문에 요소의 수가 매우 많을 경우 역행렬을 계산하기 위한 시간이 많이 소요되는 단점이 있다. 특히 히스테리시스 해석의 경우에는 주로 time-step법을 이용하여 계산하므로 가장 시간이 많이 소요되는 행렬 계산 시간을 단축함으로써 전체 계산 시간을 크게 줄일 수 있다. 최근 비선형 해석에서 TLM(Transmission Line Modeling)법이 도입되어 비선형 해석 시의 계산 시간을 크게 단축할 수 있게 되었다. 본 논문에서는 비선형 해석에 적용된 TLM법을 히스테리시스 해석에 적용하는 방법을 새로 제안한다. TLM법은 뉴튼-랍슨법과 달리 각 반복 계산 때마다 계수행렬식이 변화하지 않고 단지 구동항만 변하기 때문에 행렬의 LU를 한 번 저장해 두면 forward와 backward substitution만 시행하면 된다. 따라서 요소의 수가 증가할 경우 TLM법을 사용하면 뉴튼-랍슨법에 비해 매우 큰 계산 이득을 얻을 수 있다. 본 논문에서는 TLM법을 히스테리시스에 적용하는 방법을 기술하고 간단한 모델에 이 방법을 적용하여 뉴튼-랍슨법과의 비교를 통해 TLM법의 효용성을 보인다.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.147-153
• /
• 2007

The conjugate gradient (CG) method is one of the most efficient algorithms for solving a linear system of equations. In addition to being used as a linear equation solver, it can be applied to a least-squares problem. When the CG method is applied to large-scale three-dimensional inversion of magnetotelluric data, two approaches have been pursued; one is the linear CG inversion in which each step of the Gauss-Newton iteration is incompletely solved using a truncated CG technique, and the other is referred to as the nonlinear CG inversion in which CG is directly applied to the minimization of objective functional for a nonlinear inverse problem. In each procedure we only need to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector, significantly reducing the computational requirements needed to do large-scale inversion.