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

A numerical method for actuating performance evaluation of LIPCA proposed using a finite element method. Fully coupled formulations for piezo-electric materials were introduced and 3-dimensional eight-node incompatible element was used. After verifying the developed code with typical examples, the center deflections of LIPCA were calculated and compared with the experimental result, which were in fairly agreement.

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

A finite element model for polarization switching in electro-mechanically coupled materials is proposed and applied to predict the switching behavior of a two-dimensional ferroelectric ceramic. A complicated micro-structure existing in the material is modeled as il continuum body and a simple 3 node triangle finite element with nodal displacement and voltage degrees of freedom is used for a finite element analysis. The elements use nonlinear constitutive equations, switching criterion and kinetic relation, fur representation of material response at strong electric and stress fields. The polarization state of the material is represented by internal variables in each element, which are updated at each simulation step based on the proposed constitutive equations. The model reproduces strain and electric displacement hysteresis loops observed in the material.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.5
• /
• pp.120-127
• /
• 1998

Differential inequalities occurring in problems of obstacle contact problems are recast into variational inequalities and analyzed by finite element methods. A new a posteriori error estimator, which is essential in adaptive finite element method, is introduced to capture the errors in finite element approximations of these variational inequalities. In order to construct a posteriori error estimates, saddle point problems are introduced using Lagrange parameters and upper bounds are provided. The global upper bound is localized by a special mixed formulation, which leads to upper bounds of the element errors. A numerical experiment is performed on an obstacle contact problem to check the effectivity index both in a local and a global sense.

• Journal of Mechanical Science and Technology
• /
• v.14 no.4
• /
• pp.401-407
• /
• 2000

An elastic-plastic finite element analysis is performed to investigate detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The finite element analysis performed under plane stress using 4-node isoparametric elements can predict fatigue crack closure behavior. The mesh of constant element size along crack surface can not predict the opening level of fatigue crack. The crack opening level for the constant mesh size increases linearly from initial crack growth. The crack opening level for variable mesh size, is almost flat after crack tip has passed the monotonic plastic zone. The prediction of crack opening level using the variable mesh size proportioning the reversed plastic zone size with the opening stress intensity factors presents a good agreement with the experimental data regardless of stress ratios.

• Steel and Composite Structures
• /
• v.7 no.3
• /
• pp.241-251
• /
• 2007

The purpose of this study is to calculate the torsional rigidity of arbitrarily shaped composite sections on the basis of hybrid finite element approach. An analogy is used between the torsion problem and deformation of a plate which exhibits only shear behavior. In the analysis a simple hybrid finite element based on Hellinger-Reissner functional is presented and a set of numerical examples are performed to demonstrate and asses the performance of the developed element in practical applications.

• Journal of Welding and Joining
• /
• v.18 no.4
• /
• pp.28-37
• /
• 2000

In this paper the LNG vessel of the Moss type which is capable of lifting 15,261 tons is investigated in the view point of the pressure vessel preliminary design using the finite element method. The Pressure vessel design is based on the equivalent stress levels due to the internal pressure. The finite element model of the spherical pressure vessel is configured using 4 noded quadrilateral shell element. The finite element analysis program NASTRAN and ANSYS 5.5are implemented. The design is compared with the three kinds of the boundary condition : first, where the equator of the pressure vessel is fixed, and where the top and is fixed, and, the bottom end is fixed, respectively. A comparison is presented between the results obtained by the finite element model and by the prototype production model. Additionally just below position(case 1 & case 2) of equator ring was carried out by using ANSYS 5.5. The results show that the vessel design based on the stress is acceptable at the preliminary design.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.4
• /
• pp.90-101
• /
• 1997

A hybrid element method is presented for two-dimensional diffraction problem of water waves. In this method, only a limited fluid domain close to irregular bodies is discretized into conventional finite elements, while the remaining infinite domain is treated as one element with analytical representations of high accuracy. A finite element grid is automatically generated by using Dealunay triangulation based on the Bowyer's algorithm and a linear system of equations is approximately solved with the ILU-CGS algorithm. To validate the present scheme, Computational results are compared with the existing experimental data and other numerical solutions.

• 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 KSME Conference
• /
• 2003.11a
• /
• pp.1286-1290
• /
• 2003

This paper presents a newly implemented parallel finite element procedure for contact-impact problems. Three sub-algorithms are includes in the proposed parallel contact-impact procedure, such as a parallel Belytschko-Lin-Tsay (BLT) shell element generation, a parallel explicit time integration scheme, and a parallel contact search algorithm based on the master slave slide-line algorithm. The underlying focus of the algorithms is on its effectiveness and efficiency for inclusion in future finite element systems on parallel computers. Throughout this research, a prototype code, named GT-PARADYN, is developed on the IBM SP2, a distributed-memory computer. Some numerical examples are provided to demonstrate the timing results of the procedure, discussing the accuracy and efficiency of the code.

• Transactions of Materials Processing
• /
• v.9 no.4
• /
• pp.379-388
• /
• 2000

This study aims to examine the influence of experimental and numerical factors on the results of the test and finite element simulation to evaluate the formability of sheet metals. The stretch-forming test with a hemispherical punch is carried out to obtain the limiting dome height (LDH) and forming limit diagram (FLD) for several kinds of aluminium and steel sheet. The results of the LDH and FLD tests are analysed to find any correlation with the uniaxial tensile properties. It proves that the size of the prescribed grid has great influence on the measured value of strain. The finite element analysis of the punch stretching process is also carried out and the result is compared with the experimental data. The influence of the numerical parameters such as friction coefficient, element size and anisotropy model on the simulation results tms out to be very considerable.