(Jong Won Baek)
(Seong Jin Park)
*
Copyright © 2019 The Korean Institute of Metals and Materials
Key words(Korean)
compaction, FEM simulation, spherical powder, pressure transmission
1. INTRODUCTION
The Powder Metallurgy (PM) process is one of the more widely used methods for manufacturing
complex parts in industry. The advantages of this manufacturing technique is that
it minimizes material loss and weight [1]. The overall procedure of PM process is as shown in Fig 1.
The first stage of the PM process is mixing to prepare a compound powder with an additive
material such as a lubricant. Before mixing, sieving process is conducted to obtain
the desired powder size. Next, a compaction process is carried out to manufacture
green compacts with a desired shape. Finally, the green compacts are densified in
a sintering stage at high temperature. In the compaction stage, the final pressure
(compaction load) is determined, but the force changes continuously during the process.
The force change is influenced by the contact angle between the powder and various
other conditions. However, it is very hard to determine pressure changes using compaction
experiments, because it is impossible to measure the contact angle and contact area.
In this case, FEM simulation is a very useful method for examining pressure transmission
during the compaction process.
Most prior studies have focused on optimizing the compaction process conditions and
geometries to obtain full density [2-10]. Kwon et al. [11] and Lee et al. [12] carried out simulations on a powder compaction process for complex shapes such as
automobile parts, to predict their densities. Also, FEM has been widely used to examine
the compaction behavior of metal/ceramic materials [13-18]. However, the powders were considered to be continuum (bulk) green compacts, so
particle-particle interaction was not modelled.
The ‘particulate state’ of matter is peculiar. Powders consist of solid particles
yet exhibit fluid-like characteristics. They assume the shape of the vessel into which
they are ‘poured’ and are ‘compressible’. The consolidation of powders necessarily
involves the transmission of pressure through inter-particle contacts. There are reasons
to believe that constrained powders transmit pressure in much the same way that confined
liquids do.
When a powder is subjected to isostatic pressing, its linear shrinkage is equal in
all directions independent of the particle shape and size distribution. The obtained
compact has uniform density throughout its body. During the diecompaction of a powder,
the pressure-dependence of densification is close to that obtained during isostatic
pressing. In the absence of pressure losses, which are mainly due to die-wall friction,
one can expect the pressure-density relation to be independent of the mode of compaction.
Even under conditions of unidirectional compaction, the pores retain their shape to
a large extent [19]. During sintering, the powder compact undergoes isotropic shrinkage while retaining
of its external shape as if it were subjected to isostatic pressing.
Thus, whether pressure is applied externally (compaction) or derived from within (sintering),
it can be assumed that a hydrostatic stress state prevails during powder consolidation
that is independent of the particle shape and size distribution.
As applied to the compaction of a powder consisting of spherical particles of unequal
sizes, Equation (1) reads as follows:
where Pa is the compaction pressure; Pp is the pressure acting over the surface of a particle of radius a ; and lp is the total load transferred to the particle [20-23]. If the particles are of the same size, each one of them must carry equal load independent
of the number of nearest neighbors (coordination number):
In contrast to liquids, metal powders keep a permanent record of their response to
the applied pressure because of their ability to deform plastically. During compaction
of a metal powder, each particle attempts to indent its neighbors and, as a result,
the inter-particle contacts are converted into contact flats. At a given pressure,
the extent of flattening (total contact area) of each particle depends on the inherent
resistance of the material of the particles to indentation (hardness) and the total
load transferred to it. Thus, measuring the total contact area of the individual particles
of a compacted metal powder makes physical verification of equation (1) and (2) possible.
In this study, FEM simulation were conducted to examine the pressure transmission
(force transmission) during the compaction process. The nickel particles and cylinder
die were modelled by using the ABAQUS 6.13 tool. The simulation results were verified
by comparison experimental results. The relationships between particle arrangement
angle which affect the contact area, compaction speed, and contact force were obtained.
2. SIMULATION CONDITIONS
Figure 2 shows the initial configuration of the simulation. A total of 9 particles were modeled
and located in the cylinder shape die. Vertical speed was applied at the top of the
die, while the bottom surface is stationary and rigid. Because it does not move and
it is not deformed by the upper punch moving.
The particles were Nickel powder provided by POSCO. Experiments to investigate mechanical
properties of elasticity and plasticity were also carried out and the results are
summarized in Table 1.
Also, the boundary conditions were determined as shown in Table 2, to examine the force-displacement relationship and force-contact surface in relation
to compaction speed and different particle arrangement angles.
The compaction speed was set to 50, 100, 200 mm/s and the particle arrangement angles
were chosen to be 0, 22.5, 45 degrees. Figure 3 demonstrates the modeling of the particle arrangement angles.
In this simulation, the equation of equilibrium was applied as below:
where, σ is a stress, fi is a body force. Also, stress-strain relationship in elastic (Hooke’s law) was determined
as below:
The plastic strain rate tensor of the metal powder is defined as follows:
where, Φ and
λ
˙
are a yield function and the scalar value of the powder material, respectively. The
rate of change of the relative density can be written as follows from the mass invariant
relation.
Finally, the plastic yielding behavior of the powder can be defined by the Shima-Oyane
equation as follows [24].
where, q and p are the effective stress and hydrostatic pressure, respectively.
4. RESULTS AND DISCUSSION
4.1. Simulation verification
A compaction experiment was conducted to verify the simulation. The compaction speed
was set to 100mm/s with pressure from 20 MPa to 180 MPa increasing in 20 MPa steps
(a total 9 data points were obtained). Initially, the particle arrangement degree
was 0° to represent the ideal case. Figure 4 shows the initial/final configuration of the simulation. ABAQUS 6.13 was used in
this research with an 8-noded solid element, C3D8, with distortion, hourglass mode
control. Reduced integration was applied to decrease the computational cost. A total
of 102,536 elements were generated with an average aspect ratio of 1.05.
As shown in Fig 4, the results confirmed that the contact length between the particles or die gradually
increased during compaction. A ductile sphere subjected to free-compression under
a load lf undergoes deformation at the points of contact with the flat surfaces forming two
circular contacts of equal radius xf (refer to Fig 5). The mean contact pressure pm is given by
where
S
c
=
2
π
x
f
2
is the total contact area of the sphere. Note that the mean contact pressure does
not depend on the sphere size. Neglecting the deformation-induced increase in the
surface area of the sphere, Eq. (8) can be rewritten as
where a is the radius of the sphere; sr=sc⁄so is the dimensionless relative contact area of the sphere; and so=4πa2 is its original surface area. Equation (9) suggests that the pressure is applied over the surface of the freely compressed sphere
pf=lf⁄2πa2. Evidently, this pressure-term is independent of the material of the sphere.
Figure 6 shows the plot of relative density versus pressure for comparison. Comparing the
experimental and simulation results (Fig 6), the maximum error was within 2.5%, and the average error within 1.2%. Therefore,
the simulation was verified.
4.2. Force/displacement and contact surface area relationship with different compaction
speed
In order to study the effect of compaction speed, simulations were conducted of the
force-displacement and force-contact surface area with 0°, 22.5°, and 45° particle
arrangement angles. Figure 7 shows the force-displacement relationship with constant angles and different compaction
speed.
In 0° and 22.5° case, the force increases almost linearly at 0.11 mm displacement.
After that point (0.11 mm), the slope changes abruptly. This is caused by a change
in contact area. In other words, the particles at the upper part and lower part were
not in contact with each other until the displacement reached 0.11 mm. Therefore,
the force increased rapidly. However, as the velocity changed, there was no significant
difference observed. Even if the particles contact each other, the rate of increase
in contact surface is not sufficiently large [9].
According to the force-contact surface area plot for all cases, force increases linearly
at the transition point (interparticle contact point). After that, force also increases
but the gradient changes. Physically, the gradient (slope) is the pressure (force
/ area) [25-27]. It can be seen that the force is transmitted while the pressure is kept constant,
and the compaction process proceeds with constant pressure even after the particles
come into contact with each other [28,29]. Also, as the particle rearrangement angle increases, the contact points between
the particles gradually move away from the center. Therefore, the point of contact
is gradually delayed (the details are shown in Fig 9).
4.3. Force/displacement and contact surface area relationship with different particle
arrangement angles
When the particle arrangement angle is changed, it can be seen that a difference in
force is caused by the difference in contact area. Until the particles come into contact
with each other, force increases with almost the same slope, but the increase in force
after contact is different. Specifically, even though the force transmission among
the particles is different with constant pressure applied, pressure was constantly
transmitted considering the contact surface area and force. However, the point of
force increase was uniqued for each case due to the particle contact angle, as shown
in Fig 9 (a) and (b)
As shown in Fig 9, when the particle arrangement angle was 22.5 degrees, the powder above contacted
one powder located below. However, when the powder was placed at 45 degrees, the powder
only contacted 2 particles. This indicates, that when the powder is arranged at 45
degrees, the repulsive force between the powder after contact increases, the pressure
increases and the force is transmitted compared to the case of 22.5 degrees. In addition,
as shown in Fig 9 (a), an asymmetry in the contact between the powders appears when they are arranged at
22.5 degrees. This can cause problems, since the powder will be deformed unevenly
during the powder compaction process and finally the density distribution will not
be uniform.
5. CONCLUSION
An analysis of the pressure transmission phenomenon based on compaction speed and
powder arrangement angle was carried out during the powder compaction process. First,
to verify the developed analysis program, compressibility curve were derived through
compression process experiments and analyzed under the same condition. As a result,
it showed an error of 2.5% or less. The force-displacement, and forcecontact area
were derived with different compaction speeds. As the compaction speed increased with
the same particle arrangement angle, there was no significant difference in the change
of force. However, it was confirmed that the point at which the force increases in
relation to the particle arrangement angle changes at the same compaction speed. The
gradient of all force-contact surface areas remained constant before and after contact.
Physically, gradient means pressure in the force-area curve, and the pressure transmission
between particles remains constant as the compaction process progresses. In other
words, the force increases, but the contact area also increases accordingly, so it
can be seen that the pressure stays constant. The conclusions of this paper are summarized
below:
A compaction simulation was conducted at the particle level by using FEM and it determined
the relationship between force-displacement / contact surface area.
It was demonstrated that compaction speed, stacking, and the changes occurring in
particle coordination during the course of consolidation had no effect on the uniformity
of pressure distribution in the powders, and consequently, on the progress of consolidation.
Acknowledgements
This work was supported by the Industrial Technology Innovation Program(No. 10048358)
funded by the Ministry Of Trade, Industry & Energy(MI, Korea)
REFERENCES
Jeong J. H., Park S. J., Comp. Mater. Sci,100, 21 (2015)

Kwon Y. S., Lee H. T., Kim K. T., ASME J. Eng. Mater. Tech,199, 366 (1997)

Kwon Y. S., Chung S. H., Sanderow H. I., Kim K. T., German R. M., Proc. PM2TEC, (2003)

Khoei A. R., Shamloo A., Azami A. R., Int. J. Solids Struct,43, 5421 (2006)

Khoei A. R., Azami A. R., Azizil S., J. Mater. Proc. Technol,185, 166 (2007)

Khoei A. R., Anahid M., Shahim K., J. Mater. Proc. Technol,185, 397 (2007)

Rossi R., Alves M. K., Al-Qureshi H. A., J. Mater. Proc. Technol,182, 286 (2007)

Lee S. C., Kim K. T., Mater. Sci. Eng. A,445, 163 (2007)

Wu C. Y., Hancock B. C., Mills A., Bentham A. C., Best S. M., Elliott J. A., Powder
Technol,181, 121 (2008)

Lee S. C., Kim K. T., Int. J. Mech. Sci,44, 1295 (2002)

Kwon Y. S., Chung S. H., Kim K. T., German R. M., Sanderow H. I., Adv. Powder Metall.
Part. Mater,4, 4 (2003)

Lee M. C., Powder Metall,51, 1 (2008)

Briscoe B. J., Rough S. L., Colloid. Surface. A,89, 103 (1996)

Aydin I., Briscoe B. J., Ozkan N., MRS Bulletin,22, 45 (1997)

PM MODNET Research Group, Powder Metall,45, 335 (2002)

PM MODNET Computer Modelling Group, Powder Metall,42, 301 (1999)

Coube O., Riedel H., Powder Metall,43, 123 (2000)

Smith L. N., Powder Metall,45, 294 (2002)

Fischmeister H. F., Arzt E., Olsson L. R., Powder Metall,21, 179 (1973)

Kumar J. V., Sintered Metal-Ceramic Composites, (1984)

Kumar J. V., Technical Report No. MDN 10 1987, Mishra Dhatu Nigam Limited, (1987)

Kumar J. V., Solid State Phenom,8, 125 (1990)

Kumar J. V., Kumar C. S., Adv. Powder Metall. Part. Mater,2, 5 (1995)

Shima S., Sakamoto Y., Kotera H., Int. J. Mech. Sci,44, 1603 (2002)

Atre S. V., Park S. J., Zauner R., German R. M., Powder Metall,50, 1 (2007)

Chung S. H., Kwon Y. S., Park S. J., German R. M., Finite Element in Analysis and
Design,45, 11 (2009)

Jaeckel L. A., The Annals of Mathematical Statistics, (1972)

Lee J.-H., Park S. J., Moon J., Kang J.-Y., Park J. Y., Lee T.-H., Cho K. M., Korean
J. Met. Mater,55, 363 (2017)

Park J. Y., Park S. J., Lee J. H., Moon J., Lee T. H., Jeong K. J., Shin J. H., Korean
J. Met. Mater,55, 12 (2017)

Rubio G., Agrait N., Vieira S., Phys. Rev. Lett,76, 2302 (1996)

Dubois S., Piraux L., Gerge J. M., Ounadjela K., Duvail J. L., Fert A., Phys. Rev.
B,60, 477 (1999)

Horstemeyer M. F., Bammann D. J., Int. J. Plast,26, 1310 (2010)

Lee D. N., Kim H. S., Powder Metall,35, 275 (1992)

Figures and Tables
Fig. 1.
Overall procedure of PM process
Fig. 2.
Initial geometry of the compaction simulation
Fig. 3.
Particle arrangement angles
Fig. 4.
Compaction simulation at (a) initial configuration, (b) half time configuration, and
(c) final configuration
Fig. 5.
compression of a ductile sphere
Fig. 6.
Pressure-relative density plot for comparing experimental and simulation results
Fig. 7.
Force-displacement and force-contact surface area plots with (a) 0°, (b) 22.5°, and
(c) 45°
Fig. 9.
Force-displacement and force-contact surface area plots with (a) 50 mm/s, (b) 100
mm/s, and (c) 200 mm/s
Table 1.
Material properties of slab heater cover parts
Elasticity
|
Plasticity
|
|
Density (g/cm3)
|
Young’s modulus (GPa)
|
Poisson’s ratio
|
Yield stress (MPa)
|
Plastic strain
|
|
7.85
|
200
|
0.31
|
185
|
0
|
|
290
|
0.05
|
|
345
|
0.1
|
|
388
|
0.15
|
|
413
|
0.2
|
|
434
|
0.25
|
Table 2.
|
|
Boundary condition
|
|
Bottom
|
Fix
|
|
Particle arrangement angles
|
0°, 22.5°, 45°
|
|
Compaction speed
|
50 mm/s
|
100 mm/s
|
200 mm/s
|