• Tribology and Lubricants
• /
• v.4 no.2
• /
• pp.28-35
• /
• 1988

A general and efficient algorithm is proposed for the analysis of the frictionless elastic contact problems. It utilizes a simplex-type algorithm with a modified entry rule and incoporates finite element method to obtain flexibility matrices. The algorithmic solution is compared with the Hertzian solution for the contact problem between two cylinders to prove its accuracy and the contact problem between pin and piston rod is solved and compared with the numerical results of Frankavilla and Zienkiewicz to demonstrate the generality and effectiveness of the suggested algorithm. The contact problem between mating involute gear teeth at the worst load position is considered. The computed contact stress is smaller than the result of Hertz's theory applied to the contact between two kinematically equivalent discs and the contact area is larger than that of Hertz's theory.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1996.10a
• /
• pp.210-216
• /
• 1996

Spectral element method(SEM) has been recognized to provide accurate structural dynamic responses even at high frequency. The Duhamel's integral based SEM developed by the authors for the structures under dynamic distributed loads does not take into account the zero frequency effect and requires significant computational time due to the integration procedure of Duhamel's integral. Hence, in this paper, a new SEM algorithm is proposed. This new algorithm is based on the FEM-type procedure for the distributed loads and includes the zero frequency correction to improve the accuracy. Some numerical results are illustrated to prove the accuracy of present new SEM algorithm.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.168-175
• /
• 1996

This paper deals with the variable-node element for fluid flow and the adaptive h-version mesh refinement algorithm. The transient element has been formulated by the Galerkin approach in which the pressure term is replaced with the penalty function. The present element having variable mid-side node and is suitable for constructing a locally refined mesh avoiding the use of the highly distorted elements. A modified Gauss quadrature is needed to integrate the element matrices to solve the trouble associated with the discontinuity of derivatives of shape functions. Several numerical examples show that the proposed element can be effectively used in the h-version adapt ive mesh refinement

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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.6
• /
• pp.1773-1783
• /
• 1996

Finite element method(FEM) is one of the most popularly used method analyzing the dynamic behaviors of structures. But unless number of finite elements is large enough, the results from FEM some what different from exact analytical solutions, especially at high frequency range. On the other hand, as the spectral analysis method(SAM) deals directly with the governing equations of a structure, the results from this melthod cannot but be exact regardless of any frequency range. However, the SAM can be applied only to the case where a structure is subjected to the concentrated loads, despite a structure could be unddergone distributed loads more generally. In this paper, therefore, new spectral analysis algorithm is introduced through the spectral element method(SEM), so that it can be applied to anlystructures whether they are subjected to the concentrated loads or to the distributed loads. The results from this new SEM are compared with both the results from FEM and the exact analytical solutions. As expected, the results from new SEM algorithm are found to be almost identical to the exact analytical solutions while those from FEM are not agreed well with the exact analytical solutions as the mode number increases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.11
• /
• pp.865-873
• /
• 2013

In this paper, numerical methods for wrinkle prediction of thin membrane were studied by finite element analysis. Techniques using membrane and shell elements were applied for triangular membrane. In case of membrane element method, the wrinkling was accounted for by the wrinkle algorithm of property modification, which was implemented to ABAQUS as a user subroutine. In case of shell method, geometrically nonlinear post-buckling analysis was performed to obtain the wrinkle deformation explicitly. The wrinkling deformation was induced by seeding the mesh with a random geometric imperfection. The results were investigated focusing on the mesh convergence and the solution accuracy.

• Journal of Power System Engineering
• /
• v.19 no.1
• /
• pp.38-44
• /
• 2015

Various finite elements have been studied and developed to analyze a variety of structures in the finite element method(FEM). The transfer stiffness coefficient method(TSCM) is an effective algorithm for structural analysis but the structures which can be applied were limited. In this paper, a computational algorithm for the structural analysis of axisymmetric conical shells under axisymmetric loading is formulated using the finite element-transfer stiffness coefficient method(FE-TSCM). The basic concept of FE-TSCM is the combination of the modeling technique of FEM and the transfer technique of TSCM. The FE-TSCM has all the advantages of both FEM and TSCM. After carrying out the structural analysis of axisymmetric conical shells using FEM, FE-TSCM, and analytical method we compare the computational results of FE-TSCM with those of the other methods in terms of computational accuracy.

• Computers and Concrete
• /
• v.6 no.2
• /
• pp.109-132
• /
• 2009

This paper presents finite element impact analysis for the design of Structurally Dissipating Rock-shed (SDR), an innovative design of reinforced concrete rock-shed. By using an appropriate finite element impact algorithm, the SDR structure is modelled in a simplified but efficient way. The numerical results are firstly verified through comparisons with the results of the experiments recently realized by ESIGEC and TONELLO I.C. It is shown that, using this impact algorithm, it is possible to correctly predict the SDR structural behaviour under different rock-fall impact conditions. Moreover, the numerical results show that the slab centre is the critical impact location for reinforced concrete slab design. The impact analyses have thus been focused on the impacts at the slab centre for the SDR structural optimization. Several series of parametric studies have been carried out with respect to load cases and engineering parameters choices. These numerical results support the robustness of the new SDR concept, and serve to optimize SDR structure and improve its conventional engineering design, especially for ensuring the slab punching shear resistance.

• Journal of the Korean Magnetics Society
• /
• v.1 no.2
• /
• pp.79-84
• /
• 1991

A variational approach employing localized functional is presented to solve alternating magnetic field problems with open boundary. The functional used in the approach consists of the domain integral of finite element region only and the boundary integral of the interfacial boundary between the finite and infinite element regions. The boundary integral is obtained by transforming the infinite domain integral for the infinite element region into the interfacial boundary integral. The proposed algorithm is then applied to a simple two-dimensional problem where the analytic solutions are available. It is shown that the algorithm makes it possible to yield good agreements between the numerical and analytic solutions. and that it requires less computer storage memory and computation time than the conventional finite element method due to the reduction of the computing region.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2004.07b
• /
• pp.1222-1225
• /
• 2004

This paper presents multibit Sigma-Delta ADC using noise-shaped dynamic element matching(DEM). 5-bit flash ADC for multibit quantization in Sigma Delta modulator offers the following advantages such as lower quantization noise, more accurate white-noise level and more stability over single quantization. For the feedback paths consisting of DAC, the DAC element should have a high matching requirement in order to maintain the linearity performance which can be obtained by the modulator with a multibit quantizer. The DEM algorithm is implemented in such a way as to minimize additional delay within the feedback loop of the modulator Using this algorithm, distortion spectra from DAC linearity errors are shaped. Sigma Delta ADC achieves 82dB signal to noise ratio over 615H7z bandwidth, and 62mW power dissipation at a sampling frequency of 19.6MHz. This Sigma Delta ADC is designed to use 0.25um CMOS technology with 2.5V supply voltage and verified by HSPICE simulation.