• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.743-749
• /
• 2007

Algebraic substructuring (AS) is a state-of-the-art method in eigenvalue computations, especially for large size problems, but, originally, it was designed to calculate only the smallest eigenvalues. In this paper, an updated version of AS is proposed to calculate the interior eigenvalues over a specified range by using a shift value, which is referred to as the shifted AS. Numerical experiments demonstrate that the proposed method has better efficiency to compute numerous interior eigenvalues for the finite element models of structural problems than a Lanczos-type method.

• Communications of the Korean Mathematical Society
• /
• v.15 no.1
• /
• pp.133-142
• /
• 2000

We consider some numerical solution methods for equilibrium equations Af + E$^{T}$ λ = r, Ef = s. Algebraic problems of this form evolve from many applications such as structural optimization, fluid flow, and circuits. An important approach, called the force method, to the solution to such problems involves dimension reduction nullspace computation for E. The purpose of this paper is to investigate the substructuring method for the solution step of the force method in the context of the incompressible fluid flow. We also suggests some iterative methods based upon substructuring scheme..

• Structural Engineering and Mechanics
• /
• v.61 no.6
• /
• pp.807-816
• /
• 2017

In this paper, we introduce an efficient new model reduction method, named the automated static condensation method, which is developed for the local analysis of large finite element models. The algebraic multilevel substructuring procedure is modified appropriately, and then applied to the original static condensation method. The retained substructure, which is the local finite element model to be analyzed, is defined, and then the remaining part of the global model is automatically partitioned into many omitted substructures in an algebraic perspective. For an efficient condensation procedure, a substructural tree diagram and substructural sets are established. Using these, the omitted substructures are sequentially condensed into the retained substructure to construct the reduced model. Using several large practical engineering problems, the performance of the proposed method is demonstrated in terms of its solution accuracy and computational efficiency, compared to the original static condensation method and the superelement technique.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.135-138
• /
• 2009

Most of the studies use model order reduction for low frequency (LF) response analysis due to their high computational efficiency. In LF response analysis, one of model order reduction, algebraic substructuring (AS) retains all LF modes when using the modal superposition. However, in mid-frequency (MF) response analysis, the LF modes make very little contribution and also increase the number of retained modes, which leads to loss of computational efficiency. Therefore, MF response analysis should consider low truncated modes to improve the computational efficiency. The current work is focused on improving the computational efficiency using a AS and a frequency sweep algorithm. Finite element simulation for a MEMS resonator array showed that the performance of the presented method is superior to a conventional method.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.467-472
• /
• 2007

A modern MEMS resonator is a micro-scale structure operated over a high frequency range. In order to predict its resonant behavior in a design process, High-frequency response analysis (Hi-FRA) is demanded. Algebraic substructuring (AS) is known as a fast numerical technique to construct an eigenspace for FR and frequency sweep (FS) algorithm efficiently solves the frequency response system projected on the eigenspace. However, the existing FS algorithm using AS is developed for low-FRA, say over the range 1Hz-2KHz. In this work, we extend the FS algorithm using AS for FRA over an arbitrary frequency range. Therefore, it can be efficiently applied to systems operated at a high frequency, say over the range 230MHz-250MHz. The success of the proposed method is demonstrated by Hi-FRA of a checkerboard resonator.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.83-90
• /
• 2004

Substructuring methods are usually used in finite element structural analyses. In this study a multi-level substructuring algorithm is developed and proposed as a possible candidate for incompressible fluid solves. Finite element formulation for incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et.al.[5]. The present algorithm consists of four stages such as a gathering stage, a condensing stage, a solving stage and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At highest level, each subdomain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each subdomain has been replaced by a sequential static condensation. The global algebraic system arising feom the assembly of each subdomains is solved using Conjugate Gradient Squared(CGS) method. In this case, pre-conditioning techniques usually accompanied by iterative solvers are not needed.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Structural Engineering and Mechanics
• /
• v.65 no.4
• /
• pp.359-368
• /
• 2018

Dynamic analysis of complex and large structures may be costly from a numerical point of view. For coupled vibroacoustic finite element models, the importance of reducing the size becomes obvious because the fluid degrees of freedom must be added to the structural ones. In this paper, a component mode synthesis method is proposed for large vibroacoustic interaction problems. This method couples fluid subdomains and dynamical substructuring of Craig and Bampton type. The acoustic formulation is written in terms of the velocity potential, which implies several advantages: coupled algebraic systems remain symmetric, and a potential formulation allows a direct extension of Craig and Bampton's method to acoustics. Those properties make the proposed method easy to implement in an existing finite element code because the local numerical treatment of substructures and fluid subdomains is undifferentiated. Test cases are then presented for axisymmetric geometries. Numerical results tend to prove the validity and the efficiency of the proposed method.