FEM_BASIS_T6_DISPLAY Display a T6 Finite Element Basis Function

FEM_BASIS_T6_DISPLAY is a MATLAB program which reads data defining a quadratic triangular ("T6") mesh and displays the basis function associated with any node in the mesh.

The entire mesh is shown (which could be a problem for large meshes!). However, a finite element basis function is nonzero over only a relatively small number of elements. Over these elements, the basis function is quadratic, and the curved nature of this part of the surface is suggested by taking many plot points in this region.

The best part of the presentation is that MATLAB allows the user to turn the plot or zoom in and out, which makes it much easier to understand the information that is being presented.

The triangular mesh is defined by a node file containing the coordinates of nodes, and an element file containing sets of the six nodes that make up each element.

Note that the six nodes making up each element are to be listed in a particular order, in which the vertices come first, followed by the midside nodes. Moreover, the first midside node is the one between vertices 1 and 2, as suggested by the following diagram:

```        3
|\
| \
|  \
6   5
|    \
|     \
|      \
1---4---2
```

In general, a mesh of quadratic triangles allows the sides to be bent. However, for this application, it is assumed that the sides of each triangle are straight, that is, that each midside node lies on the line segment joining the corresponding vertices.

Usage:

fem_basis_t6_display ( 'prefix' )
where
• 'prefix' is the common filename prefix of the node file 'prefix'_nodes.txt and the element file 'prefix'_elements.txt
The program queries the user for a node index and then displays the basis function corresponding to that node.

Languages:

FEM_BASIS_T6_DISPLAY is available in a MATLAB version.

Related Programs:

BEZIER_SURFACE_DISPLAY, a MATLAB program which displays a Bezier surface.

FEM_BASIS_T3_DISPLAY, a MATLAB program which displays a basis function associated with a 3-node triangle "T3" mesh.

FEM_BASIS_T4_DISPLAY, a MATLAB program which displays a basis function associated with a 4-node triangle "T4" mesh.

FEM2D_PACK, a MATLAB library which is useful for 2D finite element calculations.

POLYGONAL_SURFACE_DISPLAY, a MATLAB program which reads two files defining a polygonal surface and displays it within MATLAB.

QUAD_SURFACE_DISPLAY, a MATLAB program which reads files defining a 3D quadrilateral mesh surface and displays it within MATLAB.

TET_MESH_DISPLAY, a MATLAB program which can read in the node and tetra files defining a tet mesh and display a wireframe image.

TRI_SURFACE_DISPLAY, a MATLAB program which reads data defining a triangular mesh of a 3D surface and displays it.

TRIANGULATION_DISPLAY_OPENGL, a C++ program which reads files defining a triangulation and displays an image using Open GL.

TRIANGULATION_ORDER6_CONTOUR, a MATLAB program which can make contours of a scalar quantity defined on the nodes of an order 6 triangulation.

TRIANGULATION_PLOT, a MATLAB program which creates an EPS image of a triangulation.

Reference:

1. Klaus-Juergen Bathe,
Finite Element Procedures,
Prentice Hall, 1996.
2. Eric Becker, Graham Carey, John Oden,
Finite Elements, An Introduction, Volume I,
Prentice-Hall, 1981.
3. Herbert Edelsbrunner,
Geometry and Topology for Mesh Generation,
Cambridge, 2001,
LC: QA377.E36,
ISBN: 0-521-79309-2.
4. Max Gunzburger,
Finite Element Methods for Viscous Incompressible Flows,
A Guide to Theory, Practice, and Algorithms,
LC: TA357.G86.
5. Gilbert Strang, George Fix,
An Analysis of the Finite Element Method,
Prentice Hall, 1973,
LC: TA335.S77.

Tests and Examples

SMALL_SQUARE is a simple mesh of 25 nodes and 8 elements over the [0,4] by [0,4] square. Test files you may copy include:

BIG_SQUARE is a simple mesh of 169 nodes and 72 elements over the [0,1] by [0,1] square. Test files you may copy include:

You can go up one level to the MATLAB source codes.

Last revised on 12 February 2010.