This preview shows page 1. Sign up to view the full content.
Unformatted text preview: face of the cube that are uniquely determined by values at their four vertices on that
face. Once again, this ensures that shape functions and U are C 0 functions on a uniform
grid of cubes or rectangular parallelepipeds. Since each shape function is the product of
onedimensional linear polynomials, the interpolant is a trilinear function of the form
U(
) = a1 + a2 + a3 + a4 + a5 + a6 + a7 + a8 :
Other approximations and transformations follow their twodimensional counterparts.
For example, triquadratic shape functions on the canonical cube are constructed by
placing 27 nodes at the vertices, midsides, midfaces, and centroid of the element (Figure
4.5.4). The shape function associated with Node (i j k) is given by (4.5.7a) with Ni ( )
given by (4.3.3bd). 4.4. ThreeDimensional Shape Functions 27 4.5.3 Hierarchical Approximations
As with the twodimensional hierarchical approximations described in Section 4.4, we use
Szabo and Babuska's 7] shape function with the representation of Shephard et al. 6].
The basis for a tetrahedral or a canonical cube begins with the vertex functions (4.5.1)
or (4.5.7), respectively. As noted in Section 4.4, higherorder shape functions are written
as products Nik (x y z) = k (
k of an entity function ) i( ) (4.5.8) and a blending function i. The entity function is de ned on a mesh entity (vertex, edge, face, or element) and
varies with the degree k of the approximation. It does not depend on the shapes
of higherdimensional entities.
The blending function distributes the entity function over higherdimensional entities. It depends on the shapes of the higherdimensional entities but not on k.
The entity functions that are used to construct shape functions for cubic and tetrahedral elements follow.
Edge functions for both cubes and tetrahedra are given by (4.4.6c) and (4.4.2e) as
p
Z
k ( ) = 2(2k ; 1)
Pk;1( )d
k2
(4.5.9a)
1; 2
;1
where 2 ;1 1] is a coordinate on the edge. The rst four edge functions are presented
in (4.4.8).
Face functions for squares are given by (4.4.3) divided by the square face blending
function (4.4.3a)
k Here, ( ( ) = P ( )P ( ) + =k;4 k 4: (4.5.9b) ) are canonical coordinates on the face. The rst six square face functions are
400 = 1 510 = 501 = 620 = 3 2;1
2
2;1
611 =
602 = 3
2:
Face functions for triangles are given by (4.4.9) divided the triangular face blending
function (4.4.9a)
k ( 1 2 3) = P ( 2 ; 1)P (2 3 ; 1) + =k;3 k 3: (4.5.9c) 28 Finite Element Approximation As with square faces, ( 1 2 3) form a canonical coordinate system on the face. The
rst six triangular face functions are
300 = 1
401 = 2
511 = ( 3;1 2 ; 1 )(2 3 ; 1) 410 =
520 = 2; 1
3( 2 ; 1)2 ; 1 2
2
5 0 2 = 3(2 3 ; 1) ; 1 :
2 Now, let's turn to the blending functions.
The tetrahedral element blending function for an edge is
ij ( 1 2 3 4 ) = i j (4.5.10a) when the edge is directed from Vertex i to Vertex j . Using either Figure 4.5.2 or Figure
4.5.3 as references, we see that the blending function ensures that the shape function
vani...
View
Full
Document
This document was uploaded on 03/16/2014 for the course CSCI 6860 at Rensselaer Polytechnic Institute.
 Spring '14
 JosephE.Flaherty
 The Land

Click to edit the document details