Case01 von Mises

 

Case 01 is based on the von-Mises elasto-plastic model with properties defined as:

 

Parameter

Value

Young's modulus, E (MPa)

200(103)

Poisson's ratio, ν

0.3

Yield stress, σy (MPa)

30(103)

Hardening slope, H (MPa)

30(103)

von Mises Material Properties

 

 

Following are descriptions of some key data and the element types investigated for the implicit and explicit analyses.

Material_data

 

Material_data defines the material properties of the model.

 

Data File

 

 

* Material_data               NUM=1

! ---------------------------------

 Material_name            "Billet"

 Grain_density                7833                

 Grain_stiffness             30000

 Porosity                     0.02

 Elastic_properties          IDM=2              

  /Young's Modulus/          200E3

 /Poisson's ratio/           0.30

 Plastic_material_type                    7

 Plastic_properties                IDM=1

 /Yield Strength/              30E3

 Hardening_type                            7

 Hardening_properties                IDM=2         JDM=2

  /Effective plastic strain/        0      1.00

  /Yield strength/                        30E3   60E3

 

1Material name is defined as "Billet".

2Grain properties (density = 7833 kg/m3 and stiffness) are defined.

3Elastic properties (Young's modulus = 200(103)MPa, Poisson's ratio = 0.30) are defined.

4Plastic material type 7, corresponding to von-Mises elastoplastic model is defined.

5The yield strength of von-Mises model is defined as 30(103) MPa, with hardening slope of 30(103) MPa.

 

 

Group_data

 

Group_data sets the group name, element type, material name, porous flow type and the associated volume entity.

 

Data File

 

 

* Group_data               NUM=1  

! -------------------------------

Group_name            "Volume1"

 Element_type            "HEX8M"

 Material_name          "Billet"

 Porous_flow_type             1

 Volumes  IDM=1

   1

1Group name and material name are defined.

2Porous flow type 1 (corresponding to dry porous media) is defined.

3Volume set 1 is registered to be associated with the current group data structure.

4Element types that are used in this example include

 

(a) Implicit analysis

Element Type

Element Descriptions

HEX8_Bbar

3D 8-noded hexahedral Bbar-averaged element

TET10

3D 10-noded tetrahedral element

HEX8M

3D 8-noded hexahedral mixed element

HEX20M

3D 20-noded hexahedral mixed element

TET4M

3D 4-noded tetrahedral mixed element

TET10M

3D 10-noded tetrahedral mixed element

 

(b) Explicit analysis

Element Type

Element Descriptions

HEX8

3D 8-noded hexahedral element with hour-glass control

 

The hourglass control is called by adding Mech_integration_type    "Belytschko_pressure" under Group_data data structure.

 

 

 

Results

 

The following figures show the development of axial, volumetric and radial strains in terms of averaged deviatoric stress. The latter is computed as the average of the whole billet. As the compression load increases, axial strain increases in negative value while the radial strain increases in the positive value due to Poisson’s effect. When the stress state reaches the von-Mises yield surface, both strain parameters increases sub-linearly. Volumetric strain is computed as εvol = εaxial + 2εradial.

 

Val_004d_03

Billet compression as the applied displacement increases

 

 

We observe volumetric compression during elastic stage because the effect of axial strain is more dominant. In plastic stage, however, the contribution of radial strain to the volumetric strain has become more dominant, and thus we observe volumetric expansion. This trend would lead to plastic collapse if the compression load continues further.

 

Val_004d_04

Development of axial, volumetric and radial strains in terms of averaged deviatoric stress for standard elements

 

 

 

Val_004d_05

Development of axial, volumetric and radial strains in terms of averaged deviatoric stress for elements with mixed formulation

 

 

The figures above are separated by the standard element group and the mixed element group mainly to ease visualisation, lest legend be overloaded by the labels. It should be mentioned that the plot of strain parameters differ slightly between the standard and mixed element groups. This difference is also highlighted in the figures where the averaged mean stress and the averaged deviatoric stress are plotted against radial displacement. Compared to the standard elements, mixed elements exhibit smaller increase in stress for each increment of radial displacement. It is expected that elements associated with mixed formulation can alleviate the averaged mean stress better than the displacement-based elements.

 

Val_004d_06

Averaged mean stress plotted against radial displacement

 

 

 

Val_004d_07

Averaged deviatoric stress plotted against radial displacement

 

 

The plastic strain and effective mean stress evolution as shown in the figure below is driven by the concurrency of both axial and radial strain development. The difference between the axial and radial stress results in distortional (shape-changing) deformation within the billet. This is reflected by the plastic strain. As the circumferential surface of the billet is progressively displaced sideways, the effective mean stress is observed to increase and also become gradually concentrated towards the lateral surface, signifying the imminence of plastic collapse.

 

Val_004d_08

Development of plastic strain and effective mean stress on HEX8M model in clipped view

 

 

 

Val_004d_09

Comparison of plastic strain and effective mean stress distribution at final time step for different element types in clipped view

 

 

References

 

[1] Lippmann,  H. (1979). Metal Forming Plasticity. Springer-Verlag, Berlin.

[2] Cheng., J.H., Kikuchi, N. (1984). An analysis of metal forming processes using large deformation elastic-plastic formulations

[3] https://abaqus-docs.mit.edu/2017/English/SIMACAEEXARefMap/simaexa-c-rezonebillet.htm