OS-V: 0840 Viscoplastic Anand Material Model

Knowledge of the material model and its parameters is essential to predict the behavior of the printed components for FEM Analysis. Anand is one of the material models used to study viscoplastic materials.

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

Before you begin, copy the file(s) used in this problem to your working directory.

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


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.


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.

1 Parameters Identification of the Anand Material Model for 3D Printed Structures’ Department of Applied Mechanics, Faculty of Mechanical Engineering, VŠB—Technical University of Ostrava,708 00 Ostrava, Czech Republic; martin.fusek@vsb.cz (M.F.); jaroslav.rojicek@vsb.cz (J.R.); frantisek.fojtik@vsb.cz (F.F.)