ABAQUS Benchmarks Manual  


1.1.9 Necking of a round tensile bar

Products: ABAQUS/Standard  ABAQUS/Explicit  

This example illustrates necking and softening of a round tensile bar. This problem has been studied by Needleman (1972), Tvergaard and Needleman (1984), Needleman and Tvergaard (1985), and Aravas (1987). The material is assumed to be a rate-independent metal in which triaxial tension stress can allow voids to nucleate and grow. The example illustrates the use of the *POROUS METAL PLASTICITY model and the *VOID NUCLEATION suboption. In ABAQUS/Explicit the *POROUS FAILURE CRITERIA suboption is used to model failure of the material after a critical void volume fraction is reached.

Problem description

We consider a long specimen with a circular cross-section. The specimen has an initial length of 2 and a radius of , with = 4. is assumed to be equal to 1 unit. Only a quarter of the specimen needs to be analyzed because of the symmetry about the and axes. Figure 1.1.9–1 shows the mesh used in the analysis. Both the geometry and the deformation are assumed to be axisymmetric. Axisymmetric elements are used, and the mesh is refined near the center of the specimen because of the expected softening and intense deformation in that region. An initial geometric imperfection is used to induce necking in the specimen analyzed with ABAQUS/Standard. In ABAQUS/Explicit the imperfection is not needed because stress wave effects induce necking at the center of the bar.

Material

The material properties used in the computation are:

Young's modulus, E:300
Poisson's ratio, :0.3
Porous material parameters: = 1.5, = 1.0, and = 2.25
Initial relative density:1.0 ( = 0.0)
Void nucleation parameters: = 0.3, = 0.1, and = 0.04
Porous failure criteria: = 0.6, = 0.59 (for ABAQUS/Explicit only)
The work hardening behavior (yield stress, , versus equivalent plastic strain, ) given under the *PLASTIC option is of the form

where = 1 is the initial yield stress, = 0.1 is the hardening parameter, and is the elastic shear modulus. Necking is expected to start when the yield stress approaches the work hardening rate, which occurs at a strain of about 10 to 12%. Hence, the work hardening behavior is described more accurately for 0.08 0.3 than for the rest of the curve.

The parameters , , and were introduced by Tvergaard (1981) to make the predictions of the Gurson model agree with numerical studies of an elastic-plastic medium containing a periodic array of voids. The parameter values used in this analysis are those chosen by Tvergaard.

The void nucleation parameters used in the material description are the same as those given by Tvergaard and Needleman (1984) and Needleman and Tvergaard (1985). These parameters describe the normal distribution of the nucleation strain (see “Porous metal plasticity,” Section 11.2.9 of the ABAQUS Analysis User's Manual). The area under the normal distribution curve represents the total volume fraction of the nucleated voids and is approximately equal to . With the normal distribution, the amount of voids nucleated between 0.2 and 0.4 is about 68% of .

Boundary conditions and loading

The kinematic boundary conditions are symmetry about (all nodes along have prescribed) and symmetry about (all nodes along have prescribed). All the nodes on the top of the specimen along 4.0, in the node set TOPSIDE, are pulled in the z-direction while being held fixed in the r-direction. In the ABAQUS/Explicit analysis the nodes in node set TOPSIDE are pulled with a prescribed velocity that increases linearly from 0 to 30 at 0.025 s and then decreases linearly from 30 to 0 at 0.05 s; in the ABAQUS/Standard analysis the displacement is applied directly to obtain the deformations desired in the two analysis steps described below.

In the ABAQUS/Standard analysis the accuracy of the implicit integration of the void nucleation and growth equation is controlled by prescribing a maximum allowable time increment in the automatic time incrementation scheme.

Results and discussion

The example problem focuses on the neck development, which is a precursor to failure in the form of “cup-cone fracture.” The formation of the neck results in a triaxial state of stress at the center of the specimen, which accelerates the growth of the nucleated voids. A detailed analysis of the cup-cone fracture can be found in Tvergaard and Needleman (1984), which predicts that void nucleation is followed by the formation of a planar crack at the center of the neck as a result of the coalescence of voids. The planar crack propagates along a zig-zag path closer to the traction-free surface, eventually leading to the formation of the well-known cup-cone fracture.

ABAQUS/Standard results

The calculations in the first step are terminated at an overall nominal strain 19%, thereby making it possible to compare the results with those of Aravas (1987). In the second step the calculations are carried on further to an overall nominal strain of 19.75%.

The results of the analysis are illustrated in Figure 1.1.9–2 to Figure 1.1.9–6. Figure 1.1.9–2 shows the computed force as a function of the overall nominal strain. The maximum load is reached at an overall nominal strain of about 10.2%. The nominal stress–nominal strain curve, as well as the contour plots of void volume fraction (Figure 1.1.9–3) and hydrostatic pressure (Figure 1.1.9–4) at an overall nominal strain of 19%, match well with the results obtained by Aravas (1987).

Figure 1.1.9–5 and Figure 1.1.9–6 show the contour plots of void volume fraction and hydrostatic pressure at an overall nominal strain of 19.75%; a comparison with Figure 1.1.9–3 and Figure 1.1.9–4 reveals a significant growth of voids and a corresponding decrease of the hydrostatic tension in the neck region, indicative of the material softening that has taken place.

ABAQUS/Explicit results

The results of the analysis are illustrated in Figure 1.1.9–2 and Figure 1.1.9–7 through Figure 1.1.9–10. Figure 1.1.9–2 shows the computed nominal stress as a function of the nominal strain. The maximum load is reached at a nominal strain of about 9%, after which the specimen softens due to coalescence of voids and eventually fractures across the neck region. Due to the relatively high speed of the loading in the ABAQUS/Explicit analysis, the void growth and coalescence and the failure propagation are coupled with dynamic effects. These dynamic effects are the source of the small differences observed in the results obtained with ABAQUS/Explicit and ABAQUS/Standard. Figure 1.1.9–7 and Figure 1.1.9–8 show total void volume fraction and pressure stress contours at 22.88% ( 0.0317 s). Figure 1.1.9–9 shows the broken tensile specimen (at 0.05 s), where only the elements whose void volume fraction is still below the ultimate failure ratio are shown. The deformed mesh is shown next to the initial mesh. Figure 1.1.9–10 shows contours of pressure in the broken bar. As was mentioned earlier, the tensile bar typically fails in a “cup-cone fracture;” because a symmetric solution was assumed in this model, a proper cup-cone fracture cannot develop in this case.

Input files

ABAQUS/Standard input files

ABAQUS/Explicit input files

neck.inp

ABAQUS/Explicit analysis.

neck_ale.inp

Model using the *ADAPTIVE MESH option.

neck_ef1.inp

External file referenced in the adaptive mesh input file.

neck_ef2.inp

External file referenced in the adaptive mesh input file.

References

Figures

Figure 1.1.9–1 Geometry and mesh for the round tensile bar.

Figure 1.1.9–2 Overall nominal stress, , vs. overall nominal strain, .

Figure 1.1.9–3 Void volume fraction at 19% (ABAQUS/Standard).

Figure 1.1.9–4 Hydrostatic pressure at 19% (ABAQUS/Standard).

Figure 1.1.9–5 Void volume fraction at 19.75% (ABAQUS/Standard).

Figure 1.1.9–6 Hydrostatic pressure at 19.75% (ABAQUS/Standard).

Figure 1.1.9–7 Total void volume fraction at 22.88% (ABAQUS/Explicit).

Figure 1.1.9–8 Hydrostatic pressure at 22.88% (ABAQUS/Explicit).

Figure 1.1.9–9 Final broken bar and its initial configuration (ABAQUS/Explicit).

Figure 1.1.9–10 Hydrostatic pressure in the broken bar (ABAQUS/Explicit).