What is shear locking?
Shear locking is an error that occurs in finite element analysis due to the linear nature of quadrilateral elements. The linear elements do not accurately model the curvature present in the actual material under bending, and a shear stress is introduced. The additional shear stress in the element (which does not occur in the actual beam) causes the element to reach equilibrium with smaller displacements, i.e., it makes the element appear to be stiffer than it actually is and gives bending displacements smaller than they should be.
How to prevent shear locking?
In areas where linear elements are loaded by in plane bending, shear locking is prevented by using preferably elemenets over the 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 Ibeam (h=200mm, t_{web}=9mm, w=200mm, t_{flange}=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):
(with E=205000 MPa and I=6611.25·10^{4} mm^{4})
Five models were run containing five different ways of meshing as illustrated in the figure below.
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.
What else is important besides preventing shear locking?
When modeling a model using linear 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.
