OS-V: 0840 Viscoplastic Anand Material Model

To analyze the rate-dependent deformation of metals at high temperatures, the Anand model was suggested. The model applies to a wide range of viscoplastic problems, including the influence of the strain rate and temperature.

Model Files

Benchmark Model

The specimen used for FEM simulation is shown below, along with its dimensions.

Figure 1. Specimen FE Model

On one side of the specimen, all the nodes are constrained to three degrees of freedom (123).

Additionally, on the other side, RBE2 is applied on the entire phase and the center of the RBE2 element is constrained to three degrees of freedom (123), along with constraint 1 DOF (1) using the SPCD entry.

Units: mm, s, Mg, N, MPa

Material

The MAT1 Bulk Data Entry is used to input material data for the model.

The Poisson’s ratio ‘υ’ is defined on the NU field. Similarly, Young's Modulus and Density are defined in the E and RHO fields, respectively.

The Anand material model is assigned through MATVP. The CTYPE field on the MATVP entry is used to identify the creep material model type. Specifying a MAT1 and a MATVP Bulk Data Entry with the same MID allows modeling creep material.

The TABLED1 Bulk Data Entry is used to apply enforced displacement in TLOAD1 and enforced temperature in TLOAD1.

Figure 2. Displacement Curve Input for 375 seconds

Figure 3. Temperature Curve Input for 80 degrees Celsius

Temperature load is applied to the component.

The DLOAD Bulk Data Entry is used to give a dynamic loading condition for the load sets defined on the TLOAD1 entry.

With the obtained final parameter values of the Anand material model, the simulations were performed for a graded tensile test at 80 degrees Celsius.

Results

Figure 4. SPCF Forces versus Time at the center SPCD node of the RBE2

Figure 5. Average True ZZ Stresses at the cross-section of the specimen

Figure 6. Force versus Time

Since Force, as an output, is not supported for solid elements, average True zz stresses at the cross-section of the specimen, along with the dynamically changing area at each timestep are used to manually calculate the force.