• Transactions of Materials Processing
• /
• v.4 no.3
• /
• pp.214-221
• /
• 1995

In the paper, we have proposed a new method to determine the initial billet for the forged products using a function approximation in the neural network. The architecture of neural network is a three-layer neural network and the back propagation algorithm is employed to train the network. By utilizing the ability of function approximation of a neural network, an optimal billet is determined by applying the nonlinear mathematical relationship between the aspect ratios in the initial billet and the final products. The amount of incomplete filling in the die is measured by the rigid-plastic finite element method. The neural network is trained with the initial billet aspect ratios and those of the unfilled volumes. After learning, the system is able to predict the filling regions which are exactly the same or slightly different to the results of finite element simulation. This new method is applied to find the optimal billet size for the plane strain rib-web product in cold forging. This would reduce the number of finite element simulation for determining the optimal billet size of forging product, further it is usefully adapted to physical modeling for the forging design.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Structural Engineering and Mechanics
• /
• v.25 no.4
• /
• pp.467-480
• /
• 2007

Two methods were developed for analyzing problems with two adjacent sub-domains modeled by different kinds of elements in finite element method. Each sub-domain can be defined independently without the consideration of equivalent division with common nodes used for the interface. These two methods employ an individual interface to accomplish the compatibility. The MLSA method uses the moving least square approximation which is the basic formulation for Element Free Galerkin Method to formulate the interface. The displacement field assumed by this method does not pass through nodes on the common boundary. Therefore, nodes can be chosen freely for this method. The results show that the MLSA method has better approximation than traditional methods.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.183-189
• /
• 2014

In this paper, a polyhedral element is presented to solve three-dimensional problems by developing shape functions based on Wachspress coordinates and moving least square approximation. A subdivision of polyhedrons into tetrahedral domains is performed for the construction of shape functions of polyhedral elements, and numerical integration of the weak form is carried out consistently over the tetrahedral domains. The weight functions for moving least square approximation are defined by solving Laplace equation with boundary values based on Wachspress coordinates on polyhedral element faces. Polyhedral elements presented in this paper have similar properties to conventional finite element regarding the continuity, the completeness, the node-element connectivity and the inter-element compatibility. Numerical examples show the effectiveness of the present method for solving three-dimensional problems using polyhedral elements.

• Kyungpook Mathematical Journal
• /
• v.47 no.4
• /
• pp.529-548
• /
• 2007

The purpose of this article is to find a relation between the finite difference method and the boundary element method, and propose a new approach deriving a discrete approximation formula as like that of the finite difference method for harmonic functions. We develop a discrete approximation formula on a uniform grid based on the boundary integral formulations. We consider three different boundary integral formulations and derive one discrete approximation formula on the uniform grid for the harmonic function. We show that the proposed discrete approximation formula has the same computational molecules with that of the finite difference formula for the Laplace operator ${\nabla}^2$.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.1
• /
• pp.67-77
• /
• 2000

In this paper we consider finite element mathods for nonlinear parabolic problems defined in ${\Omega}{\subset}R^d$ ($d{\leq}4$). A new initial approximation is taken. Optimal order error estimates in $L_p$ for $2{\leq}p{\leq}{\infty}$ are established for arbitrary order finite element. One order superconvergence in $W^{1,p}$ for $2{\leq}q{\leq}{\infty}$ are demonstrated as well.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.2
• /
• pp.85-97
• /
• 2002

In this paper, finite volume element methods for nonlinear parabolic problems are proposed and analyzed. Optimal order error estimates in $W^{1,p}$ and $L_p$ are derived for $2{\leq}p{\leq}{\infty}$. In addition, superconvergence for the error between the approximation solution and the generalized elliptic projection of the exact solution (or and the finite element solution) is also obtained.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.1
• /
• pp.1-13
• /
• 2012

In this paper we present an a posteriori error estimator for the stabilized P1 nonconforming finite element method of the linear elasticity problem based on a nonsymmetric H(div)-conforming approximation of the stress tensor in the first-order Raviart-Thomas space. By combining the equilibrated residual method and the hypercircle method, it is shown that the error estimator gives a fully computable upper bound on the actual error. Numerical results are provided to confirm the theory and illustrate the effectiveness of our error estimator.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.4
• /
• pp.249-260
• /
• 2008

In this work we propose a mixed finite volume method for the Signorini problem which are based on the idea of Keller's finite volume box method. The triangulation may consist of both triangles and quadrilaterals. We choose the first-order nonconforming space for the scalar approximation and the lowest-order Raviart-Thomas vector space for the vector approximation. It will be shown that our mixed finite volume method is equivalent to the standard nonconforming finite element method for the scalar variable with a slightly modified right-hand side, which are crucially used in a priori error analysis.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.324-329
• /
• 2001

A novel method for non-matching interfaces on the boundaries of the finite elements in partitioned domains is presented by introducing interface elements in this paper. The interface element method (IEM) satisfies the continuity conditions exactly through interfaces without recourse to the Lagrange multiplier technique. The moving least square (MLS) approximation in the present study is implemented to construct the shape functions of the interface elements. Alignment of the boundaries of sub-domains in the MLS approximation and integration domains provides a consistent numerical integration due to one form of rational functions in an integration domain. The compatibility of displacements on the boundaries of the finite elements and the interface elements is always preserved in this method, and the completeness of the shape functions of the interface elements guarantees the convergence of numerical solutions. The numerical examples show that the interface element method is a useful tool for the analysis of a partitioned system and for a global-local analysis.