• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

• Journal of IKEEE
• /
• v.2 no.2 s.3
• /
• pp.263-269
• /
• 1998

We use the slope bounded saturation nonlinear feedback element instead of relay to find ultimate gain and period of linear plant. Saturation nonlinear element reduces the high harmonics of plant output. The reduction of high harmonics improve the accuracy of describing function method used to find ultimate gain and period. We give a simple procedure to find ultimate gain and period with saturation nonlinear element. A PID controller design method with known time delay element is also given, which is very useful when oscillation is not occurred with nonlinear element.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.647-660
• /
• 1995

Given the tetrahedron-based octree approximation of a solid as described in part(I) of this thesis, in this part(II) a systematic procedure of 'boundary moving' is developed for the fully automatic generation of 3D finite element meshes. The algorithm moves some vertices of the octants near the boundary onto the exact surface of a solid without transforming the topology of octree leaf elements. As a result, the inner octree leaf elements can be used as exact tetrahedral finite element meshes. In addition, as a quality measure of a tetrahedral element, 'shape value' is propopsed and used for the generation of better finite elements during the boundary moving process.

• Structural Engineering and Mechanics
• /
• v.28 no.4
• /
• pp.461-477
• /
• 2008

An adaptive finite element method for analyzing two-dimensional and axisymmetric nonlinear elastic fracture mechanics problems with cracks is presented. The J-integral is used as a parameter to characterize the severity of stresses and deformation near crack tips. The domain integral technique, for which all relevant quantities are integrated over any arbitrary element areas around the crack tips, is utilized as the J-integral solution scheme with 9-node degenerated crack tip elements. The solution accuracy is further improved by incorporating an error estimation procedure onto a remeshing algorithm with a solution mapping scheme to resume the analysis at a particular load level after the adaptive remeshing technique has been applied. Several benchmark problems are analyzed to evaluate the efficiency of the combined domain integral technique and the adaptive finite element method.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.617-624
• /
• 1997

This Paper Presents an effective and simple procedure for the simulation of the motion of the solid-liquid interfacial boundary and the transient temperature field during phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual reciprocity boundary element method. The dual reciprocity boundary element approach provided in this paper is much simpler than the usual boundary element method applying a reciprocity principle and an available technique for dealing with domain integral of boundary element formulation simultaneously. The effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of an example with its semi-analytical or other numerical solutions where available.

• Transactions of Materials Processing
• /
• v.9 no.2
• /
• pp.128-139
• /
• 2000

Wrinkling is one of the major defects in sheet metal products and may be also attributable to the wear of the tool. The initiation and growth of wrinkles are influenced by many factors such as stress state, mechanical properties of the sheet material, geometry of the body, and contact condition. It is difficult to analyze the wrinkling initiation and growth considering the factors because the effects of the factors are very complex and the wrinkling behavior may show a wide variation for small deviations of the factors. In this study, the bifurcation theory is introduced for the finite element analysis of wrinkling initiation and growth. All the above mentioned factors are conveniently considered by the finite element method. The finite element formulation is based on the incremental deformation theory and elastic-plastic elements considering the planar anisotropy of the sheet metal. The proposed method is verified by employing a column buckling problem. And then, the initiation and growth of wrinkling in deep drawing of cylindrical cup are analyzed.

• Transactions of Materials Processing
• /
• v.9 no.6
• /
• pp.590-597
• /
• 2000

A finite element formulation lot the simulation of axisymmetric sheet hydroforming is proposed, and an implicit program is coded. In order to describe normal anisotropy of steel sheet, Hill's non-quadratic yield function (Hill, 1979) is employed. Frictional contacts among sheet surface, rigid tool surface, and flexible hydrostatic pressure are considered using mesh normal vectors based on finite element of the sheet. Applied hydraulic pressure is also considered as a function of forming rate and time and treated as an external loading. The complete set of the governing relations comprising equilibrium and interfacial equations is approximately linearized for Newton-Raphson algorithm. In order to verify the validity of the developed finite element formulation, the axisymmetric bulge test is simulated. Simulation results are compared with other FEM results and experimental measurements and showed good agreements. In axisymmetric hydroforming processes of a disk cover, formability changes are observed according to the hydraulic pressure curve changes.

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

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.61 no.8
• /
• pp.1123-1129
• /
• 2012

This paper presents the fast steady state analysis using time-stepping finite element method for a high power three-phase transformer. The high power transformer spends huge computational cost of the time-stepping finite element method. It is because that the high power transformer requires a lot of time to reach steady state by its large inductance component. In order to reduce computational cost, in this paper, the adaptive time-step control algorithm combined with the embedded 2nd 4th singly diagonally implicit Runge-Kutta method and the analysis strategy using variation of the winding resistance are studied, and their numerical results are compared with those from the typical time-stepping finite element method.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.37 no.1
• /
• pp.18-24
• /
• 1988

Up to date, C-0 Finite Element Method which is the means for analyzing electric machinery can not be got the precision magnetic flux density because the magnetic flus density has the discontinuity in the interelement. To supplement this defect, we propose the C-1 finite element method of 9 D.O.F. in this paper. In this method, the vector potential and the magnetic flux density are continuous on the interelement and direction derivative of potential would be an unknown value. We developed the algorithm to apply this method. For examining the utility, we applied this method to analytic model and compared with the result of C-0 Finite Element Method using linear element.