Nonlinear finite elements/Homework 11/Solutions/Problem 2

Problem 2: Billet upset forging

Given:

• Compression of cylindrical billet with rigid dies.
• No-slip contact between billet and dies.
• $r_0$ = 10 mm; $l_0$ = 30 mm.
• $\mu$ = 384.6 MPa; $K$ = 833.3 MPa.
• Yield stress given by

${\sigma }_y={\sigma }_0+H{\epsilon }_p.$

• ${\sigma }_0$ = 1 MPa; $H$ = 3 MPa.
• Compressive force of 1 KN applied to die.

Solution

Part 1

This problem is solved using ANSYS 8.0. Plane182 (2-D 4-node axisymmetric element) is used to model a quarter of the billet and the die. Contac48 (2-D contact element) is used for the node to surface contact. The final shape of the billet is shown in Fig 4.

File:NFE HW11 p2 def.png Figure 4. Deformed shape of the billet at the final load step.

Part 2

The die moves less in the ANSYS simulation than it does in the plot shown in Simo and Hughes (p. 325, Fig. 9.8 (b)). In Simo and Hughes, the die stroke is around 20.5 mm. In our simulation, the die stroke is around 19 mm.

The top of the billet shows an undulation that indicates that contact is not being done correctly. The billet material appears to be moving through the die during deformation. Simo and Hughes do not show such an effect.

Part 3

A plot of the die force vs. the displacement at the center of the billet is shown in Fig 5. The die stroke is (30 - $u$) mm where $u$ is the displacement shown in the plot.

File:NFE HW11 p2 fd.png Figure 5. Billet force-displacement plot.

Part 4

See Fig 9.8 (c) in Simo and Hughes. Our simulation shows a smaller displacement for the same load when compared with Simo and Hughes.

Part 5

The smaller displacement observed in our simulation may be an indication of volumetric locking when fully integrated PLANE182 elements are used. Underintegrated elements or assumed strain elements may alternatively be used to confirm that possibility.