• Journal of the korean Society of Automotive Engineers
• /
• v.12 no.6
• /
• pp.48-59
• /
• 1990

A finite element program was developed using line elements for simulating the stretch/draw forming operation of an arbitrarily-shaped plane-strain section. An implicit, incremental, updated Lagrangian formulation is employed, introducing a minimum plastic work path assumption for each time step. Geometric and material nonlinearities are also considered within each time step. The finite element equation is based on the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The membrane approximation is adopted under the plane stress assumption. The sheet material is assumed to obey a rigid-viscoplastic constitutive law. The developed program was tested in the die-tryout of typical automotive inner panels. In order to determine a single friction coefficient and boundary length, FEM results and measurements of thinning for a stretched section of final die were compared. After finding analysis parameters, the sheet forming operations of original and final die designs were simulated. Excellent agreement between measured and computed thickness strains was obtained and the developed program was able to identify die designs which were rejected during die tryout.

• Journal of applied mathematics & informatics
• /
• v.37 no.5_6
• /
• pp.357-382
• /
• 2019

In this paper, we consider boundary value problem for a weakly coupled system of two singularly perturbed differential equations of convection diffusion type with discontinuous source term. In general, solution of this type of problems exhibits interior and boundary layers. A numerical method based on streamline diffusiom finite element and Shishkin meshes is presented. We derive an error estimate of order $O(N^{-2}\;{\ln}^2\;N$) in the maximum norm with respect to the perturbation parameters. Numerical experiments are also presented to support our theoritical results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.11
• /
• pp.7-14
• /
• 2005

A subdomain-interface finite element method is suggested to solve a class of fully- coupled thermomechanical problems with contact boundaries. The penalty method is used for connecting subdomains that satisfy interface compatibility conditions. As a result, effective stiffness matrices are always positive definite, and computational efficiency can be improved to a considerable degree. Moreover, any complex-shaped domain can be divided into independently modeled subdomains without considering the conformity of meshes on interfaces. Using a computer code based on the present method, these advantageous features are shown through a set of numerical studies.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.5
• /
• pp.404-411
• /
• 2007

new domain/boundary decomposition method is suggested to perform efficient finite element analyses of contact problems. A penalty method is used for connecting an interface or contact interfaces with neighboring subdomains that satisfy continuity conditions. As a result, the derived effective stiffness matrices are always positive definite, and computational efficiency can be improved to a considerable degree. Moreover, any complex-shaped domain can be divided into independently modeled subdomains without considering the conformity of meshes along the interface. Using a computer code based on the present method, these advantageous features are confirmed through a set of numerical examples.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.5
• /
• pp.55-62
• /
• 2003

Subdomain-interface variational formulation is presented to solve a class of dynamic contact problems of composite structures. The penalty method is used for imposing inequality constraints on contact surfaces and for connecting finite element subdomains that satisfy interface compatibility conditions. As a result, any complex-shaped domain can be easily divided into independently modeled subdomains without considering the conformity of meshes on interfaces. Some advantageous features of the present method are shown through a set a numerical studies with a developed computer code.

• Computational Structural Engineering
• /
• v.2 no.4
• /
• pp.69-77
• /
• 1989

Until recently post-processing of finite element model has been heavily relied on expensive graphic peripheral devices. With the aid of inexpensive microcomputers, very economical post-processor graphics program called MICRO-POST has been developed. Model geometry or results of analysis for the unlimited meshes can be easily presented in a number of low-cost graphic devices. The paper presents the procedure obtaining the device-independent graphics, and the structure and functions of the program. It also describes efficient I/O scheme to overcome the memory limitation, and dialogue-type input technique to control the plot operation in an interactive manner. Through the post processing examples for the general purpose finite element programs, it demonstrates the usefulness of the program.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.3-17
• /
• 2003

This paper deals with the problems and their possible solutions in the development of finite element for analysis of shell. Based on these solution schemes, a series of flat shell elements are established which show no signs of membrane locking and other defects even though the coarse meshes are used. In the element formulation, non-conforming displacement modes are extensively used for improvement of element behaviors. A number of numerical tests are performed to prove the validity of the solutions to the problems involved in establishing a series of high performance flat shell elements. The test results reveal among others that the high accuracy and fast convergence characteristics of the elements are obtainable by the use of various non-conforming modes and that the ‘Direct Modification Method’ is a very useful tool for non-conforming elements to pass the patch tests. Furthermore, hierarchical and higher order non-conforming modes are proved to be very efficient not only to make an element insensitive to the mesh distortion but also to remove the membrane locking. Some numerical examples are solved to demonstrate the validity and applicability of the presented elements to practical engineering shell problems.

• Journal of Computational Design and Engineering
• /
• v.1 no.2
• /
• pp.88-95
• /
• 2014

Because the cost of performance testing using actual products is expensive, manufacturers use lower-cost computer-aided design simulations for this function. In this paper, we propose using hexahedral meshes, which are more accurate than tetrahedral meshes, for finite element analysis. We propose automatic hexahedral mesh generation with sharp features to precisely represent the corresponding features of a target shape. Our hexahedral mesh is generated using a voxel-based algorithm. In our previous works, we fit the surface of the voxels to the target surface using Laplacian energy minimization. We used normal vectors in the fitting to preserve sharp features. However, this method could not represent concave sharp features precisely. In this proposal, we improve our previous Laplacian energy minimization by adding a term that depends on multi-normal vectors instead of using normal vectors. Furthermore, we accentuate a convex/concave surface subset to represent concave sharp features.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.57-70
• /
• 2008

A high-order accurate Euler flow solver based on a discontinuous Galerkin finite-element method has been developed for the numerical simulations of blade-vortex interaction phenomena on unstructured meshes. A free vortex in freestream was investigated to assess the vortex-preserving property and the accuracy of the present flow solver. Blade-vortex interaction problems in subsonic and transonic freestreams were simulated by adopting a multi-level solution-adaptive dynamic mesh refinement/coarsening technique. The results were compared with those of other numerical and experimental methods. It was shown that the present discontinuous Galerkin flow solver can preserve the vortex structure for significantly longer vortex convection time and can accurately capture the complex unsteady blade-vortex interaction flows, including generation and propagation of acoustic waves.