• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1205-1219
• /
• 2011

In this work, a stabilized characteristic finite volume method for the time-dependent Navier-Stokes equations is investigated based on the lowest equal-order finite element pair. The temporal differentiation and advection term are dealt with by characteristic scheme. Stability of the numerical solution is derived under some regularity assumptions. Optimal error estimates of the velocity and pressure are obtained by using the relationship between the finite volume and finite element methods.

• Geotechnical Engineering
• /
• v.6 no.1
• /
• pp.25-34
• /
• 1990

The finite element and boundary element methods of underground analysis are both well established numerical techniques for determination of stress and displacement distributions at underground excavation. The finite element method presents antithetical advantages and limitations. Complex constitutive behaviour may be modelled, at the expense of numerical efficiency and, for infinite domain, inadequate representation of remote bounadry conditions. The inherent advantages of the boundary element method are the ease with which infinite domain problems may be analysed, and the efficiency of analysis typically associated with a boundary value solution procedure. Application of the method is limited by the requirements linear constitutive behaviour for the medium. A combined of the finite element and boundary element methods of underground analysis is shown to preserve the advantages of each procedure, and, eliminates their individual disadvantages. Procedures employed in this papers described combined FEBEM algorithm. Solutions of underground excavation verifying the performance of combined FEBEM code are compared with theoretical solution, boundary element solution and finite element solution.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.16 no.1
• /
• pp.29-36
• /
• 2008

In this paper, the orbital forming analysis of an automotive hub bearing was studied to predict forming performances using the explicit finite element method. To find an efficient solution technique for the orbital forming, axisymmetric finite element models and 3D solid element models were solved and numerically compared. The time scaling and mass scaling techniques were introduced to reduce the excessive computational time caused by small element size in case of the explicit finite element method. It was found from the numerical simulations on the orbital forming that the axisymmetric element models showed the similar results to the 3D solid element models in forming loads whereas the deformations at the inner race of bearing were quite different. Finally the strains at the inner race of bearing and the forming forces to the peen were measured for the same product of the numerical model by test, and were compared with the 3D solid element results. It was founded that the test results were in good agreements with the numerical ones.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.4
• /
• pp.199-206
• /
• 2018

In this study, numerical analysis is applied to a two - dimensional model for verifying the general finite element program, Abaqus' s extended finite element method(XFEM). The cohesive element model used in the existing research has a limitation in simulating the actual crack because of the disadvantage that the crack path should be predicted and the element should be inserted. For this reason, the extended finite element method(XFEM), which predicts the path of cracks based on the directionality and specificity of stress, is emerging as a new solution in crack analysis. The validity of the XFEM application was confirmed by comparing the cohesive element analysis with the XFEM analysis by applying the crack path to the self - evident two - dimensional model. Numerical analysis confirms stress distribution and stress specificity immediately before crack initiation and compares it with actual crack initiation path. Based on this study, it is expected that cracks can be simulated by performing actual crack propagation analysis of complex models.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.226-233
• /
• 1999

When we numerically model the bridge under seismic condition, the full model combining the super-structure and the sub-structure is considered for the more accurate results than the separate model. In this case, the super-structure is connected with the sub-structure by the elastic pad shoe that is difficult to model, because it has the three translational elastic stiffness and the three rotational elastic stiffness. The two-node General Link element is derived in finite element equation representing such a pad shoe, and it is verified by comparing the one General Link element model with the corresponding three legacy spring element model. It is easy to model the pad shoe, if the General Link finite element is used. And the seismic analysis result of the bridge full model structure, which is modeled with the General Link element, has been compared with the one of the separate model structure. The present study gives. more conservative result than that of the separate model, which does not consider the dynamic behaviour of the sub-structure.

• Structural Engineering and Mechanics
• /
• v.6 no.8
• /
• pp.883-899
• /
• 1998

The strength, deformation and buckling of a large engineering structure consisting of four ellipsoidal shells, two cylindrical shells with stiffening ribs and large holes, one conical shell and three pairs of large flanges under external pressure, self weight and heat sinks have been analysed by using two kinds of five different finite elements - four assumed displacement finite elements (shell element with curved surfaces, axisymmetric conical shell element with variable thickness, three dimensional eccentric beam element, axisymmetric solid revolutionary element) and an assumed stress hybrid element (a 3-dimensional special element developed by authors). The compatibility between different elements is enforced. The strength analyses of the top cover and the main vessel are described in the paper.

• Structural Engineering and Mechanics
• /
• v.26 no.1
• /
• pp.69-86
• /
• 2007

A four-noded plate bending quadrilateral (PBQ4) and an eight-noded plate bending quadrilateral (PBQ8) element based on Mindlin plate theory have been adopted for modeling the thick plates on elastic foundations using Winkler model. Transverse shear deformations have been included, and the stiffness matrices of the plate elements and the Winkler foundation stiffness matrices are developed using Finite Element Method based on thick plate theory. A computer program is coded for this purpose. Various loading and boundary conditions are considered, and examples from the literature are solved for comparison. Shear locking problem in the PBQ4 element is observed for small value of subgrade reaction and plate thickness. It is noted that prevention of shear locking problem in the analysis of the thin plate is generally possible by using element PBQ8. It can be concluded that, the element PBQ8 is more effective and reliable than element PBQ4 for solving problems of thin and thick plates on elastic foundations.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1994.10a
• /
• pp.229-234
• /
• 1994

The major objective of the present paper is to estabilish finite element simulation technique in order to further analyze the shearing process. For this the ductile fracture criterion and element kill method are used in the present work. It is postulated that a fracture initiation is based on the magnitude of local effective strain. The features of sheared surface are easily observed by the element kill method. The simulation results are compard with existing experimental results. It is found that the results of the present work are close agreement with the existing results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.2
• /
• pp.271-277
• /
• 1994

The fatigue crack opening behavior is analyzed using finite element method. Because extremely fine mesh subdivision is required when using constant stress constant strain triangular element, this study uses conventional two dimensional eight node isoparametric elements. Since plasitc zone size is similar to crack propagating length per each load cycle because of relatively large element size, a new analysis model that a crack propagates every two load cycle is suggested. the opening load and crack opening displacement can be obtained accurately by this method.

• Structural Engineering and Mechanics
• /
• v.37 no.6
• /
• pp.633-648
• /
• 2011

This paper investigates the aspect-ratio locking of the isoparametric 8-node quadrilateral (QUAD8) element. An important finding is that, if finite element solution is carried out with in exact arithmetic (i.e., with no truncation and round off errors), the locking tendency of the element is completely avoided even for aspect-ratios as high as 100000. The current finite element codes mostly use floating point arithmetic. Thus, they can only avoid this locking for aspect-ratios up to 100 or 1000. A novel method is proposed in the paper to avoid aspect-ratio locking in floating point computations. In this method, the offending terms of the strain-displacement matrix (i.e., $\mathbf{B}$-matrix) are multiplied by suitable scaling factors to avoid ill-conditioning of stiffness matrix. Numerical examples are presented to demonstrate the efficacy of the method. The examples reveal that aspect-ratio locking is avoided even for aspect-ratios as high as 100000.