VM318

VM318
Simulation of Necking of a Circular Bar

Overview

Reference:

Simo, J. C. (1992). Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering, (99), 61-112.

Analysis Type(s): Static (ANTYPE,0)
Element Type(s):

2-D Structural Elements (PLANE182)

3-D Structural Elements (SOLID185)

Input Listing: vm318.dat

Test Case

A bar with a circular cross section is under uniaxial tension. The bar is geometrically imperfect: The radius reduces linearly from the bar end R to the bar middle which equals 0.982*R. The bar, however, is symmetrical along the axis. When the bar elongates under the uniaxial tension, the cross section in the middle starts to shrink. Displacement at the middle section is obtained and compared to the reference results.

Material Properties Geometric Properties Loading

Neo-Hookean Model:

μ = 80.1938 MPa
d = 2/164.21 (1/MPa)

Nonlinear isotropic hardening (VOCE):

σ0 = 0.45 MPa
R0= 0.12924 MPa
R = 0.265 MPa
b = 16.93

Bar is circular with radius linearly reduced from the bar end R to the bar middle 0.982*R.

R = 6.413 mm

The length L of the bar:

L = 53.334 mm
Bar elongates 7 mm at one end.

Analysis Assumptions and Modeling Notes

The problem is solved in two ways:

  1. A 2D model using axisymmetric PLANE182

  2. A 3D model using SOLID185

The material is defined by TB,HYPER,,,,NEO and TB, NLISO,,,,VOCE. Because the bar is symmetrical along the axis, a half-bar is simulated. NLGEOM,ON must be set in the solution setting. After the bar elongates 7 mm, the necking of the bar is plotted by PLNSOL in /POST1 (see Figure 580: Final Deformation of the Circular Bar) and obtained by *GET in /POST1 and NSOL in /POST26 (see Figure 581: Displacement at the Middle Section Versus Total Elongation).

Figure 580: Final Deformation of the Circular Bar

Final Deformation of the Circular Bar

Figure 581: Displacement at the Middle Section Versus Total Elongation

Displacement at the Middle Section Versus Total Elongation

Stress Results Comparison

 TargetMechanical APDLRatio
U_X (2D mode)-3.740-3.7360.9989
U_X (3D mode)-3.740-3.8011.0163