Remarks on the Use of Incompatible Modes
In the analysis of solids, the 4-node plane stress/strain and axisymmetric elements
The incompatible mode elements can be thought of — and can be formulated — as
The incompatible mode elements are well known to be most powerful when undistorted, e.g. when rectangular in shape in 2D analyses.
However, it is known that the use of enhanced strain elements can provide difficulties in large strain (almost incompressible) analysis, see e.g. ref. , and the same holds for elements based on incompatible modes. For such analyses, the displacement/pressure (or u/p) elements available in ADINA are more reliable and effective.
It seems less well known that the elements based on incompatible modes can show what may be considered a surprising behavior in linear analysis, as illustrated in the following example.
A simply-supported axisymmetric plate is subjected to pressure loading. A schematic of the model considered is given above. We also show a typical finite element discretization.
The use of
On the other hand, the incompatible mode element results are only acceptable as long as the elements are not too thin. Indeed, if only results up to N = 100 are considered, the solution has converged. In practice, surely, no more elements will be used.
However, merely as a convergence study, if we increase the number of elements further, an increasing displacement at the center of the plate is seen. The elements near the support cause the displacement to increase unphysically, due to the incompatibility between the elements. Of course, the increase only occurs when the elements are very thin and then the increase is also small, but it is good to know about this phenomenon.
Repeating the analysis with the fully built-in support condition gives the results in the following figure. The divergent behavior using the incompatible mode elements is not seen.
Of course, point-load conditions, as seen in the