FEM logo

Generate nice meshes at areas of interest and prevent shear locking.

When modeling a model using shell (plate) elements, the model should comply with the following guidelines to increase the reliability of your results:

The code checker works with the element stresses. Try to reduce the use of triangular elements as much as possible. In The NX Nastran manual the following direction is given:

“You should use the CQUAD4 element when the surfaces you are meshing are reasonably flat and the geometry is nearly rectangular. For these conditions, the quadrilateral elements eliminate the modeling bias associated with the use of triangular elements, and the quadrilaterals give more accurate results for the same mesh size. If the surfaces are highly warped, curved or swept, you should use triangular elements.”

For ANSYS users the same rule applies in case of using linear plate elements.

Shear locking

In areas where elements are loaded by in plane bending, prevent shear locking by using more than one element over the web height. This is illustrated in the example below. In the area of interest ensure that the elements are as rectangular as possible (preferably square), to give the most accurate results.

The results shown below are calculated with NX Nastran version 8.5 (CQUAD4 elements), but calculations with ANSYS 14 show similar results (SHELL63, SHELL181, BEAM44 and BEAM188 elements). The following example shows an I-beam (h=200mm, tweb=9mm, w=200mm, tflange=15mm) with length L=8m. The end load is W=21.36kN. The beam is clamped at the other side. The theoretical maximum deflection is (according to Euler Bernoulli):

equation
(with E=205000 MPa and I=6611.25·104 mm4)

Five models were run containing five different ways of meshing as illustrated in the figure below.

Beam analysis results with different element shapes

The result of deflection and maximum stresses due to the different element shapes is given in the table below:

Element shape Deflection δ
[mm]
Max stress
[MPa]*
Error in deflection
[%]
Trapezoidal elements 83.1 144.4 69.1
Parallelogram 185.1 209.4 31.2
Rectangle 268.3 259.7 0.26
Square 270.2 279.2 -0.44
Square 3 elements in height 270.2 277.4 -0.44
*Note that these stresses include bending of the flange element.  Therefore this value may be calculated with the theoretical equation σ=W·L·½h/I.  The average value (average max and min of the combined stresses is 259MPa for the square element, which is correct).

To get accurate results with the buckling code checker, when in plane bending occurs, only the last model gives reliable results. The reason is that due to bending the stress differs from the neutral axis to the flange of the beam.  With only 1 element in height, this element stress is not accurately calculated with linear elements in the model.