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.
Model solved; simply-supported condition; linear analysis
Typical finite element mesh; 4-node incompatible mode elements; one element through thickness and N equal elements along length; in the solutions using 9-node elements, each 4-node element is replaced by a 9-node element
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.
Results for simply-supported condition, linear analysis
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.
Typical finite element mesh with fully built-in support condition
Results for fully built-in support condition, linear analysis
Of course, point-load conditions, as seen in the
- K. J. Bathe, Finite Element Procedures, Prentice Hall, 1996.
- D. Pantuso and K. J. Bathe, “On the Stability of Mixed Finite Elements in Large Strain Analysis of Incompressible Solids”, Finite Elements in Analysis and Design, 28, 83-104, 1997.