Products: ABAQUS/Standard ABAQUS/Explicit
This example illustrates coupled temperature-displacement analysis in a metal forming application. The case studied is an extension of the standard test case that is defined in Lippmann (1979); thus, some verification of the results is available by comparison with the numerical results presented in that reference. The example is that of a small, circular billet of metal that is reduced in length by 60%. Here the problem is analyzed as a viscoplastic case, including heating of the billet by plastic work. Such analysis is often important in manufacturing processes, especially when significant temperature rises degrade the material. The problem is also analyzed in ABAQUS/Standard using a porous metal material model. The same problem is used in Upsetting of a cylindrical billet: quasi-static analysis with mesh-to-mesh solution mapping (ABAQUS/Standard) and adaptive meshing (ABAQUS/Explicit), Section 1.3.1, to illustrate mesh rezoning in ABAQUS/Standard and adaptive meshing in ABAQUS/Explicit.
The specimen is shown in Figure 1.3.161: a circular billet, 30 mm long, with a radius of 10 mm, compressed between flat, rough, rigid dies. All surfaces of the billet are assumed to be fully insulated: this thermal boundary condition is chosen to maximize the temperature rise.
The finite element model is axisymmetric and includes the top half of the billet only since the middle surface of the billet is a plane of symmetry. In ABAQUS/Standard elements of type CAX8RT, 8-node quadrilaterals with reduced integration that allow for fully coupled temperature-displacement analysis, are used. A regular mesh with six elements in each direction is used, as shown in Figure 1.3.161. In addition, the billet is modeled with CAX4RT elements in a 12 × 12 mesh for both ABAQUS/Standard and ABAQUS/Explicit analyses.
The contact between the top and the lateral exterior surfaces of the billet and the rigid die is modeled with the *CONTACT PAIR option. The billet surface is defined by means of the *SURFACE option. The rigid die is modeled as an analytical rigid surface or as an element-based rigid surface, using the *RIGID BODY option in conjunction with the *SURFACE option. The mechanical interaction between the contact surfaces is assumed to be nonintermittent, rough frictional contact in ABAQUS/Standard. Therefore, two options are used in conjunction with the *SURFACE INTERACTION property option: the *FRICTION, ROUGH option to enforce a no-slip constraint between the two surfaces and the *SURFACE BEHAVIOR, NO SEPARATION option to ensure that separation does not occur once contact has been established. In ABAQUS/Explicit the friction coefficient between the billet and the rigid die is 1.0.
The problem is also solved in ABAQUS/Standard with the first-order fully coupled temperature-displacement CAX4T elements in a 12 × 12 mesh. Similarly, the problem is solved using CAX8RT elements and user subroutines UMAT and UMATHT to illustrate the use of these subroutines.
No mesh convergence studies have been performed, but the comparison with results given in Lippmann (1979) suggests that these meshes provide accuracy similar to the best of those analyses.
The ABAQUS/Explicit simulations are performed both with and without adaptive meshing.
The material definition is basically that given in Lippmann (1979), except that the metal is assumed to be rate dependent. The thermal properties are added, with values that correspond to a typical steel, as well as the data for the porous metal plasticity model. The material properties are then as follows:
Young's modulus: | 200 GPa |
Poisson's ratio: | 0.3 |
Thermal expansion coefficient: | 1.2×105 per °C |
Initial static yield stress: | 700 MPa |
Work hardening rate: | 300 MPa |
Strain rate dependence: | ; /s, |
Specific heat: | 586 J/(kg°C) |
Density: | 7833 kg/m3 |
Conductivity: | 52 J/(m-s-°C) |
Porous material parameters: | |
Initial relative density: | 0.95 ( 0.05) |
Since the problem definition in ABAQUS/Standard assumes that the dies are completely rough, no tangential slipping is allowed wherever the metal contacts the die.
The kinematic boundary conditions are symmetry on the axis (nodes at 0, in node set AXIS, have 0 prescribed) and symmetry about 0 (all nodes at 0, in node set MIDDLE, have 0 prescribed). To avoid overconstraint, the node on the top surface of the billet that lies on the symmetry axis is not part of the node set AXIS: the radial motion of this node is already constrained by a no-slip frictional constraint (see Common difficulties associated with contact modeling in ABAQUS/Standard, Section 29.2.11 of the ABAQUS Analysis User's Manual, and Common difficulties associated with contact modeling using the contact pair algorithm in ABAQUS/Explicit, Section 29.4.6 of the ABAQUS Analysis User's Manual). The rigid body reference node for the rigid surface that defines the die is constrained to have no rotation or -displacement, and its -displacement is prescribed to move 9 mm down the axis at constant velocity. The reaction force at the rigid reference node corresponds to the total force applied by the die.
The thermal boundary conditions are that all external surfaces are insulated (no heat flux allowed). This condition is chosen because it is the most extreme case: it must provide the largest temperature rises possible, since no heat can be removed from the specimen.
DELTMX is the limit on the maximum temperature change allowed to occur in any increment and is one of the controls for the automatic time incrementation scheme in ABAQUS/Standard. It is set to 100°C, which is a large value and indicates that we are not restricting the time increments because of accuracy considerations in integrating the heat transfer equations. In fact, the automatic time incrementation scheme will choose fairly small increments because of the severe nonlinearity present in the problem and the resultant need for several iterations per increment even with a relatively large number of increments. The large value is used for DELTMX to obtain a reasonable solution at low cost.
In ABAQUS/Explicit the automatic time incrementation scheme is used to ensure numerical stability and to advance the solution in time. Mass scaling is used to reduce the computational cost of the analysis.
The AMPLITUDE=RAMP parameter is included because the default amplitude variation for a transient, coupled temperature-displacement analysis is a step function, but here we want the die to move down at a constant velocity.
Two versions of the analysis are run: a slow upsetting, where the upsetting occurs in 100 seconds, and a fast upsetting, where the event takes 0.1 second. Both versions are analyzed with the coupled temperature-displacement procedure. The fast upsetting is also run in ABAQUS/Standard as an adiabatic static stress analysis. The time period values are specified on the data line associated with the *COUPLED TEMPERATURE-DISPLACEMENT procedure, *DYNAMIC TEMPERATURE-DISPLACEMENT procedure, and the *STATIC procedure options. The adiabatic stress analysis is performed in the same time frame as the fast upsetting case. In all cases analyzed with ABAQUS/Standard an initial time increment of 1.5% of the time period is used; that is, 1.5 seconds in the slow case and 0.0015 second in the fast case. This value is chosen because it will result in a nominal axial strain of about 1% per increment, and experience suggests that such increment sizes are generally suitable for cases like this.
The results of the ABAQUS/Standard simulations are discussed first, beginning with the results for the viscoplastic fully dense material. The results of the slow upsetting are illustrated in Figure 1.3.162 to Figure 1.3.164. The results for the fast upsetting coupled temperature-displacement analysis are illustrated in Figure 1.3.165 to Figure 1.3.167; those for the adiabatic static stress analysis are shown in Figure 1.3.168 and Figure 1.3.169. Figure 1.3.162 and Figure 1.3.165 show the configuration that is predicted at 60% upsetting. The configuration for the adiabatic analysis is not shown since it is almost identical to the fast upsetting coupled case. Both the slow and the fast upsetting cases show the folding of the top outside surface of the billet onto the die, as well as the severe straining of the middle of the specimen. The second figure in each series (Figure 1.3.163 for the slow case, Figure 1.3.166 for the fast case, and Figure 1.3.168 for the adiabatic case) shows the equivalent plastic strain in the billet. Peak strains of around 180% occur in the center of the specimen. The third figure in each series (Figure 1.3.164 for the slow case, Figure 1.3.167 for the fast case, and Figure 1.3.169 for the adiabatic case) shows the temperature distributions, which are noticeably different between the slow and fast upsetting cases. In the slow case there is time for the heat to diffuse (the 60% upsetting takes place in 100 sec, on a specimen where a typical length is 10 mm), so the temperature distribution at 100 sec is quite uniform, varying only between 180°C and 185°C through the billet. In contrast, the fast upsetting occurs too quickly for the heat to diffuse. In this case the middle of the top surface of the specimen remains at 0°C at the end of the event, while the center of the specimen heats up to almost 600°C. There is no significant difference in temperatures between the fast coupled case and the adiabatic case. In the outer top section of the billet there are differences that are a result of the severe distortion of the elements in that region and the lack of dissipation of generated heat. The temperature in the rest of the billet compares well. This example illustrates the advantage of an adiabatic analysis, since a good representation of the results is obtained in about 60% of the computer time required for the fully coupled analysis.
The results of the slow and fast upsetting of the billet modeled with the porous metal plasticity model are shown in Figure 1.3.1610 to Figure 1.3.1615. The deformed configuration is identical to that of Figure 1.3.162 and Figure 1.3.165. The extent of growth/closure of the voids in the specimen at the end of the analysis is shown in Figure 1.3.1610 and Figure 1.3.1613. The porous material is almost fully compacted near the center of the billet because of the compressive nature of the stress field in that region; on the other hand, the corner element is folded up and stretched out near the outer top portion of the billet, increasing the void volume fraction to almost 0.1 (or 10%) and indicating that tearing of the material is likely. The equivalent plastic strain is shown in Figure 1.3.1611 (slow upsetting) and Figure 1.3.1614 (fast upsetting) for the porous material; Figure 1.3.1612 and Figure 1.3.1615 show the temperature distribution for the slow and the fast upsetting of the porous metal. The porous metal needs less external work to achieve the same deformation compared to a fully dense metal. Consequently, there is less plastic work being dissipated as heat; hence, the temperature increase is not as much as that of fully dense metal. This effect is more pronounced in the fast upsetting problem, where the specimen heats up to only 510°C, compared to about 600°C for fully dense metal.
Figure 1.3.1616 to Figure 1.3.1618 show predictions of total upsetting force versus displacement of the die. In Figure 1.3.1616 the slow upsetting viscoplastic and porous plasticity results are compared with several elastic-plastic and rigid-plastic results that were collected by Lippmann (1979) and slow viscoplastic results obtained by Taylor (1981). There is general agreement between all the rate-independent results, and these correspond to the slow viscoplastic results of the present example and of those found by Taylor (1981). In Figure 1.3.1617 rate dependence of the yield stress is investigated. The fast viscoplastic and porous plasticity results show significantly higher force values throughout the event than the slow results. This effect can be estimated easily. A nominal strain rate of 6 sec is maintained throughout the event. With the viscoplastic model that is used, this effect increases the yield stress by 68%. This factor is very close to the load amplification factor that appears in Figure 1.3.1617. Figure 1.3.1618 shows that the force versus displacement prediction of the fast viscoplastic adiabatic analysis agrees well with the fully coupled results.
Two cases using an element-based rigid surface to model the die are also considered in ABAQUS/Standard. To define the element-based rigid surface, the elements are assigned to rigid bodies using the *RIGID BODY, ISOTHERMAL=YES option. The results agree very well with the case when the analytical rigid surface is used.
The automatic load incrementation results suggest that overall nominal strain increments of about 2% per increment were obtained, which is slightly better than what was anticipated in the initial time increment suggestion. These values are typical for problems of this class and are useful guidelines for estimating the computational effort required for such cases.
The results obtained with ABAQUS/Explicit compare well with those obtained with ABAQUS/Standard, as illustrated in Figure 1.3.1619, which compares the results obtained with ABAQUS/Explicit (without adaptive meshing) for the total upsetting force versus the displacement of the die against the same results obtained with ABAQUS/Standard. The agreement between the two solutions is excellent. Similar agreement is obtained with the results obtained from the ABAQUS/Explicit simulation using adaptive meshing. The mesh distortion is significantly reduced in this case, as illustrated in Figure 1.3.1620.
Slow upsetting case with 144 CAX4T elements, using the fully dense material.
Fast upsetting case with 144 CAX4T elements, using the fully dense material.
Slow upsetting case with 144 CAX4RT elements, using the fully dense material.
Fast upsetting case with 144 CAX4RT elements, using the fully dense material.
Slow upsetting case with CAX8RT elements, using the fully dense material.
Slow upsetting case with CAX8RT elements, using the fully dense material and an element-based rigid surface for the die.
Fast upsetting case with CAX8RT elements, using the fully dense material.
Slow upsetting case with CAX8RT elements, using the porous material.
Fast upsetting case with CAX8RT elements, using the porous material.
Slow upsetting case with 144 CGAX4T elements, using the fully dense material.
Fast upsetting case with 144 CGAX4T elements, using the fully dense material.
Fast upsetting case with 144 CGAX4T elements, using the fully dense material and an element-based rigid surface for the die.
Fast upsetting case with 144 CGAX4T elements, using the fully dense material and an element-based rigid surface for the die with surface-to-surface contact.
Slow upsetting case with CGAX8RT elements, using the fully dense material.
Fast upsetting case with CGAX8RT elements, using the fully dense material.
Adiabatic static analysis with fully dense material modeled with C3D10M elements.
Adiabatic static analysis with fully dense material modeled with C3D10M elements using surface-to-surface contact.
*POST OUTPUT analysis of cylbillet_c3d10m_adiab_dense.inp.
Adiabatic static analysis with fully dense material modeled with CAX6M elements.
Adiabatic static analysis with fully dense material modeled with CAX8R elements.
*POST OUTPUT analysis using the fully dense material.
Slow upsetting case with fully dense material modeled with CAX4RT elements and without adaptive meshing; kinematic mechanical contact.
Fast upsetting case with fully dense material modeled with CAX4RT elements and without adaptive meshing; kinematic mechanical contact.
Slow upsetting case with fully dense material modeled with CAX4RT elements and with adaptive meshing; kinematic mechanical contact.
Fast upsetting case with fully dense material modeled with CAX4RT elements and with adaptive meshing; kinematic mechanical contact.
Fast upsetting case with fully dense material modeled with CAX4RT elements and without adaptive meshing; penalty mechanical contact.
Lippmann, H., Metal Forming Plasticity, Springer-Verlag, Berlin, 1979.
Taylor, L. M., A Finite Element Analysis for Large Deformation Metal Forming Problems Involving Contact and Friction, Ph.D. Thesis, U. of Texas at Austin, 1981.
Figure 1.3.162 Deformed configuration at 60% upsetting: slow case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.163 Plastic strain at 60% upsetting: slow case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.164 Temperature at 60% upsetting: slow case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.165 Deformed configuration at 60% upsetting: fast case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.166 Plastic strain at 60% upsetting: fast case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.167 Temperature at 60% upsetting: fast case, coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.168 Plastic strain at 60% upsetting: fast case, adiabatic stress analysis, ABAQUS/Standard.
Figure 1.3.169 Temperature at 60% upsetting: fast case, adiabatic stress analysis, ABAQUS/Standard.
Figure 1.3.1610 Void volume fraction at 60% upsetting: porous material, slow coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.1611 Plastic strain at 60% upsetting: porous material, slow coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.1612 Temperature at 60% upsetting: porous material, slow coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.1613 Void volume fraction at 60% upsetting: porous material, fast coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.1614 Plastic strain at 60% upsetting: porous material, fast coupled temperature-displacement analysis, ABAQUS/Standard.
Figure 1.3.1615 Temperature at 60% upsetting: porous material, fast coupled temperature-displacement analysis, ABAQUS/Standard.