- Home
- Documents
*Finite element analysis of compressible viscoplasticity using ... - 1 - Finite Element Analysis of...*

prev

next

out of 37

View

0Download

0

Embed Size (px)

HAL Id: hal-00574862 https://hal-mines-paristech.archives-ouvertes.fr/hal-00574862

Submitted on 10 Mar 2011

HAL is a multi-disciplinary open access archive for the deposit and dissemination of sci- entific research documents, whether they are pub- lished or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L’archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d’enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Finite element analysis of compressible viscoplasticity using a three-field formulation: Application to metal

powder hot compaction Michel Bellet

To cite this version: Michel Bellet. Finite element analysis of compressible viscoplasticity using a three-field formulation: Application to metal powder hot compaction. Computer Methods in Applied Mechanics and Engi- neering, Elsevier, 1999, 175 (1-2), pp.Pages 19-40. �10.1016/S0045-7825(98)00317-X�. �hal-00574862�

https://hal-mines-paristech.archives-ouvertes.fr/hal-00574862 https://hal.archives-ouvertes.fr

- 1 -

Finite Element Analysis

of Compressible Viscoplasticity

Using a Three-Field Formulation.

Application to Metal Powder Hot Compaction

Michel BELLET

Ecole des Mines de Paris

Centre de Mise en Forme des Matériaux (CEMEF)

UMR CNRS 7635

BP 207

06904 Sophia Antipolis Cedex, France

Abstract In the present study, a finite element model has been formulated to

simulate the hot forging stage in powder metallurgy manufacturing route. The

compacted material is assumed to obey a purely viscoplastic compressible flow rule. A

three-field formulation (velocity, volumetric strain rate and pressure) has been

developed. The associated three-dimensional finite element discretization is detailed. In

order to take advantage of an automatic remeshing procedure for linear tetrahedra, the

compatible P1+/P1/P1 element is used (4-node element plus additional degrees of

freedom and bubble interpolation for velocity). The complete model includes

thermomechanical coupling and friction. The formulation is validated versus an

analytic solution of uniaxial free compaction and applied to the hot forging of an

automotive connecting rod preform.

Keywords powder metallurgy - powder compaction - compressible viscoplasticity -

simulation - finite element method - mixed formulation - three-field formulation -

connecting rod

Comput. Methods Appl. Mech. Engrg. 175 (1999) 19-40

- 2 -

1 INTRODUCTION

Compaction of powdered materials has experienced great development in many

industries for production of various parts. The reasons for manufacturing a product

from powder are either economic or technical. They can be divided into three classes:

ability to elaborate materials which are difficult or even impossible to melt;

achievement of particular structures or properties (porosity, fine grains, isotropy,

purity); simplification of manufacturing routes, providing near net shape components

and raw material savings.

Densification of metal powders is obtained by compression and/or sintering. A typical

process is outlined as follows:

- Mixing of powder with some binder-lubricant. - Cold compaction step. This can be achieved either by punch compaction or by cold

isostatic pressing, typically up to 75 to 85% of full density.

- Dewaxing (binder removal) and optional sintering under specific conditions, to increase particle bonding.

- Hot compaction step: hot forging or hot isostatic pressing, in order to achieve full density.

Mathematical models of such a manufacturing process can be useful since they can

provide the production engineers with indication of possible underfill or final porosity,

evolution of density distribution and material flow throughout the process, tooling load

estimates, press size requirements, etc. They can also allow accurate and inexpensive

parametric studies of process variables.

The present work is essentially focused on the application of the finite element method

(FEM) to the analysis of hot forging of cold pre-compacted preforms. Hence we adress

here large viscoplastic deformations through generalized Newtonian fluid flow models,

including compressibility effects and coupling with heat transfer. Several authors have

already contributed to this research field. Among them, let us quote Im and Kobayashi

[20] who developed a FEM formulation for powdered metal forging an implemented it

into the bulk forming computation code DEFORM ® , for treating two-dimensional

axisymmetric and plane strain problems. The material is considered as a compressible

rigid-viscoplastic continuous medium. The authors have used a penalty-like

formulation. Such an approach penalizes the terms connected with volumetric strain in

both the porous and the fully dense regions of the preform. The only unknown of the

discretized problem is the nodal velocity field. A similar approach has been developed

by Barata Marques and Martins [4] who have clearly shown that some terms of the

finite element tangent stiffness matrix tends to become infinite as the material

- 3 -

approaches the dense state. A cut-off value of the relative density is then used in order

to limit the value of those terms and transform them into classical penalty terms

enforcing the incompressibility constraint, as in flow formulation. In order to avoid

locking, they have used a reduced Gauss integration. However, as pointed out by Jinka

and Lewis [22], such a formulation yields a wrong estimation of pressure, and

consequently stresses in the compacted material. These authors preferred using a mixed

velocity-pressure formulation and applied it successfully to the two-dimensional

analysis of hot isostatic pressing.

All these works were limited to two-dimensional analysis of compaction processes. The

first three-dimensional FEM simulation of powder forming has been developed by

Chenot et al. [7], but using a similar penalty-like flow formulation as [20, 4] and 8-node

linear hexahedral elements. Following this preliminary attempt to model metal

compaction in three-dimensions, the objective of the present study is to set out a new

three-dimensional mixed formulation, which would be easily implemented in

FORGE3 ® , a three-dimensional finite element code initially developed for non steady-

state large transformations of pure viscoplastic material [9]. In FORGE3 ® , the non

linear equilibrium equations are solved for the primitive variables velocity and pressure

using tetrahedral elements of P1+/P1 type. This permits the use of an automatic

remeshing procedure [10] and iterative solvers such as the preconditioned general

minimum residual method. Such solvers give rise to efficient parallelization [11].

Recently, elastic-viscoplastic constitutive equations have been implemented in the code

[2]. As part of these developments, it is required that the aimed compressible

formulation for hot powder compaction should be cast in the same type of element in

order to be consistent with the general software environment.

Hence, in the present paper, we shall first recall the governing equations of

compressible viscoplastic continuous media. The finite element formulation of the

mechanical problem is discussed, regarding particularly the choice of primitive

variables. We justify the choice of a three-field formulation and give details about its

finite element discretization using tetrahedral elements. The heat transfer problem and

its coupling with the mechanical one are also presented. This new formulation is

implemented in the computation code FORGE3 ® and is validated by comparison with

analytical solutions derived in the case of uniaxial free compression test. Finally, an

example of application to the hot forging of a connecting rod preform is presented.

2 GOVERNING EQUATIONS

2.1 Material Constitutive Equations

- 4 -

The powdered metal is considered as an average continuous medium. It is characterised

by its local relative density rρ , defined as the ratio of the apparent specific mass by the specific mass of the dense metal. Assuming material isotropy, the theoretical concepts

of plasticity can be extended to such an average medium, as proposed by Green [17], by

defining the equivalent stress as:

2/1

2)(Tr: 2

3

+=σ σss fc (1)

where σσσσ is the Cauchy stress tensor, s is the deviatoric stress tensor and, using the Einstein summation convention,

ijij ss=ss : iiσ=σTr (2)

Coefficients c and f are functions of the relative density rρ . They verify c(1) = 1 and f(1) = 0 (fig. 1), so that the classical definition of the von Mises equivalent stress is

obtained for the dense material.

Assuming an associated plasticity flow rule, i