It is by now well established that the computational analysis of significant problems in structural and continuum mechanics by the matrix displacement method often requires elements of higher sophistication than used in the past. This refers, in particular, to regions of steep stress gradients, which are frequently associated with marked changes in geometry, involving rapid variations of the radius of curvature. The philosophy underlying the idealisation of such configurations into finite elements was discussed in broad terms in ref. 1. It was emphasised that the so successful, constant strain, two-dimensional TRIM 3 and three-dimensional TET 4 elements do not, in general, prove the best choice. For this reason elements with a linear variation of strain like TRIM 6 and TET 10 were originally evolved and followed up with the quadratic strain elements TRIM 15, TRIA 4 (two-dimensional) and TET 20, TEA 8 (three-dimensional) of ref. 2. However, all these elements are characterised by straight edges and necessitate a polygonisation or polyhedrisation in the idealisation process. This may not be critical in many problems, but is sometimes of doubtful validity in the immediate neighbourhood of a curved boundary, where stress concentrations are most pronounced. To overcome this difficulty with a significant (local) increase of elements does not always yield the most economical and technically satisfactory solution. Moreover, there arises another inevitable shortcoming when dealing with TRIM and TET elements with a linear or quadratic variation of strain. Indeed, while TRIM 3 and TET 4 elements permit a very elegant extension into the realm of large displacements, this is not possible for the higher order TRIM and TET elements. This is simply due to the fact that TRIM 3 and TET 4 elements, by virtue of their specification, always remain straight under any magnitude of strain, but this is not so for the triangular and tetrahedron elements of higher sophistication.
Programming for solving plane rigid frame based on MATLAB
