• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1581-1593
• /
• 2001

An efficient variable-reordering method for finite element meshes is used and the effect of variable-reordering is investigated. For the element renumbering of unstructured meshes, Cuthill-McKee ordering is adopted. The newsy reordered global matrix has a much narrower bandwidth than the original one, making the ILU preconditioner perform bolter. The effect of variable reordering on the convergence behaviour of saddle point type matrix it studied, which results from P2/P1 element discretization of the Navier-Stokes equations. We also propose and test 'level(0) preconditioner'and 'level(2) ILU preconditioner', which are another versions of the existing 'level(1) ILU preconditioner', for the global matrix generated by P2/P1 finite element method of incompressible Navier-Stokes equations. We show that 'level(2) ILU preconditioner'performs much better than the others only with a little extra computations.

• Structural Engineering and Mechanics
• /
• v.52 no.6
• /
• pp.1099-1120
• /
• 2014

In this paper, reliability-based design optimization (RBDO) of structures is addressed. For this purpose, the global search and optimization capabilities of genetic algorithm (GA) are combined with the efficiency and reasonable accuracy of an advanced moment-based finite element reliability method. For performing RBDO, three variants of GA including a real-coded, a binary-coded and an improved binary-coded GA are developed. In these methods, GA performs (finite element) reliability analyses to evaluate reliability constraints. For truss structures which include finite element modeling, reliability constraints are evaluated using finite element reliability analysis. Response sensitivity required for finite element reliability analysis is obtained by direct differentiation method (DDM) rather than finite difference method (FDM). The proposed methods are examined within four standard test examples and real-world design problems. The results illustrate the superiority and efficiency of the improved binary-coded GA. Results also illustrate that DDM significantly reduces the computational cost and improves the efficiency of the optimization procedure.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.38 no.1
• /
• pp.22-28
• /
• 1989

This paper describes error estimate mothod for adaptive finite element analysis of two dimensional electrostatic and magnetostatic field problems. To estimate the local errors, divergence theorem is used for electrostatic field and Ampere's circuital law for magnetostatic field. To confirm the effectiveness of the proposed error estimators, adaptive finite element computations are performed using the proposed error estimators. The rates of convergence of global errors are comparable with those of existing adaptive finite element schemes which make use of field continuity conditions between element boundaries. This algorithm of error estimate can be easily implemented because of its simplicity. Especially, when the value of charge in electrostatic field and the value of current in magnetostatic field are to be figured out, this method is considerded to be preferable to other approaches.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.11
• /
• pp.1785-1795
• /
• 2001

A multi-layered brick element fur the finite element method is developed for analyzing the three-dim-ensionally layered composite structures subjected to both thermal and mechanical boundary conditions. The element has eight nodes with one degree of freedom for the temperature and three for the display-ements at each node, and can contain arbitrary number of layers with different material properties with-in the element; the conventional element should contain one material within an element. Thus the total number of nodes and elements, which are needed to analyze the multi-layered composite structures, can be tremendously reduced. In solving the global equation, a partitioning technique is used to obtain the temperature and the displacements which are caused by both the mechanical boundary conditions and temperature distributions. The results by using the developed element are compared wish the commercial package, ANSYS and the conventional finite element methods, and they are in good agreement. It is also shown that the Number of nodes and elements can be tremendously reduced using the element without losing the numerical accuracies.

• Structural Engineering and Mechanics
• /
• v.4 no.1
• /
• pp.1-8
• /
• 1996

A new error measure for a static finite element analysis is proposed. This error measure is a modification to the Zienkiewicz and Zhu energy norm. The new error estimator is a global error measure for the analysis and is independent of finite element model size and internal stresses, hence it is readily transportable to other error calculations. It is shown in this paper the the new error estimator also produces conservative error measurements, making it a suitable procedure to adopt in commerical packages.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1022-1027
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered metal bearing and flexible structures by using the finite element method and the mode superposition method considering the asymmetry of the gyroscopic effect and sintered metal bearing. The eigenvalues and eigenvectors are calculated by solving the eigenvalue problem and the adjoint eigenvalue problem by using the restarted Arnoldi iteration method. The decoupled equations of motion can be obtained from global finite element motion equations by using the orthogonal relation between the right eigenvectors and left eigenvectors. The decoupled equations of motion are used to analyze the unbalance response of a high speed polygon mirror scanner motor. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.10
• /
• pp.226-232
• /
• 2003

Difficulties in the finite element modeling of human spine are evaded by using a stiffness matrix element whose properties can be characterized from experimentally measured stiffness of functional spinal units. Relative easiness is in that inter-vertebral discs, ligaments, and soft tissues connecting vertebrae do not need to be modeled as they are. The remarkable coupling effect between distinct degrees of freedom induced by the geometric complexity can be accommodated without much effort. An idealized block model with simple geometry for vertebra is employed to assess the feasibility of this method. Analyses are performed in both levels of motion segment and spinal column, and the result is compared with that from detail model. As far as the global behavior of spine is concerned, the simplification is found not to aggravate inaccuracy only if sufficient experimental data is provided and interpreted properly.

• Structural Engineering and Mechanics
• /
• v.21 no.6
• /
• pp.677-698
• /
• 2005

In engineering practice, tubular connections are usually assumed pinned or rigid. Recent research showed that tubular joints may exhibit non-rigid behavior under axial or bending loads. This paper is concerned with establishing a new classification for tubular joints and investigating the effect of joint rigidity on the global behavior of CHS (Circular Hollow Section) lattice girders. Parametric formulae for predicting tubular joint rigidities are proposed, which are based on the finite element analyses through systematic variation of the main geometric parameters. Comparison with test results proves the reliability of these formulae. By considering the deformation patterns of respective parts of Vierendeel lattice girders, the boundary between rigid and semirigid tubular connections is built in terms of joint bending rigidity. In order to include characteristics of joint rigidity in the global structural analysis, a type of semirigid element which can effectively reflect the interaction of two braces in K joints is introduced and validated. The numerical example of a Warren lattice girder with different joint models shows the great effect of tubular joint rigidities on the internal forces, deformation and secondary stresses.

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

Tire belt durability is characterized by the initiation of a crack at the belt edge region, and the propagation of the crack until the ply is separated. Experimental methods have been used to analysis of the belt durability in the cord-reinforced rubber composite tires, but it takes much cost and time to make experiments. In this paper, a finite element method to analyze the fatigue life of a crack at the belt edge of tires is presented. The fatigue life is analyzed by using a three-dimensional Finite Element Modeling. This method includes a global-local finite element analysis to provide the detail necessary to model explicitly an internal crack and use of the J-integral for energy release rate evaluation.