Unpublished,

Shear Locking and Hourglassing in MSC Nastran, ABAQUS, and ANSYS

.
(2006)

Abstract

A solid beam and a composite beam were used to compare how MSC Nastran, ABAQUS, and ANSYS handled the numerical difficulties of shear locking and hourglassing. Their tip displacements and first modes were computed, normalized, and listed in multiple tables under various situations. It was found that fully integrated first order solid elements in these three finite element codes exhibited similar shear locking. It is thus recommended that one should avoid using this type of element in bending applications and modal analysis. There was, however, no such shear locking with fully integrated second order solid elements. Reduced integration first order solid elements in ABAQUS and ANSYS suffered from hourglassing when a mesh was coarse. If there was only one layer of elements, the reported first mode of the beam examples from ABAQUS and ANSYS was excessively smaller than the converged solutions due to hourglassing. At least four layers of elements should, therefore, be used in ABAQUS and ANSYS. MSC Nastran outperformed ABAQUS and ANSYS by virtually eliminating the annoying hourglassing of reduced integration first order 3D solid elements because it employed bubble functions to control the propagation of non-physical zero-energy modes. Even if there was only one layer of such elements, MSC Nastran could still manage to produce reasonably accurate results. This is very convenient because it is much less prone toA solid beam and a composite beam were used to compare how MSC Nastran, ABAQUS, and ANSYS handled the numerical difficulties of shear locking and hourglassing. Their tip displacements and first modes were computed, normalized, and listed in multiple tables under various situations. It was found that fully integrated first order solid elements in these three finite element codes exhibited similar shear locking. It is thus recommended that one should avoid using this type of element in bending applications and modal analysis. There was, however, no such shear locking with fully integrated second order solid elements. Reduced integration first order solid elements in ABAQUS and ANSYS suffered from hourglassing when a mesh was coarse. If there was only one layer of elements, the reported first mode of the beam examples from ABAQUS and ANSYS was excessively smaller than the converged solutions due to hourglassing. At least four layers of elements should, therefore, be used in ABAQUS and ANSYS. MSC Nastran outperformed ABAQUS and ANSYS by virtually eliminating the annoying hourglassing of reduced integration first order 3D solid elements because it employed bubble functions to control the propagation of non-physical zero-energy modes. Even if there was only one layer of such elements, MSC Nastran could still manage to produce reasonably accurate results. This is very convenient because it is much less prone to errors when using reduced integration first order 3D solid elements in MSC Nastran.

Tags

Users

  • @gdmcbain

Comments and Reviews