• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.3
• /
• pp.32-41
• /
• 1995

In the analysis of metal forming processes by the finite element method, there are many numerical instabilities such as element locking, hourglass mode and shear locking. These instabilities may have a bad effect upon accuracy and convergence. The present work is concerned with improvement of stability and efficiency in two-dimensional rigid-plastic finite element method using various type of elemenmts and numerical intergration schemes. As metal forming examples, upsetting and backward extrusion are taken for comparison among the methods: various element types and numerical integration schemes. Comparison is made in terms of stability and efficiency in element behavior and computational efficiency and a new scheme of adaptive directional reduced integration is introduced. As a result, the finite element computation has been stabilized from the viewpoint of computational time, convergency, and numerical instability.

• Structural Engineering and Mechanics
• /
• v.29 no.5
• /
• pp.469-488
• /
• 2008

Finite element methods have often been used for structural analyses of various mechanical problems. When finite element analyses are utilized to resolve mechanical systems, numerical uncertainties in the initial data such as structural parameters and loading conditions may result in uncertainties in the structural responses. Therefore the initial data have to be as accurate as possible in order to obtain reliable structural analysis results. The typical finite element method may not properly represent discrete systems when using uncertain data, since all input data of material properties and applied loads are defined by nominal values. An interval finite element analysis, which uses the interval arithmetic as introduced by Moore (1966) is proposed as a non-stochastic method in this study and serves a new numerical tool for evaluating the uncertainties of the initial data in structural analyses. According to this method, the element stiffness matrix includes interval terms of the lower and upper bounds of the structural parameters, and interval change functions are devised. Numerical uncertainties in the initial data are described as a tolerance error and tree graphs of uncertain data are constructed by numerical uncertainty combinations of each parameter. The structural responses calculated by all uncertainty cases can be easily estimated so that structural safety can be included in the design. Numerical applications of truss and frame structures demonstrate the efficiency of the present method with respect to numerical analyses of structural uncertainties.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.1
• /
• pp.22-26
• /
• 2002

In this paper, a finite element method is applied to the numerical calculation of the harbor calmness. The mild stop equation as the basic equation is used. The key of this model is that the bottom friction and boundary absorption are imposed. A numerical result is presented and compared with the results obtained from the other numerical analysis. These results are in very well agreement. This method calculating the calmness can be broadly utilized making the new design of harbor and fishing port in the future.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.32 no.E
• /
• pp.59-66
• /
• 1990

Abstract The formulation of the probabilistic finite element method was briefly reviewed. The method was implemented into a computer program for frame analysis which has the same analogy as finite element analysis. Another program for Monte Carlo simulation of finite element analysis was written. Two sample structures were assumed and analized. The characteristics of the second moment statistics obtained by the probabilistic finite element method was examined through numerical studies. The applicability and limitation of the method were also evaluated in comparison with the data generated by Monte Carlo simulation.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.5
• /
• pp.141-153
• /
• 1999

A finite element procedure for the analysis of rubber-like hyperelastic material is developed. The volumetric incompressiblity conditions of the rubber deformation is included in the formulation by using penalty method. In this paper, the behavior of the rubber deformation is represented by hyperelastic constitutive relations based on a generalized Mooney-Rivlin model. The principle of virtual work is used to derive nonlinear finite element equation for the large displacement problem and presented in total-Lagrangian description. The finite element procedure using analytic differentiation resulted in very close solution to the result of the well known commercial packages NISAII AND ABAQUS. Numerical tests show that the results from the numerical differentiation method coincide very well with those from the analytic method and the well known commercial packages in static analysis. The convergency of rubber usingν iteration method is also discussed.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.4
• /
• pp.19-29
• /
• 1982

Galerkin's finite element method is applied to a two-dimensional heat convection-diffusion problem arising in the hydrodynamic lubrication of thrust bearings used in naval vessels. A parabolized thermal energy equation for the lubricant, and thermal diffusion equations for both bearing pad and the collar are treated together, with proper juncture conditions on the interface boundaries. it has been known that a numerical instability arises when the classical Galerkin's method, which is equivalent to a centered difference approximation, is applied to a parabolic-type partial differential equation. Probably the simplest remedy for this instability is to use a one-sided finite difference formula for the first derivative term in the finite difference method. However, in the present coupled heat convection-diffusion problem in which the governing equation is parabolized in a subdomain(Lubricant), uniformly stable numerical solutions for a wide range of the Peclet number are obtained in the numerical test based on Galerkin's classical finite element method. In the present numerical convergence errors in several error norms are presented in the first model problem. Additional numerical results for a more realistic bearing lubrication problem are presented for a second numerical model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• pp.3-10
• /
• 1997

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.10
• /
• pp.1590-1597
• /
• 2004

Three methods of design sensitivity analysis for structures such as numerical method, analytical method and semi-analytical method have been developed for the last three decades. Although analytical design sensitivity analysis can provide very exact result, it is difficult to implement into practical design problems. Therefore, numerical method such as finite difference method is widely used to simply obtain the design sensitivity in most cases. The numerical differentiation is sufficiently accurate and reliable fur most linear problems. However, it turns out that the numerical differentiation is inefficient and inaccurate in nonlinear design sensitivity analysis because its computational cost depends on the number of design variables and large numerical errors can be included. Thus the semi-analytical method is more suitable for complicated design problems. Moreover, semi-analytical method is easy to be performed in design procedure, which can be coupled with an analysis solver such as commercial finite element package. In this paper, implementation procedure fur the semi-analytical design sensitivity analysis outside of the commercial finite element package is studied and the computational technique is proposed for evaluating the partial differentiation of internal nodal force, so called pseudo-load. Numerical examples coupled with commercial finite element package are shown to verify usefulness of proposed semi-analytical sensitivity analysis procedure and computational technique for pseudo-load.

• Structural Engineering and Mechanics
• /
• v.20 no.4
• /
• pp.405-420
• /
• 2005

In this paper, a hybrid/mixed nonlinear shell element is developed in polar coordinate system based on Hellinger/Reissner variational principle and the large-deflection theory of plate. A numerical solution scheme is formulated using the hybrid/mixed finite element method (HMFEM), in which the nodal values of bending moments and the deflection are the unknown discrete parameters. Stability of the present element is studied. The large-deflection analyses are performed for simple supported and clamped circular plates under uniformly distributed and concentrated loads using HMFEM and the traditional displacement finite element method. A parametric study is also conducted in the research. The accuracy of the shell element is investigated using numerical computations. Comparisons of numerical solutions are made with theoretical results, finite element analysis and the available numerical results. Excellent agreements are shown.