#include <ElmDefs.h>
Inheritance diagram for ElmT10n3D::
Public Methods | |
| ElmT10n3D () | |
| ~ElmT10n3D () | |
| virtual Ptv(real) | setCentroid () const |
| virtual bool | canBeInside (const Ptv(real) &globpt, const Mat(real) &coor) const |
| virtual bool | isInside (const Ptv(real) &locpt) const |
| virtual bool | isStrictlyInside (const Ptv(real) &locpt) const |
| virtual NumItgDomain | getSideShape (int) const |
| virtual void | geomtFunc (Vec(real) &N, const Ptv(real) &loc_pt) const |
| virtual void | basisFunc (Vec(real) &N, const Ptv(real) &loc_pt) const |
| virtual void | dLocGeomtFunc (Mat(real) &dNloc, const Ptv(real) &loc_pt) const |
| virtual void | dLocBasisFunc (Mat(real) &dNloc, const Ptv(real) &loc_pt) const |
| virtual void | d2LocBasisFunc (ArrayGen(real) &d2Nloc, const Ptv(real) &loc_pt) const |
| virtual void | drawElement (const Mat(real) &coor, GraphBasics &plotter, int resolution=1) const |
| virtual void | drawSide (int side, const Mat(real) &coor, GraphBasics &plotter, int resolution=1) const |
| virtual int | getNoNodesOnEdge () const |
| virtual Ptv(real) | getVertexPt2DGrid (int i, int j) const |
| virtual Ptv(real) | getCenterPt2DGrid (int i, int j, RectanglePartition ind) const |
| virtual Ptv(real) | getVertexPt2DGrid (int, int, int) const |
| virtual Ptv(real) | getCenterPt2DGrid (int, int, RectanglePartition, int) const |
Static Public Methods | |
| void | shape (Vec(real) &N, const Ptv(real) &loc_pt) |
| void | dshape (Mat(real) &dNloc, const Ptv(real) &loc_pt) |
| void | setUpNodeCoor (MatSimple(real) &geomt_coor_) |
| void | setUpNodesOnSides (MatSimple(int) &geomt_nodes_sides_) |
NAME: ElmT10n3D - definition of 3D quadratic tetrahedral finite element
DESCRIPTION:
The class is derived from class "ElmDef" and implements the standard three dimensional 10 noded quadratic tetrahedral element with nodes at the four corners and on the six edges. Compare the relationship between "ElmT3n2D" and "ElmT6n2D".
Most of the operations on this element type are defined through virtual functions or member functions of the base class "ElmDef".
The following numbering is used: the corner nodes, 1, 2, 3 and 4, are as for the "ElmT4n3D" element, node 5 is between node 1 and 2, node 6 is between node 2 and 3, node 7 is between node 3 and 1, node 8 is between node 1 and 4, node 9 is between node 2 and 4, and node 10 is between node 3 and 4.
|
|
The only constructor performs all the initialization. |
|
|
|
|
|
calculates the basis (trial) functions over the element. In other words, the routine computes the interpolation functions for functions defined over the element. Functions entering differential equations etc should be interpolated according to these basis functions. The return array "N" can be empty (not ok) at input, the "basisFunc" will redimension the array properly. The length of "N" equals the number of basis function nodes in the element. For isoparametric elements "basisFunc" and "geomtFunc" are identical. However, for mixed finite elements, one may have a bilinear geometry function (four nodes) and, e.g., a constant basis function. The basis function has then only one node (for example in the centroid of the element) and the length of "N" will be 1. Confer subclasses and the FAQ for information on conventions regarding basis function nodes and the ordering. It should be mentioned that class "BasisFuncGrid" is used to assign basis nodes to a grid (class "GridFE" is meant to contain the geometry of the grid, including geometry nodes). Reimplemented from ElmDef. Reimplemented in ElmT10gn4bn3D, and ElmT10gn1bn3D. |
|
|
returns a true value if a point in the global (physical) coordinate system can be inside the element. The algorithm for detecting whether the point is inside or outside is only approximative. To achieve an exact algorithm, one can use "findLocPt" and then test the local point by using the function "isInside". Reimplemented from ElmDef. |
|
|
the second order derivatives in the local element coordinates of the functions defined by "basisFunc". Reimplemented from ElmDef. Reimplemented in ElmT10gn4bn3D, and ElmT10gn1bn3D. |
|
|
the derivatives in the local element coordinates of the functions defined by "basisFunc." As for the "geomtFunc" and "basisFunc" functions the return array can be empty or have wrong dimension at input. The function will redimension the array properly. Reimplemented from ElmDef. Reimplemented in ElmT10gn4bn3D, and ElmT10gn1bn3D. |
|
|
the derivatives in the local element coordinates of the functions defined by "geomtFunc." As for the "geomtFunc" and "basisFunc" functions the return array can be empty or have wrong dimension at input. The function will redimension the array properly. Reimplemented from ElmDef. |
|
|
function for drawing an element in the global, physical coordinate system. The drawing consists of geometric primitives (line segments) in the "GraphBasics" format. Note that this function does not flush or close the "GraphBasics" file so the physical file may be incomplete after a call to this function. An explicit call to "GraphBasics close" will complete the file and make it ready for reading. The function is used in drawing routines in class "DrawFE". Reimplemented from ElmDef. |
|
|
function for drawing a specified side of an element in the global, physical coordinate system. See "drawElement" for comments. Reimplemented from ElmDef. |
|
|
implements the first derivatives of the basis functions for this element type. Primarily for internal use. |
|
|
calculates the geometry functions on the element (see the description above). These functions are used for the isoparametric mapping associated with the element (these functions determine the geometric shape of the element in the global, physical coordinate system). Recall that the geometry functions are defined with respect to the local element coordinate system and the local node numbering of the element. The return array "N" can be empty (not ok) at input; the "geomtFunc" will redimension the array properly. The "N" array contains the values of the geometry functions at the evaluation point. The length of the "N" array equals the number of (geometry) nodes in the element, and the ordering of the entries in "N" follow the ordering of the local nodes Reimplemented from ElmDef. |
|
|
See documentation of one of the overloaded functions. Reimplemented from ElmDef. |
|
|
as "getVertexPt2DGrid" but the center point of the box or triangle with index (i,j) is returned. The lower, leftmost box has index (1,1). If triangles are desired, the user must indicate whether the center of the lower or the upper triangle is desired. The "enum" "RectanglePartition" takes the names "LOWER", "UPPER" or "RECTANGLE" indicating lower triangle, upper triangle or four-sided box, respectively. Reimplemented from ElmDef. |
|
|
returns the number of nodes on an edge (an edge is the intersection between two adjacent sides in a 3D element). The function has no meaning for 1D and 2D elements. Reimplemented from ElmDef. |
|
|
returns the shape of the sides of the element. This is the counterpart to "getDomainShape" for example for numerical integration over the sides. Reimplemented from ElmDef. |
|
|
See documentation of one of the overloaded functions. Reimplemented from ElmDef. |
|
|
returns a "Ptv(real)" containing the coordinates of a vertex point (i,j) in the 2D grid over a 2D element or over a side in a 3D element (see function "setResolution2DGrid)." A vertex point is a corner of the boxes/triangles of the 2D grid. The lower, leftmost vertex has coordinates (1,1). Reimplemented from ElmDef. |
|
|
returns a true value if the local coordinates of a point implies that the point is inside the element. The function is mainly used to check if "findLocPt" has found a point that is really inside the element. Reimplemented from ElmDef. |
|
|
returns a true value if the local coordinates of a point implies that the point is -strictly- inside the element, i.e. inside but not on the boundary. Reimplemented from ElmDef. |
|
|
Reimplemented from ElmDef. |
|
|
load the "MatSimple" object with coordinate values for each node in the element. The row index refers to local node numbering, and the column index refers to the spatial direction. |
|
|
similar to the "setUpNodeCoor", this function fills a "MatSimple" object with node numbers for each side. The row index refers to side numbers, and the column index refers to local node numbers. |
|
|
implements the basis functions for this element type. Primarily for internal use. |